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The properties of equilibrium boundary layers 


3y A. A. TOWNSEND 


Emmanuel College, Cambridge 
(Received 7 June 1956) 


SUMMARY 

The recent work on equilibrium (i.e. self-preserving) turbulent 
boundary layers in adverse pressure gradients is compared with 
theoretical predictions based on current generalizations about the 
turbulent shear flow. Using only assumptions of similarity, it is 
possible to show that an equilibrium layer can exist only if the free 
stream velocity varies as a power of distance downstream with an 
exponent greater than — } and if the velocity defect from the free 
stream is small. Assuming further that the effective eddy viscosity 
is independent of distance from the wall over the outer part of the 
layer, most of the properties of equilibrium layers may be computed 
from the known behaviour of layers in zero pressure gradient. ‘The 
predicted values of skin friction and the predicted shape and 
magnitude of the mean velocity distribution are in fair agreement 
with the observations of Clauser. Finally, the modifications that 
are necessary if the velocity defect is not small are discussed briefly. 


1. INTRODUCTION 

Much of the recently published work on turbulent boundary layers has 
abandoned the long-standing tradition of purely empirical approach based 
on direct analogy with the behaviour of laminar layers and has begun to use 
generalizations about the characteristics of turbulent shear flow. ‘lhe 
movement began, perhaps, with the work of Ludweig & ‘lillmann (1949) 
who established experimentally that the distribution of mean velocity near 
a smooth wall is independent of the pressure gradient in the free stream, and 
is determined by the shear stress at the wall and the fluid viscosity. ‘This 
conclusion, which has since received abundant experimental confirmation 
(e.g. Coles 1956), is more than an empirical generalization, being a necessary 
consequence of the mixing-length and similarity theories of turbulence and 
indeed of any theory of turbulence that assumes net turbulent transfer of 
energy to be small within the wall region of substantially constant shear 
stress. At much the same time it was being realized that the apparent 
variation in shape of the velocity profiles for a constant pressure layer is an 
illusion, produced by measuring velocities relative to the stationary wall 
rather than to the free stream. Profiles of velocity defect from the free 
stream velocity are accurately similar in shape except within the viscous 
layer next the wall (e.g. Clauser 1954). Again, this ‘defect law’ is not 
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only an empirical result but is a consequence of two widely accepted generali- 
zations about turbulent flows, that the general motion of a fully developed 
turbulent flow is independent of the fluid viscosity*, and that developing 
flows tend to a condition of moving equilibrium or self-preservation. 

Although these two results are in a sense experimental rediscoveries of 
much earlier theoretical work (see, for example v. Karman (1930) and 
Goldstein (1938)), they are none the less valuable, and their application 
to boundary layers in pressure gradients is of great interest. A funda- 
mental advance in this direction was made by Clauser (1954) who, by trial 
and error adjustment of the external pressure gradient, constructed 
‘equilibrium’ boundary layers satisfying the defect law and who discussed 
and analysed their behaviour. It is possible, using only the generalizations 
employed in inferring the existence of a universal velocity profile near a wall 
and the defect law for constant pressure layers, to determine the conditions 
necessary for the existence of an equilibrium boundary layer in a pressure 
gradient (Townsend 1956). The purpose of this paper is to summarize 
the results of this analysis, including the consequences of assuming an 
effective eddy viscosity independent of distance from the wall in the outer 
layer, and then to compare the predictions of this theory with the results and 
generalizations of Clauser (1956). 


2. NOTATION 
Consider a boundary layer on a smooth flat plate whose surface is the 
plane, y = 0. Oz is the direction of homogeneity at right angles to the two- 
dimensional mean flow, and Ox is parallel to the surface and in the general 
direction of the flow. ‘Then, 


U,V are the components of the mean velocity parallel to 
Ox and Oy respectively, 
u,v, w are the components of the turbulent velocity 
fluctuation parallel to Ox, Oy, Oz, 
P is the pressure, 
v is the kinematic viscosity, 
u. &, are the Ox component of the mean velocity and the 


pressure in the free stream just outside the layer, 
T is the shear stress at the wall, 
K, A are the constants in the universal logarithmic 
distribution of mean velocity near a smooth wall 
(equation 3.12), 


Cy, = 275/U; is the local coefficient of skin friction, 
y = Th*/KU, is a non-dimensional quantity specifying skin 
friction, 
uy = Th*/K is the velocity scale of the boundary layer, 


* This statement is not inconsistent with the dependence of the velocity profile 
near the wall on fluid viscosity. 'The magnitude of the viscosity merely determines a 
velocity of translation for the fully turbulent part of the flow, and velocity differences 
within it are unaffected. 
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§ = vir¥exp(- -A)is the length scale (or ‘thickness’) of the layer, 





7 denotes y/6, 
f(y) = (U—U,)/uy is a universal function specifying the mean velocity 
distribution, 
215 Lo» £19 are functions specifying the distributions of tur- 
bulent intensities and stresses, 
= —)J finan, Te= | (fC)? @n, 
# ( ~ 0 
‘ sale i 1 . | I, 
} = K*hytexp( —4), W.-Y] 14 (247 ah, 
Y \ I, J 
“U,dx U,x/v. 
R. = ; = . is the Reynolds number, 
Jov : 
Vp is the effective eddy viscosity in the outer layer, 
R, = uy 8/vp is a non-dimensional constant relating the effective 
eddy viscosity to the velocity distribution, 
s is a non-dimensional variable of position, 
6 is a constant determining the scale of the velocity 
defect. 


It should be noted that the pressures and stresses are ‘ kinematic’, that is, 

they are the ordinary mechanical values divided by the fluid density. Also, 
since the pressure in the free stream is related to the velocity by 

ou, oP, 

Uu,~— =-=> 


Ox Ox 


(2.1) 





’ 


its gradient is specified by the variation of free stream velocity. For the 
most part, we consider a free stream velocity variation 


U, =a Be (« >, (2.2) 


which represents an adverse pressure gradient if the exponent is negative. 


3. THE CONDITIONS FOR SELF-PRESERVING DEVELOPMENT OF A 
BOUNDARY LAYER 
The defect law for the mean velocity distribution in a boundary layer is 
equivalent to the assumption of self-preserving profiles in a wake, i.e. 

U = U,+u f(y/), (3.1) 
where u, and 6 are scales of velocity and length depending on x, and the 
function is independent of x. For the establishment of a self-preserving 
flow it is necessary that the distributions of shear stress and turbulent 
intensities should be expressible in terms of the same scales by other functions 
independent of distance downstream, say 


we = usgey/5), uw = uhgi(y/d), = uh galy/S). (3.2) 
Naturally, the distributions must satisfy the equations of mean motion, 
which, for the boundary layer, and to the approximation usually found 
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sufficiently accurate, can be reduced to the single equation, 





ye _ aU , Guv Auto), OU, BU ae 
ee ee eee 7. ae es ai 1 ae (3.3) 
Ox oy cy ON x oy” 


In terms of the self-preserving functions of equations (3.1), and (3.2) this 
may be written 
d(u,U,) . uy d(U,5) ,, Bile , ra | 
ae f_ —— nf + Uo ay (f? + 2(81-82))— 
uy (uy) ., { - updd on, us _ Vly ce : 
8 ax J | J dn — 5 Te 81 —82)+ "§ O12 = xl ? (3.4) 





where dashes indicate differentiation with respect to 7. 

This equation may be satisfied exactly by self-preserving functions 
independent of x if the ratios of the non-zero coefficients are independent 
of x. For laminar flow, the functions representing the fluctuating motion 
are identically zero and self-preserving flow is possible either if 

ty = Uy x, 8 oc xft-ai2, (3.5) 
or if 
ua=U, «ce, bac, (2c > 0), (3.6) 
(Goldstein 1939), If 6 is now defined by 
—-—=— = —— =f say 
v dx l+a si 
substitution in equation (3.4) gives an equation for f(»), 


anf (x : lf a) +nf? = f”. (3.7) 


~ 


(The exponential flow (3.6) corresponds to nm = 2.) When the non- 
dimensional stream-function defined by 
o'(n) = 14+f(n), 


is introduced, this becomes the Hartree form of the Falkner-Skan equation, 


ue tas” = n(us’* 1). (3.8) 
It is known that acceptable solutions satisfying the boundary conditions, 
b(0) = £'(0) = 0, (0) = —n, Y'( 2%) = 1, exist only if n —(-199 or 
a -()-11 (Hartree 1937). ‘This sets a lower limit to the possible values of 


the exponent and to the severity of the adverse pressure gradient. 
In turbulent flow, a self-preserving flow is possible only if 
uy « U, x (x,—x)}, 0 © (X)—x), (Xe 2), (3.9) 
corresponding to the accelerated boundary layer within a wedge. ‘The 
equation may be satisfied approximately if the velocity defect is small, that 


is, if 
Uy f < U;, (3.10) 
unless 7» < 5. ‘Then equation (3.4) may be approximated by 
dull.) teGU,d) ., 8, as 
abso “0 at C ek " f gt () (3.11) 


dx ° > ae >= 
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except very close to the wall where direct viscous stresses are appreciable. 

The existence of a wall layer within which the velocity distribution is given 

by 2 1/2 
i . TO TY 

Us =| log—~ +A (3.12 

ri = * (3.12) 


implies that we may choose as scales 


K 
and s= =r exP( as! -4). (3.13) 


It may then be shown that self-preserving flow is only possible for free 
stream velocity variations of three types, 

(2) U,=Bx* (x>0,a> —}), 

(5) U, = B(-x)* («x >0,a < -}), 

(c) U, = Be (p> 0). 
Of these only the first can represent retarded flow (for 0 > a > — 4) and 
correspond with an equilibrium layer of the type studied by Clauser. For 
this type of velocity variation, the relation between wall stress and position 


is given by aoe 
Y = K-J, yr texp(- -A) = 2 


a+l1 »v 
with the approximation (3.10). Here y = 7//*/KU, is a friction parameter, 


Il 





(3.14) 


andJ, = —| f(») dy isa constant of the flow anda function of a. It should 
~ 0 
be emphasized that the wall stress is uniquely related to the velocity in the 
free stream and so to the position. ‘The existence of the proper variation 
of free steam velocity (or pressure) alone is not sufficient to ensure self- 
preserving development; the boundary layer must be matched to the 
pressure gradient. 
Clauser has observed that the parameter 


5* d. “© 
Il = 6* dP, where 6* = | (1- w) ay 











To dx 7 JO U, 
is independent of x in an equilibrium layer. In the present notation, 
I, uj d dU a 2 
n= -+3 >= - . (3.15) 
Kuz dx 1+3a 


using equations (3.13) and (3.14). The constancy of this parameter is a 
simple consequence of this very general theory of boundary layers with a 
small velocity defect. It is also the condition that the terms in the equation 
for the momentum integral should preserve a constant ratio, 


4. THE VELOCITY DISTRIBUTION IN THE OUTER PART OF AN EQUILIBRIUM 
LAYER 

Experience with other self-preserving free turbulent flow suggests that 

an acceptable description of the mean velocity distribution may be obtained 

by assuming an effective eddy viscosity constant over the outer part of the 








566 A. A. Townsend 


layer, the mean velocity distribution due to this being joined smoothly to 
the logarithmic distribution that obtains in the inner layer (Clauser 1956, 
Townsend 1956). Clauser has done this without assuming the velocity 
defect to be small. His results can be derived assuming the relation (3.14) 
between wall stress and position, although this relation is not accurate for 
finite velocity defect. In the outer layer, the terms of equation (3.4) that 
involve the turbulent intensities are usually negligible, and then, putting 
= ou 


—yv = Vr : (4.1) 


a 
where uy d/vp = R,, a constant depending on a, the equation of mean motion 
becomes 


ee _ l+a i a : 
” a ; ae _— ae 
RR! 1+ 3at 1+ 3al E | Y | J an | 1+3a Yh 0. (4.2) 


: a , K?R, 1+a\}? — 
In terms of a new position variable s = 7, this is 
I, 1+3a 


d*f nf + 7 (s - y| Fs ds) — Inyf? = 0. (4.3) 











ds* ds 
The non-dimensional stream function defined by 


Ui(s) = 1+ yfls) 
Ui" + yp" = An(y’2— 1). (4.4) 


This is the Hartree equation with two differences. The first and obvious 
one is that the coefficient outside the bracket is half its value in the equation 
for laminar flow (3.8), and the second is that the boundary condition at s = 0 
is no longer (0) = 0 but ¥’(0) = 1+yf(0). This boundary condition is 
not independent of x, and does not conform to the initial assumption of self- 
preservation. 

For very small velocity defect in the outer layer, the equation for f(s) 


becomes 
af df .— ; 
ds2 Bi 4 —nf = (), (4.5) 


satisfies the equation 


The appropriate solution of this equation is a multiple of the function 

Hh,(s), defined by 

x t” : 
Hh,(s) = | = exp[—4(t+s)?] dt. (4.6) 
Jo n! , 
The complete variation of mean velocity is obtained by adjusting the 
multiplying constant of the function describing the variation in the outer 
layer to give a smooth junction with the velocity distribution in the inner 
layer where 

f(y) = log». (+.7) 
At high Reynolds numbers, part of the region of constant stress is within 
the region of small velocity defect, so that the mean velocity distribution 
within this region of overlap must be both self-preserving and identical 
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with the universal logarithmic profile. The conditions for this have been 
given as equations (3.13), and this process of matching the velocity profile 
for constant eddy viscosity to the logarithmic profile is no more than a 
working approximation consistent with these conditions and the assumptions 
of self-preservation. 


5. THE HYPOTHESIS OF LARGE EDDY EQUILIBRIUM 

There is a considerable body of evidence showing that the rate of spread 
of a turbulent flow into the ambient undisturbed fluid depends on the 
presence of a system of large eddies which are also responsible for the 
observed intermittency of turbulent flow near the free stream boundary. 
‘The author has postulated that these large eddies are an integral part of the 
mechanism that regulates the intensity of the turbulent motion and the rate 
of entrainment of non-turbulent fluid (Townsend 1951, 1956). By assum- 
ing a plausible form for these eddies, it is possible to obtain a complete 
solution for the boundary layer problem in terms of the universal constants, 
K and A, of equation (3.12) and of two numbers specifying the position and 
scale of the large eddies with respect to the mean velocity distribution. 
Estimates of these numbers may be made from a consideration of the origin 
and dynamics of the large eddies, or they may be determined by analysing 
experimental measurements of the layers. If it is assumed that these 
numbers are independent of pressure gradient, all the properties of equi- 
librium boundary layers can be predicted from the known properties of 
constant-pressure layers without ambiguity. For details, reference should 
be made to the full account (Townsend 1956). 

To these approximations, the mean velocity distribution consists of the 
distribution 


f(n) = logy (5.1) 
in the inner layer (say for 7 < yp) and the distribution 
f(y) = —CHh,(s), (5.2) 
where 
K?R, \'2 Ud 
‘= (ra =) * oD 


in the outer layer (y > ). As Clauser points out, the experimental deter- 
mination of mean velocity profiles is not sufficiently accurate to make profit- 
able a comparison of the shapes of observed and predicted profiles, and the 
test of a theory must be the prediction of the scales of variation. 

If the velocity distribution in the outer layer is extrapolated to the wall, 
the extrapolated velocity there, U, should be given by 


7 U,- U, 


K = CHh,(0) = c VG) (5.3) 


Fe | jn) ! . 
The scale of velocity variation is so defined by the constant C, which is given 


by a 
’ i | 22, ~, (l+n 22+" (/1+n\ 2 7 
i ~« [ 1-7 (> i- (+ yy], (5.4) 


lin| v7 rr 


1/2 
0 








= a ra 
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the constants, 22 and 1.7, being selected to conform to observations of 
constant-pressure layers. If (1+) is small, 


6°66 = 
, —— So) 
. 1 att ‘ ( 
Clauser uses two scales of distance from the wall. The first is 
oU,-—U [ sd 
A= | ay = Bo (5.6) 


Distances expressed in this scale are related to distances expressed in the 
scale of the variable s (equation (4.3)) by 


% = (Fe): = [k(1+n)]*s. ( 


The second scale is 59, the total thickness of the layer, presumably the 
distance from the wall beyond which the mean velocity variation is im- 
perceptible. ‘This scale will be roughly three times the scale of s. 
Clauser expresses the eddy viscosity in the outer layer in terms of the 
non-dimensional constant, k, which is given by 
1 U,5* | ie mwik* 12° 
k V7 hii 22 1-7 (—")1- le (): \? (5.8) 


V7 


mn 


.7) 

















The comparison of these predictions with Clauser’s measurements is 
shown in figures 1 and 2, using the values of m in table 1 computed from 
equation (3.15). This procedure gives a definite value for a, and should 
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| 0-318 | 
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0-20 
00-0667 | 


Distribution | 
Distribution II | 
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Table 1 


make some allowance for the effect of a finite velocity defect. ‘hese values 
of a are in fair agreement with the observed distributions of free stream 
velocity (figure 3). 

The computed variation of C agrees fairly well with observation, the 
difference being most marked for the lesser pressure gradient. The theory 
also predicts a nearly constant value for k in agreement with Clauser’s 
findings, although the computed values are 20-25%, less than the observed 
values. This is caused by the use of the velocity profile (5.2) to compute J,. 
By overestimating the velocity defect near the free stream, this leads to 
values of J, about 10° greater than the true values. The constancy of k 


r OO 


is not unexpected as U, 8* = | (U,—U) dy is nearly the simplest possible 


expression for the product of length and velocity scales for distributions of 
nearly similar shapes. 
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Figure 3. Comparison of Clauser’s pressure distributions with the appropriate 
power laws. 


Comparison of predicted and observed scales of distance may be made 
by computing the position at which the velocity defect is one-half of the 
extrapolated velocity defect, U,— Up, or the position at which the velocity 
defect is small. Some results of this comparison are shown in table 2. 

Lastly, the skin friction may be computed from the friction equation 

Y = K-, ytexp( = -A) - _—_ a ll (5.9) 


Vv 





using values of J, given by 


I, = ((1-+n)C*+ Texp| 1 | (5.10) 


"™(In) | 
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If 1+7 is small, this is nearly 


ee 6°66 \12 e 
I, = 7-66 exp| 1- (7) / (5.11) 


54/2 So 
A do Vi /2 do ee see Hh, (so) 
Obs. | Calc. | Obs. 





3-06} 0-010 


-26 1-18 
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Figure 4+. Observed and theoretical distributions of friction coefficient. 


Figure 4 compares the predicted values of local friction coefficient with the 
experimental measurements. ‘The agreement is good only for the larger 
pressure gradient. 
6. THE CRITICAL VALUE OF THE POWER LAW EXPONENT 
The theory for small velocity defect in the outer layer, which is valid at 
extremely high Reynolds numbers (equation (3.14) shows y to decrease in- 
definitely with increasing Reynolds number), predicts that equilibrium 
layers exist only if a > —}4. ‘The basis of this prediction is the necessity 
for a positive value of the friction function Y. If the effective eddy viscosity 
is constant within the outer layer, equation (4.6) shows that self-consistent 
velocity profiles exist, the critical profile for a = — } being 
f(s) = —Ce™*, (6.1) 


which has a zero gradient at the wall indicating zero wall stress. 
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‘There are two simple ways of looking at the consequences of finite 
velocity defect, viz. by including second-order terms in the derivation of the 
stress-position relation (assuming self-preservation of velocity profiles and 
wall stress) or by finding the range of a for which the modified Hartree 
equation (4.4) has acceptable solutions (assuming that the coefficients of 
equation (3.4) may be computed from the relation (3.14)). Ideally, both 
should be investigated but the result is a modified form of equation (4.4) 
whose solutions have not been studied. 

The equation for the momentum integral is 


1 . ky 
aR | th — (21, + I,)y}y *exp(- a 4) - 


v dU 1 
Saeed ae we aes ee 6.2 
t OR de (21, -1, v)y exp(- 4) K%, (6.2) 


where terms of order less than y have been omitted. With the use of friction 
function 


1 
Y, = KL, -(21,+ 1,)yly *exp(- -A) , (6.3) 
it becomes 
ay, 2a 7 
G2 IO en (ay gee ar) (ee ee ) 
dR, Il1+a [tate t tel R, , 6-4) 


to the same approximation. Since y is a slowly varying quantity, an 
approximate solution may be obtained by neglecting its variation. This is 
l+a 


Y= Rely 3a4 2a tl)’ my 





which has meaning only if 
1 
——-342ay(24+40,/h)° 
The quantity y/,//, is nearly proportional to the velocity defect ratio 
(U,—U,)/U,. For the critical velocity profile defined by equation (6.1), 
RC. 1 U,-U, 
"a 


a 





(6.6) 





(6.7) 


Alternatively, we may suppose the development of the layer to be des- 
cribed by the theory for small velocity defect, and look at the conditions for 
acceptable solutions of the modified Hartree equation (4.4) for various 
defect ratios. Numerical solutions quoted by Clauser (1956) show that 
the critical value of the exponent varies with defect ratio in the way shown 
in figure 5. A comparison with critical values obtained from equations 
(6.6) and (6.7) shows fair agreement between the two estimates only for 
small velocity defect ratios. 

Neither of these arguments is altogether satisfactory and further work 
is desirable. 
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Figure 5. Critical values of the power-law exponent. 


7. CONCLUDING REMARKS 


vi 
“I 
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The properties of self-preserving or equilibrium boundary layers may 


well prove to be of considerable importance in the understandi 


ng of the 


general behaviour of boundary layers in pressure gradients. ‘he theory 
based on the assumptions of small velocity defect and large eddy equilibrium 


fits the observations of Clauser as well as can be expected in vi 


ew of the 


comparative crudity of the theory and the absence of disposable constants. 
It should be noted that the hypothesis of large eddy equilibrium is very 
nearly equivalent to assuming eddy viscosity to be a universal multiple of 


the integral of the velocity defect. 
The condition that the velocity defect is small is satisfied very 


weakly in 


the experimenta: layers, and better agreement between theory and experi- 


ment might be obtained by considering the effects of finite defec 
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The generation of edge waves by moving pressure 
distributions 


By H. P. GREENSPAN 


Harvard University, Cambridge, Massachusetts 


(Received 1 August 1956) 


SUMMARY 

An analytical study is made of the resurgent wave motion 
induced by pressure distributions moving parallel to a straight 
coast line. ‘The resurgence is shown to consist of an infinite 
number of edge wave modes; expressions for these modes are 
given and the wavelengths, frequencies and amplitudes are shown 
to be in agreement with experimental results. The effects of a 
Gaussian pressure distribution are analysed. For large-scale 
disturbances off the east coast of the United States only the funda- 
mental mode is excited. Conditions, such as storm size, speed 
and distance from shore, for maximum induced wave amplitudes 
are derived. 


1. INTRODUCTION 

It has been observed that hurricanes and other similar phenomena 
travelling approximately parallel to a neighbouring coast line sometimes 
induce a resurgent wave motion. ‘These waves are characterized by the 
facts that they propagate along the shore line, are essentially confined within 
a distance of one wavelength from the coast (hence the name ‘ edge waves’), 
have wavelengths measured in miles, amplitudes measured in feet, and 
periods measured in hours. 

Munk, Snodgrass & Carrier (1956) conjectured that this resurgence 
was primarily due to the deviation of the storm pressure distribution from 
normal atmospheric pressure, and was neither wind-generated nor a 
consequence of other effects. Using a simple mathematical model for the 
pressure deviation, with the assumption of quasi-steady state conditions 
(that the disturbance has been in existence and in motion for an infinite 
time interval), they obtained results in agreement with experimental evidence. 
Computed periods, and amplitudes agreed with known values. Their 
original conjecture concerning the cause of the resurgence was thereby 
substantiated. 

In this paper, the transient problem (i.e. with the disturbance originating 
at a finite time) is considered and applied to a more precise physical model. 


2. ‘THE FUNDAMENTAL EQUATION AND ITS SOLUTION 


In order to motivate the choice of an appropriate model on which to 
base our calculations, we anticipate (as indicated in the introduction) that 
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the hurricane generated surface waves will display an amplitude 7 of the 
order of 3 feet, a wavelength A of about 200 miles, and a period of about 
6 hours. The extremely long wavelengths and small wave heights suggest 
that the linear shallow water theory may be applied. In order that this 
theory may provide a valid approximation, the wavelength of a progressive 
wave propagation must be large compared to the ocean depth, the ocean 
depth must be large compared to the wave height, and the gradients of the 
wave height and depth must be small compared to unity. Simple estimates 
of the orders of magnitude involved show that in the phenomena in question 
these conditions are satisfied. 

The linearized theory essentially eliminates the depth variable z from the 
exact analysis by averaging the velocity components over the depth. The 
pressure is then hydrostatic (see Lamb 1945 for a more complete analysis). 








Figure 1. Fluid with a fixed boundary and a free surface. 


With the meaning of the symbols shown in figure 1, and with uniform 
density p, the conservation equations of mass and momentum are 








on(x, y, t) ee ts e(uh) # e(vh) F (1) 

ot Ox dy ” 
eu(sy) _ _ 1 aP(nyat) oe _ _ 1LaPlunat) — 
at _ w ° “es a. -—- ‘ 


where P(x, y, 2,t) = P(x,y,0,t)+gp(y—2), and P(x,y,0,t) is the applied 
surface pressure. We may eliminate u and v to obtain 
1 d2y 1 ‘ 
V .(AVn) 7 oe po (AV?P+Vh.VP). (3) 
Here and in the following, P(x, y,0,t) is replaced by P or P(x,y,t). We 
can also define a potential function by the relations u = dd/dx, v = dd/ey, 
and eliminate 7 to obtain an equation for ¢: 


16 1 @P 


V.AV8)— 5 Se = oe ap (4) 
the wave height being 
P 10¢ : 
ee eee ee a 
b pe got ©) 


It remains to specify the ocean depth h, and to linearize the topography. 
We consider only a straight coast line, and take the ocean depth to be a 








ui 
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linear function of distance from shore y, i.e. h = ay. For the constant «, 
we choose the average slope of the ocean depth over the first 120 miles. 
Munk, Snodgrass & Carrier (1956) have shown that since the wave motion 
under investigation is sensibly zero farther than A/27 from shore, it is essen- 
tially independent of the manner in which the ocean depth varies at these 
large distances. Since the phenomenon in which we are primarily interested 
is essentially confined to the continental shelf, the approximation of constant 
slope is adequate. ‘The total linearization of the actual topography is then 
a straight coast line and a constant slope depth (figure 2). 
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Figure 2. Linearized ocean depth and coast line topography. 


With A = xy, equation (3) becomes 


Cr 1 Cn 1 C P 
yV2n+ — => -9,(x,y,t) = ——(— +yV?P}, (6) 
of ag ol” ai pg\oy ~ 
and equation (+) becomes 
C db ] C “4 1 c I ? 
vV2d +  — seats go(X, y; t) ars: (eS (/ ) 
. ct ag of° ‘ pon ¢ f 


Given that P(x, y,¢) = 0 for t < 0, the correct boundary conditions are: 
n(x, y, 0) = 0, en(x,¥,0)/et = 0, n(x, y, t) + 0 as |x| -- & for all positive y, 
and also as y ~ for all x (the phenomenon 1s confined to the shelf), and 
finally that the wave height is finite at the shore line. ‘This finiteness condition 
is the most liberal physical requirement one can apply in the coastal region 
where the linearized differential equations are not a valid description. 

Since equations (6) and (7) differ only in the form of the forcing function, 
we need only solve the former; the solution to the latter is obtained by 
replacing » by ¢ and q,(x, vy, t) by g.(x,y,t). We can now take Fourier and 
Laplace transforms of equation (6). Defining 


f(R,¥,8) | e' dx | e~* n(x, y, t) dt 
| ; (8) 
O,(k, y, 8) e“* dx | e—*q(x,y,t)dt, 
where 7 = 1,2, we obtain , 
‘ d*f df {s? ae oe 
(— +k) )f = —O,(k,y,5). (9) 


~ dy* dy \ xg A 
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The boundary conditions become |f(k,0,s)| < oo and f(k,y,s)>0 as 
yo. Anticipating difficulties with the inverse Fourier transform, we 
shall solve, instead of equation (9), 


df d 2 
yi + +o (= +hy)f = ~0,(k, 55), (10) 


where Re = (R? +67)!2, (11) 


subject to the same boundary conditions. After taking the inverse Laplace 
and Fourier transforms of the solution of equation (10), we shall let ¢ tend 
to zero to obtain the solution of equation (6). If it is convenient and 
proper, we can also perform this limiting process at some intermediate 
step. ‘The function ky is not single-valued in the complex k plane and we 
restrict its possible values by two branch lines originating at +7e and 
extending to +700 respectively. The reasons for this particular choice 
will become evident shortly. We choose the Riemann sheet for which 
ky > 0 on the entire real axis. 

A few simple substitutions will enable us to rewrite equation (10) in a 
more convenient form. Let 


f(k¥8) = exp(—he ¥) A099), | 
u = 2key, (12) 
Pilk, u,s) = p(k, du/Re, 5), | 


and equation (10) becomes 


. . 
up; +(1 _ u)p; = (5+ ra) = ~ SPs) su) O i(R, 4 yu/ /Re, 5). (13) 


The boundary conditions become exp(—}u)p,(u)->0 as u—oo and 
|p,(k, 0,s)| < co*. 'To solve (13), we expand 7,(k, y, s) and 


_ exp(3 u) 
2Ry 





O,(R, 4u/Re, s) 


in terms of Laguerre polynomials (see Courant & Hilbert 1953), and 
then use the differential equation to determine the unknown coefficients. 


The expansions are 





pile) = pilkuss) = SALLA(u) 14 
yy 1 oO 
— 2B) 6,(%, Julke,s) = 3B, Ll) (15) 
Se n=0 
where L,(u) = e*( Fn) (ure* q (16) 


* We can consider k and s to be real numbers until it is necessary to invert the 
Fourier and Laplace transforms. 
F.M. 2Q 
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and ; \¢* 
A, = iaoe | e~“p,(u)L,(u) du, | 
heen F PA) oy ile, AE sie “ 
n (n!)? Jy 2k, pmo ola I Saal a n ° J 


Substituting these expansions into equation (13) and using the fact that 
L,,(u) satisfies Laguerre’s differential equation 


uL’(u)+(1—u)L,(u)+nL,(u) = 0, (18) 
we obtain 
5 4,(-n—-1- po) Lalu) = & By Lalu) (19) 
0 2ughy 0 
We use the orthogonality relation 
| e*L,(u)L y(t) du = (n!)Bnn (20) 
/0 
to equate coefficients, and hence determine that 
—-B 
A,, = Bn (21) 





(n+ E+ /2aghy) 


This implies that 








“A B,, L,(u) - 

= = : : - 
pi(u) pilk, u, s) 2 (n+ 1 4. 52 2ugks) (22) 

(p,(u) satisfies the boundary conditions, since each L,,(u) does) or 
, | L,(u)exp(—1u)O,(k, 3u/ky, s) du 
hw we conlnwkaak S 22 

Fby8) = exp(— hey), a Fa F[lagh, akan al?#ed)- 

(23) 


The wave height is then obtained by inverting the Fourier and Laplace 
transforms; i.e. 
> oO 7. p= + & 
Py ’ _ —Ar 
n(x, y,t) = rar e dk | 


ae 
tJ — co “¥ -10+8 








e"'f(k, 2, s) ds. (24) 


77 


The path of integration in the & plane for the inverse Fourier transform 
must lie within a strip containing the real axis. ‘The function ky has been 
defined with this in mind; i.e. the inversion contour can pass between the 
two branch points +7ze without crossing a branch line. 

There are other ways to solve equation (13); methods involving the 
solutions of the homogeneous equation (i.e. confluent hypergeometric 
functions) lead to extremely difficult integrations. A Green’s function 
technique for solving (6), although of little practical value, reduces to an 
interesting eigenvalue problem, one with both a discrete and a continuous 
spectrum; this will be the subject of a future paper. The advantage of 
the method used lies in the fact that we shall be able to consider separately 
the effects of a pressure distribution on the individual edge wave modes, 
easil’ eliminating extraneous parts of the solution. The wave motion 
can be determined with little difficulty. 
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3. APPLICATION TO THE MODEL OF Munk, SNopGrRAss & CARRIER 

Let us now consider this model, adjusted so that it arises at some finite 

time. ‘The surface pressure is 
Pya(y +a) 

P(x, y,t) = . , H(t 25 

(#95) (x— Ut)? +(y+a)/ (), (25) 

where H(t) is the Heaviside function, U is the velocity at which the disturb- 

ance moves parallel to the coast, and a is the half-pressure radius. ‘The 

resemblance to; actual hurricanes is poor, but the Fourier transform is 

simple. We shall solve for the wave height directly. ‘The forcing 
function is, from equation (6), 

a Pia (x+Ut)?—(y+a)? 
1% Vs 0) = pg [(x— Ut)? +(y+a)*]?’ 


and its transform is 

















exp(3u) 7P,a k, exp(—ak,) 
Oil, 3u/Re, 5) ) SIRO YOR , 

which is constant, so that 
aP,a k, exp(—ak,) 
pg 2k,(s—ikU) ’ 
Thus the expansion of equation (22) reduces to one in which only the 
‘first’ Laguerre polynomial, L,(u)=1, appears. This implies, as we 
shall see, that this pressure distribution can excite only the fundamental 

edge wave mode. The Fourier Laplace transform of 7 is 


BR, = and B, = 0 forn <0. 








—) _ -yky exp{—ke(y +4)} 
ae (s—tkU)(s?+agky) ’ (26) 
where 
amP, a 
=-——. 27 
y : (27) 


The wave height is determined by taking inverse Laplace and Fourier 
transforms as in equation (24). The function f(k, y,s) is not meromorphic 
in k, but is in the variable s, having simple poles at s = kU, s = +i(aghs)"*. 
It is, therefore, much simpler to compute first the inverse Laplace trans- 


form. We set 
1 4x +, 


r(k, y,t) = =| f(k, y, s)e* ds, 


2mt) —iw+e, 


and upon performing the integration we obtain for t > 0 


exp(tkUt) | 
r(k,y,t) = —yke exp(—hely+0)}{ Sta + 


4 explilaghe)} expla ilaghe)!t} | og 
+ Taghe)"8iagky)=—RU} * Taghe)ogk,y=+koyy °°) 


and r(k, y,t) = 0 for t < 0. 





2Q2 
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As a function of k, r(k, y, t) has no poles on the real axis. The solution of 
the problem is then 








(x,y, t) = - |. e*r(h, 95) dh, (29) 
TIT 
or 
n(x,y,t) = — (I+ 11+), (30) 
where 
r r exp{—k,(y+a)} exp —2k(x— Ut)} 
= k= : dk 
' J J. (ky U? —ag) ‘ 
we _ ( exp[—itkx—(aghy)'*t}]ky exp{—ka(y+)} 
see J fu = J, 2(aghy)**{(agky)K®— RUS oe 


ifkx+ (aghky)'?t}]ky exp{—ky(v+a)} dk 
2(aghy)'*{(aghky)!? +RU} 





III = | fan dk = jae 


Although the original integrand has no poles, and the path I’ may be the 
real axis, each individual term does have one or more poles; I’ must therefore 
be suitably altered. We can easily verify that (30) is a solution of equation (6) 
by formally carrying out the differentiations. 

Our task is now to evaluate and interpret these integrals. J is the 
easiest and we evaluate it first. This term represents the quasi-steady 
state solution derived by Munk, Snodgrass & Carrier (1956). (The terms 
due to a finite starting time are therefore the last two integrals.) The best 
contour to take is one that minimizes the effect in front of the pressure spot 
and produces a reaction behind it. Any contour will finally give a correct 
result for 7 provided we are consistent; however, we wish to give a definite 
physical meaning to J, independent of JJ and III. ‘The poles of f; are 
located at kg U?—ag = 0, ie. 


k 


r 

(32) 

k= -k = - 7 | 

for « = 0; k, is located in the right half plane, and &; in the left half plane. 
The best contour is one that goes under both of these singularities as in 
figure 3. For x > Ut, we may close the contour I in the lower half plane by 
the quadrants, T',, I}, the two sides of the branch line Ty, and the circular 
path about the branch point, as shown; for x < Ut we may similarly close 
the contour in the upper half plane. The poles contribute values for 
x < Ut only. The integrals about the branch lines are easily computed. 
The integrals over connecting quadrants tend to zero as |k| > oo, and the 





50) 
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integrals round the small circles about branch points tend to zero as the 
radius |k| tends to zero. The result is 


2 
I = [Re ¥(1,1, ikel]x— Ut] +i(y+a)]) x < Ut, 
4 ; 
I = 7exp{—ho(y +4)} sin ko(x— Ut) + (33) 


2 
+ paReF(1, 1, iko[|x — Ut|+i(y+a)]) x < Ut 


where V is the confluent hypergeometric function of the second kind (see 
Erdelyi 1954). Asymptotically, for large values of (x—Ut)?+(y+a)?, 





a to 
“ft, NS 
[ wit Ge ot | 
, = P ae o. 





Figure 3. Fourier inversion contours for X < or > 0. 


Y is proportional to the surface pressure P(x,y,t). The result consists 
of aresurgence for x < Ut, together with a motion restricted to the immediate 
neighbourhood of the pressure spot for which the wave height is the direct 
effect of the distribution on the ocean surface. 

We now show that JJ and JII essentially yield a cancelling wave, 
— (4n/U?) exp{—k,(y+a)}sink,(x— Ut) for x <}4Ut, implying that the 
resurgence is given asymptotically by 


27P, ak 
=_— - . exp{—k,(y+a)}sink,(x— Ut) for }Ut<x< Ut 


7 ~ 0 elsewhere, 


where k, is «g/U? as in (32). 
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The integral to be evaluated is 
r expl —ikx—(agk,)'*t}-—ka(yta 
q. | Sea 
J, Qagke)*{(agk,)®— kU} 
To obtain results which are applicable when x and Ut are both large compared 
to y +a, we use a steepest descent procedure; it is noted that if 





exp{—ke(y+a)} 
were not present, we could integrate explicitly. ‘The method of steepest 
descents is used to remove exp{—ky(y+a)} from the integrand, and the 
remaining integral is then evaluated. 


Let x <0, t= 4/B, then II = { exp{xf(k)W(k) dt, 
IT 


where 
(agk,)'? exp{—k,(y+a)}k 
= — * Ree cet es < K(k = = — Ty? 
f(t) = i(k ES), and Mh) = serie Sgke AO] 
We are interested in x and Ut large compared to y+a; accordingly 
exp{—s(y+a)} is included in the slowly-varying function %(k). The 
saddle point, determined from the condition f’(s)) = 0 is 











ag 
n= %. (34) 
The path of steepest descent is (if k = €+iu, and « = 0) the parabola 
ag \2 org . 
&2 analy Mh ees => — y 35 
ee 2B =) 








Figure +. Contour deformation of the inversion path. 


By considering the properties of the integrand f,,;, and, specifically, 
the behaviour of the factor 
exp[ — thx + (agky)'?—ky(y +)] 


tends to infinity in the third quadrant, we can show that 
| fudk =| frdk 
JT JF, 


(see figure +). ‘Thus far it has been essential to keep one branch line in 
the upper half plane, and the other in the lower half plane. The branch 
line from te to 700 can now be turned into the third quadrant (figure 5); 
and, for computing purposes, it will be taken to be the ray @ = 37, where 
the angle @ has its usual meaning. We use the contour I; about the branch 


k 
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line and the paths I’, and I’;, which are the arcs of a circle of radius R, 
to form a closed loop (in the limit as R-> oo) with the path of steepest 
descent. Cauchy’s theorem states that 


$ fu dk = 2ni > res f 77, (36) 


and thereby relates [. f1,dk to the integral over the path of steepest descent 
J In dk (the eveeiiaun of integration are shown in figure 5). The integrals 
S.D 


over the connecting arcs tend to zero as R tends to infinity, with the result 
that 


| ; fiz dk = | - fra dk + | _firdk+ 2ni ¥ resfin (37) 


where the residues are evaluated only at those singularities contained within 
the closed contour. The pole of f;; (k, = ky = «g/U?) is within the contour 
only when s) > k,, i.e. for 8B <}U orx <}Ut. For x large, the integral 








ao, 








Figure 5. Conversion from the fundamental contour to the path of the steepest descent. 


about the branch line is negligible compared with the integral along the 
path of steepest descent in the neighbourhood of the saddle point, and 
with the residue. (This is the advantage of the ‘deformed branch line- 
steepest descent’ method over other more direct estimations.) Hence we 


obtain 
on | fin dk (38) 
S.D. 


2z1 - : 
— Fa expl — ih x+i(aghe)!t— ho(y + @)JH(SU- B), 


where H is the Heaviside unit function. If the integral over the path of 
steepest descent were evaluated by conventional means as detailed by 
Morse & Feshbach (1953), the result would be 


THY 1/2 1 
| OFT ——__—__~ exp[ — dag(y+ 24 livox/B2—Ilin] — 
” (=) 2xg(1 — U/28) exp[— }ag(y + a)/B? + fiagx/B? — fiz] 


2zt 


od Ta exPl — tho x+i(aghy)'*t—ko(y +a)]H(}U—8). (39) 
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In the neighbourhood of s, = k, (8 = 3U or x = 3Ut?), this approximation 
fails, since the pole of the integrand is too near the saddle point, i.e. the 
‘slowly’ varying function is the dominant factor. We can avoid this 
difficulty by slightly modifying the conventional technique. Instead of 
evaluating the entire integrand at the saddle point, we evaluate only part 
of it, i.e. exp[—Ae(y+a)]. We obtain 

k, exp[ —tkx + i(agk,)!*t—k,(y +.)] dk 

J.D. 2aghs)'*{(agky)'?— kU} 
r= ok exp[—thkx +i(agk,y)"*t] dk 
~ exp[—5(y +4)] | 2(agh,)"2\(agh,)!2 RU} 








“ OW. 


= exp[-s(y+a)] | grrdky say. (40) 


The last integral is equivalent to an integral over the contour I’; (figure 6), 
provided we again add the appropriate residue term when s) < kp or 
2 <43Ut; Le. 


Qe 





J gp Se = | gu dk+ Tyexpl—ihox + i(aghy)*t}HU—p). (41) 
So 
2 Ky 
ae i 
y i Sg as et 
r, S.0: 





Figure 6. Equivalent inversion path. 


The integral over ' may be evaluated by using formula 809 of Foster & 
Campbell (1954): 
exp[ —tkx +1(xhk)! 72] 


_o 1-Ujagy= “ 








= Fa exp[—ihy(x— Ut)] erfe[(ag/U*ix)!*(3 Ut—x)] + 


ri 
+ (ag/imU?x)!? exp[i(gat?/ky x)] — = exp[—ik,(x—Ut)]. (42) 
Asymptotically, the final result becomes 


_ 
II~ 7 exp[—59( y + a) —tko(x — Ut)] {3 erfe[(ag/U*tx)"*(4 Ut —x)] —1}+ 


+ Fa H(QU—B)expl—ih(x— UN] {expt — sy +.a)] —expl— Ay +a)]} (43) 
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where s, = ag/48?, ky = ag/U?, t = x/B, and 


2 ru 
erfou = 1— —| exp(—v?) dv. 
V7) 9 

For x >4Ut, we have JJ ~ 0; for 0 <x < 3Ut, 

II ~ —(2mi/U*) exp[—ko(x— Ut)—ko(y +)]. 
Thus, for x < }Ut, IJ just cancels that contribution of J which arises from 
the pole at k,. It may also be shown that JJ = JJ* (complex conjugate). 
The wave height is therefore 


7 = — i (I+2Rell); (44) 

and the wave height of the solneaiebe aibieaiaaii is given asymptotically by 
Nw ~ 9, x> Ut; ) 

Nw ~ es exp[—)(y +a)] sinky(x— Ut), 4Ut <x < Ut; } (45) 
Nw ~ 9, | x < $Ut. |} 


The other contributions to the solution 7 do not contribute significantly 
to the resurgence. 














Name CAROL EDNA — — 

Date 30 Aug.— 11-12 Sept. | 14-15 Sept. | 21-22 Sept. 
1 Sept. 1954 1954 1944 1938 

Velocity U (knots) 32-34 32 33 40 





Wave Period T (hours) 


comp.: « = 5-0 x 10-3 5-8-6°1 5:8 6-0 43 
a = 4:2x10-3 6:9-7-2 6-9 7-1 8-6 

obs.: Atlantic City 5°5 6-0 5-6 — 
Sandy Hook 7:0 7-0 72 8-0 





Duration 


comp. 16-24 17-24 11-12 9 
obs.: Atlantic City 20 23 23 —— 
Sandy Hook 26 _ 30 16 























Table 1. Periods and durations of sea level disturbance due to four hurricanes. 


The resurgent wave motion induced by this simple pressure distribution, 
equation (25), is essentially restricted to the interval }Ut < x < Ut, and 
is characterized by a single wave number k, = ag/U? (that of the fundamental 
edge wave mode) or the corresponding period T = 27U/(ag). The dura- 
tion D, which is the time taken for the resurgence to pass a given point, is 
equal to the total time ¢ for which the disturbance has been in existence 
and in motion. ‘The periods and durations of the resurgences of several 
hurricanes have been computed by Munk, Snodgrass & Carrier; the 
observation points were Atlantic City, where « = 5-0 x 10-%, and Sandy 
Hook, where « = 4:2x10-*. Their results are presented in table 1. 





586 H. P. Greenspan 


Although the model of equation (25) bears little resemblance to an actual 
disturbance, the results are in agreement because, evidently, hurricanes 
also excite only the fundamental edge wave mode (see §4). ‘The computed 
periods of the fundamental mode are in excellent agreement with the 
observed periods. 

Mathematically, the existence of a resurgence is due to the fact that the 
integrand of the Fourier inversion integral has poles at k, and k, The 
analytic form of the waves is obtained by evaluating the residues of the 
integrand at these poles. ‘These results may be generalized to the case of 
an arbitrary pressure distribution which affects all the edge wave modes. 
The analysis of each integral arising from the inverse Laplace transform 
of the infinite series, equation (23), is substantially the same as that already 
performed in equation (26) et seq. The integrals involving the mth Laguerre 
polynomial give rise to the mth edge-wave mode. The analytic form of the 
resurgence is again obtained from an evaluation of the residues of the 
Fourier inversion integral. ‘The final result is that, for disturbances 
originating at ¢ = 0 and travelling parallel to the coast with constant velocity 
U, the resurgence, which isasum of all the edge-wave modes, is asymptotically 
given by 


Nw ~ 9, a> 42; 


Nw ™ y A,(k,,) exp(—k, y)L,(2k,,y) sink,(x— Ut), }Ut <x < Ut; } (46) 


0 

Ny 0, x < }Ut, 
where exp(—h,,y)L,(2,,y)sink,(x— Ut) is the nth mode, its amplitude 
being A,(k,,), k,, = (2n+1)ag/U? = (2n+1)ky. The resurgence is confined 
to the interval }Ut < x < Ut; the front moves with the speed of the storm 
U, and the rear moves with group velocity }U. ‘The duration is the same 
as that previously calculated. 

We may now apply these results to a more exact mathematical model, 
no new techniques being needed. 


4. ‘THE GAUSSIAN PRESSURE DISTRIBUTION 
'l'o investigate how the various modes are stimulated by a real disturbance, 
we need a more realistic model. We consider, therefore, the effects of a 
travelling pressure taken to have a Gaussian distribution given by 


P(x,y,t) = Py exp| - — ea), 


a 
where P, is the maximum pressure deviation, and a = 1:19r,, ry being the 
half-pressure radius. It is advantageous to use the potential-function 
formulation of equation (4). The wave height is then given by equation (5). 
On transforming to a moving coordinate system (x’ = x— Ut, y’ = y, t' = t) 
and then dropping the primes, equation (4) becomes 


(SS o? ) ob «61 eh QU Oh U oP 





(47) 











+= — — = + ——_ a cae 48 
oy ag of? ag ox? pga cx (48) 


se 5) 
ox dy? 
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In this coordinate system, the pressure is independent of ¢. By taking 
Fourier and Laplace transforms, we obtain 


12 ) s+z2Uk) 
EP, Af pe, C+iUR 
de dy |? xg 
Utk ere 2L9/A 
= P, exp[ — (vy —9)?/a?] avz exp( — a?k?/4), (49) 


apgs 





where 


r a0 roo -_) 
® = | etdx | etdbdt=| e$(k,y,t) dt. 
— “0 “ —o 
The solution of this equation (we actually solve a modified equation with 
yk%, in place of yk?) is, from equation (13), 
2vmtU Py ak, 


-. A/(R)L,(2Rg V 
rt) SOS exp(—ky) > nh )L sf ky) 
p 0 


n!)*s{(s+iUk)? +(2n+ I)ky ag} 





where 
Aj(k) = kexp(—a?k?/4)| L,,(2R« y)exp( — ke y)exp{—(v —¥)?/a?} dy. 
J0 
(51) 
The coefficient A/(k) has no poles in the finite k-plane; its branch line and 
branch points coincide with those of ky. ‘Taking an inverse Laplace 
transform, we find that 
“* 2v7iUP, ak, exp(—ky vy) & AL(R)L,, 
O(k, yt) = — : a 


Pp 1) 
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We have written only the first term of three in the brackets, because, as we 
have mentioned, it may be shown that this term gives the wave motion for 
x < Ut; the other terms add to yield a cancelling wave in the region 
x <1Ut. 'To determine the induced modal amplitudes, we need compute 
only the wave function arising from the first term, since all additional terms 
yield the same wave. 

To compute the amplitudes of the waves constituting the resurgence, 
we perform the inverse Fourier transform and evaluate the residues arising 
from the poles, which are at 

k= +k, = +(2n+1)ky = +(2n+ 1)ag/U?. (53) 


In the original coordinates, the result is 








$(x,94t) = S4(x,9,t)+ | Fdk,  4Ut <x < Ut; | 


$(x,Vv,t) = | F dk, elsewhere, | 
where ¢,, is the wave motion of the mth mode, determined, as we have said, 
by evaluating the residues of the integrand at the poles k = +k,; thus 


2 2 v7UP,a\ ,, , 
¢,(x,y,t) = — TAA — Ai(k,,) cosk,(x— Ut)L,(2k,, y)exp(—k,,y)e 
(55) 
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The integral in (54) represents all the other terms, which do not contribute 
tothe resurgence. Since the wave height is related to the potential function 
by equation (5), the amplitude of the mth edge-wave mode, as defined in 
equation (45), is 

2 


sa 2 2v7r UP, : 
A,{Rn) — > a(- et) Ail) (56) 


The amplitude ratio of the nth mode to the mth is 








A, (m!)?k 7 exp( —a°k’, /4) 
A. ~ exp( ~@hk? /4) x 
[ “6 L,(2k,, y’)exp(—R,, v’)exp[ —(9’ —¥)?/a?] dy’ 
7 . CF 


| Ly(2km ¥’)exp(—Rm Y’)expl — (y" —¥0)?/a%] dy’ 
0 


The amplitude ratios A,/Ay and A,/A, are plotted in figures (7) and (8) 
as functions of kya for yy = 0 (storm centred on the coast line). 
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Figure 7. Amplitude ratio of the first harmonic to the fundamental vs kp a. 
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We shall interpret these results in terms of the half-pressure radius rp 
and the speed of propagation U. For y=0, P= Pyexp(—r*/a’), or 
r = 0-837a. For hurricanes travelling parallel to the east coast of the 
United States* at 32-34 knots, which is typical, k>' lies in the range 
30 < k=’ < 40 miles. The half-pressure radius of these disturbances is 








3.0 




















Figure 8. Amplitude ratio of the second harmonic to the fundamental vs kya. 


of the order of 100 miles, which means that a = 120 miles, or that 
3 <k,a<4. From figures 7 and 8 it is immediately apparent that for 
these values of kya only the fundamental edge wave, 


Lo(2k 3) exp(—Ry 3) sin hy(x— Ut), 
is excited. All other modes are negligibly small. These results also show 
that for a disturbance of hurricane proportions to be an effective generator 


* Henceforth we shall refer exclusively to hurricanes in this location. In analysing 
the effects of disturbances off other coasts, the appropriate values of «, U, a, etc. 
must be used. 
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of the first harmonic, it must travel at 70 knots at least. ‘To excite the 
second harmonic, it must have speeds of 120 knots or more. We may 
draw the conclusion that hurricanes will never effectively stimulate any 
other mode except the fundamental edge wave. For smaller disturbances 
the higher modes may be generated. A storm with a half-pressure radius of 
approximately 25 miles and a speed of 34 knots has a value of kya between 
0:75 and 1; k, is thus located exactly in the region where higher harmonics 
are stimulated. ‘To induce still higher harmonics, disturbances must 
travel faster or be smaller in size or both. 

Let us investigate the fundamental mode, and determine the conditions 
(size, distance from shore, etc.) for optimum edge-wave generation. The 
wave form of the fundamental mode is given by 


2nP : . , 
No = | = ° (ky a)? exp( — ky a) erfe( — _* 5) exp —kyy)sinky(x— Ut). 
(58) 





The bracketed expression is the wave amplitude. To find the storm 
position y, that produces optimunt amplitudes (with constant kya), we 
must determine the value of yp, say yf“, which maximizes 


; ¥ ; ko 
fo) = exp(—Fo ¥o) ere( — = =), 


a < 


that is, for which f’(v,)) = 0. Hence, we need to solve 


>) . 9 
2 exp(—u?) 
— “ aoe where u = —y,/a+}kya. 
i eal to ¢ i 7 


For k, a = 2, 2:5, 3-0, 3-5, we find u = 0-6, 0-9, 1-2, 1-5. In general these 
results yield maximum wave heights for y) in the range 25 to 35 miles. 
Table 2 summarizes these results, giving the maximum wave height in 
terms of the coefficient 27P)/(pg). For a maximum pressure difference of 
one pound per square inch, the resulting wave heights, in the hurricane 
region k,a= 3-5, are approximately 2-5 ft. It is seen that the wave 
height measured at the coast varies slightly with the distance of the storm 
centre from shore, being 2-3 ft. for vy) = 0, and rising to a maximum of 
25 ft. 

Knowing that maximum wave heights are induced by hurricanes 
24 to 35 miles from shore, we can now determine what size disturbance 
in this region will produce the largest effects. We maximize the wave 
amplitude with respect to the variable a, holding ky and y, fixed. This 
leads to the equation 


vr erfc(—yo/a+ 43k) a) = (¥vo/a+ 4k, a)exp[ —(yo/a— 3k, a)?]. 


For ky = 0-025 and y) = 30, a = 68 miles, and the half-pressure radius is 
r, == 54 miles. Our conclusion is that hurricanes are too big to produce 
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maximum wave heights; ‘smaller’ storms (ry == 50 miles) 25 to 35 miles 
from shore produce the largest waves (for a given pressure deviation P,). 




















Pe Maximum one Storm centre 
me amplitude —- (ymax) 
20 | 0-712 27Py || 1/30 & 60 22°8 
pe || 1/40 & 80 30-4 
2 . FE Le 
2-5 | 0-529(,,)4 2,2 te 
1/40 & 100 34-0 
1/30 & 90 27:0 
3-0 | 0-328 pies a 
, 325 (5) 1 1740 & 120 36-0 
= ia 1/30 & 105 26-2 
a5 0-170:¢ .. : sts 
' OC) 1/40 & 142 37-0 
In particular, if Po 1 lb/in.*, and pg 


62-5 lb/ft.3, then 27Po/(pg) 14-5, and the 
following wave heights are computed for 
various distances V9. 





] | 
kya ¥o 0 Vo = ymax Vo 100 
2-0 10-4 — 
2°5 - 7:65 — 
3-0 4-44 4°75 2:7-3°7 
3°5 2:3 2:46 1:65-2:02 




















Table 2. Maximum wave heights in fundamental mode : 
a and yp are in miles, and wave heights in feet. 


5. CONCLUSION 

We can now predict the effects of any particular distribution on the 
resurgence from a knowledge of the maximum pressure difference, half- 
pressure radius, and speed. One should hence be able to calculate the 
positions of highest tide along the east coast following hurricanes. While 
the amplitude of the resurgence is not in itself remarkable, the extremely 
long periods enable these waves to be effectively superimposed upon an 
incoming tide; this is a dangerous situation, especially in those areas where 
there are resonant phenomena. A two foot or three foot crest coming 
several hours after the passage of a storm (and therefore somewhat unex- 
pectedly), and superimposed upon the incoming tide, may greatly enlarge 
existing resonant wave motions (in bays, estuaries, etc.). While these 
situations cannot be avoided, it may be possible to reduce their effects by 
predicting their occurrence. 








592 H. P. Greenspan 


It should be noted that in order to simplify the mathematical analysis 
the effect of the earth’s rotation has been neglected. ‘This approximation 
is quite proper for the study of wave motion whose period is a few hours 
or less. However in the application of the theory to hurricane resurgences 
it is estimated that errors of 15°, are probable. That we were able to 
obtain such excellent and consistent agreement between observation and 
theory with so glaring an omission is, to say the least, peculiar. At 
present those additional effects which presumably nullify the rather large 
contribution from rotation are unknown. 


The author wishes to thank Professor George Carrier for suggesting this 
problem and for his assistance in preparing this article. ‘The paper is 
part of a thesis submitted in partial fulfillment of the requirements for the 
degree of Doctor of Philosophy in Applied Mathematics at Harvard 
University. ‘The work was sponsored by the Scripps Institution of Oceano- 
graphy under contract from the Office of Naval Research. 
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Irrotational motion associated with free turbulent flows 
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SUMMARY 

The irrotational motions induced by, but outside of, the self- 
preserving turbulent wake and jet are examined. It is found that 
there is a mean flow towards the centre of the jet, although there is 
no such flow in the case of the wake. Phillips’s (1955) results on 
the nature of fluctuating irrotational flows are found to be largely 
unaffected by the introduction of simple inhomogeneities into the 
boundary conditions, but another mode of fluctuation with move- 
ment along the lines of mean flow is also shown to be likely. It is 
pointed out that, contrary to a statement of Corrsin & Kistler 
(1954), it is possible for the vorticity-free fluid between bulges of 
turbulent fluid to partake of the mean velocity of the turbulent 
fluid. 


1. INTRODUCTION 

Experimental researches in recent years have established beyond doubt 
that in all those flows, such as wakes, jets and boundary layers, where 
turbulence is expanding into non-turbulent regions of fluid, there exists a 
relatively sharp, although violently contorted, boundary separating the 
fully turbulent region and an external region. ‘This external region, 
although not stationary, is characterized by very much smaller velocity 
gradients than is the turbulent region. Corrsin & Kistler (1954) and 
Townsend (1956) have pointed out that since vorticity can be introduced 
into non-vortical fluid only by viscosity, which is a very short range force, 
the motion in this external region is probably irrotational, and therefore 
can be derived from a scalar potential. 

Phillips (1955) has considered a simple case of such irrotational motion, 
taking boundary conditions such that the velocity component u, is a 
stationary random function of the cartesian coordinates x,, and x3 on the 
plane x, = 0* and that the velocity u vanishes as x, > 0. It is the purpose 
of this paper to extend Phillips’s results to some situations in which the 
statistical properties of uw, are not independent of both x, and x3, and to 
offer some comments and speculations on the phenomena. 

* We have changed the order of the indices from that used by Phillips in order 
to conform to customary usage in wakes and jets. wx, is the principal direction of 
flow and x, is parallel to any direction of symmetry. 


F.M. 2R 
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We begin with a review of Phillips’s analysis. Given that 
(u2),, =0 = | e&-* dA(k), 
where k is a (two-dimensional) vector wave-number with components 
(k,,0,k,) and magnitude k, x is the position vector with components 
(x,, X., ¥3) and the integration is over all k-space (two-dimensional); and 
defining 


where the bar indicates an average over the plane x, = 0 and the asterisk 
denotes the complex conjugate; he finds that, for all x, > 0, 


uz = | O(k)e* dk, 


a= - R2 " (1.1) 
=| 72 ke ahr dik (¢ = 1, 3). | 
It is shown that, since 6(k) is an even function, 
— - | a 
UU, = | = (ee 2 dk = 0 = uy us. (1.2) 


He also states that u,u, = 0, and indicates that a similar proof would 
establish this relation. ‘This is not the most general case, however, for in 
fact 
"kk; 9 

k? 
which is not generally zero. Since u,u; is a symmetric tensor there will 
always, of course, be a possible choice of the x, and x3-coordinate directions 
for which u,v, = 0. In the ordinary cases of the two-dimensional wake, 
jet or boundary layer, this choice corresponds to the ‘natural’ one of x, as 
the downstream direction. For the wake of a yawed cylinder or the boundary 
layer over a plate whose leading edge is not perpendicular to the flow, 
however, u, u, ~ 0 if the x,-direction is chosen to be the principal direction 
of flow. ‘There is no need to exclude such examples, and we may continue 
the analysis in the manner of Phillips. 

Assuming that 6(k) may be expanded near k = 0 in the form 

Ok) = 04+k,0,+k,k,0,,+k,k;k, 8,4... 
where 6, 4, 0,;, etc. are cartesian tensors of the indicated order and are 
independent of the &,, Phillips shows that @=6,=0. At large values 
of x, the values of 6(k) contributing significantly to the integrals in equations 
(1.1) to (1.3) are given approximately by 


(k)e-2*"> dik, (1.3) 





Uj;us = | 


A(k) ~ kk; 9,;- (1.4) 
Putting (1.4) in (1.3) we get 
Uy Us ~ | 1 hy iyo? dl 
= | (6540) 2d (1.5) 


R2 
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since terms in 6,, and 6,, evidently vanish for reasons of symmetry. ‘There 


is no loss of generality if we take 6, to be symmetric, so 
> d a] a , 
- R2k2 


uu; ~ 2 | ey 6,,6°7** dk 
$5 770.5.X6 *. (1.6) 
We also write down for reference relations derived by Phillips: 
us ~ 370), +453)xz 4, 
uy ~ 37(30,, + O33)x3 4, f (1.7) 


os 357(O14 v 3653)x3 4. J 


~ 


~ 


u’; 
Thus we see that in general there zs a Reynolds shear stress due to the 
irrotational motion, with its direction parallel to the plane x, = 0. 

It is evident from (1.1) that 


uz = ur +u2, (1.8) 


so that the velocity fluctuations are not isotropic. ‘This is of some conse- 
quence, for in the following sections we shall consider situations where u, 
is not statistically stationary on the plane x, = 0, and it can be shown that 
this lack of istropy then gives rise to stress gradients. ‘Thus, in the absence 
of viscosity the total cartesian stress tensor in an incompressible fluid is 


7; = P8;+uu; (1.9) 


7; = 
where P is the pressure, 6,; is the Kronecker delta, and uj, U; is the Reynolds 
stress. If the fluctuating velocities are distributed isotropically, this 
becomes 
7; = (P+u*)6,;, (1.10) 

and it is possible for gradients in u2 to be compensated by gradients in P, 
so that no net stress gradients (or accelerations) need result. In the absence 
of istropy this compensation becomes impossible and variations in turbulent 
intensity are necessarily accompanied by stress gradients. 

In addition, Corrsin & Kistler (1954) have shown that in the special 
case of irrotational flow the relation 


Ou; OU; 

mae ith isgeel (1.11) 
OX; OX 
imposes conditions on the Reynolds stresses which have some interesting 
implications. For example, expression (1.8) can be derived very simply 
from the boundary conditions and the relation (1.11). ‘Thus the Reynolds 


stress gradients muxi satisfy the relations 


i 








e 
=— Ul; = u 
Ox; 








Ous 
: Aa 
Cx; 


But the fluctuations are stationary in the x,- and x;-directions. Hence 
CU, Uy = OU U: 
+? = += = 0, 
Ox, OX 
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and SO Cc is 1 C 
SS 
OX, 2 0: 2 : 
Then, since u vanishes as x, > ©, 
mS 2 
us = tu;u; = uy+u 


2. CYLINDER WAKE 

The predictions of Phillips’s theory have been remarkably confirmed by 
comparison with measurements, made by A. A. ‘Townsend in the wake of a 
circular cylinder. It is probably worth examining the potential flow 
associated with such a wake in more detail. 

We shall consider only a self-preserving flow which Townsend (1956) 
has shown to be possible only if the turbulent velocities are small compared 
with the mean flow velocity U,. In this case the characteristic scale L 
varies as x!* when the x,-direction is the principal direction of flow and 
the origin of x, is some virtual origin not too far from the actual iocation of 
the cylinder. It is also shown that, for the transport of momentum to be 
independent of x,, U, « xj", where U, is some characteristic velocity 
which may be taken as U, = U,,—(U)min for any given x,. Root-mean- 
square turbulent velocities vary as U5. 


— Uo 









> 
————— OD 


o(L+dL) 





: 
8 


Figure 1. The self-preserving wake. 
The scale in the x,-direction is greatly expanded. 


In figure 1 the total mean transport of fluid across the section AA per 
unit length in the x-direction is 
al(U,—KU,)+adL U,, 
where AU, is the average value of U,,— U, across the wake, the total extent 
of which is aL, and K and a are constant. The total mean transport of 
fluid across the section BB is 


a(L+dL){U,,—K(U,+dU,)}. 
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Neglecting second order differentials, the transport across AB in the 
x,-direction is therefore 
laK(LdU,+U,dL) = 4aKd(U, L). 
Hence ee 
U, «< —(U,L) = 0. (2.1) 
Thus we find that there is no mean flow from the x-direction into a self- 
preserving two dimensional wake. We shall see later that this 1s not true 


of a jet. 
y 


3. FLUCTUATING MOTIONS OUTSIDE A TURBULENT WAKE 
Returning to the case considered by Phillips, let us assume that it is 
possible to define 6(k, x,) for every value of x, on some plane* x, = 0 
outside, but close to, the main body of the wake. We assume that, although 
6(k, x,) is inhomogeneous in x,, Phillips’ main conclusions are not materially 
altered. 
Our similarity conditions require tha 


O(k, x,) = x(kx,"*) = x(e), (3.1) 


oF 


where x is a function of € = kx;* only, in order that the characteristic 
scales should vary as x} '~, and that the energy per unit mass in the fluctuations, 


O(k, x,) dk x; "| y(e) de, 


should vary as x7’. ‘The expansion 
O(k, x.) = R,R,O.(x,) HR, RR Ry, 8 stm (X1) +o 
may then be written 
O(k,x,) = x RR; xj + x7 RRR Rm Xijtin t (3.2 
where y,,;, etc. are absolute constants for the particular flow. 
At sufficiently large values of x,, we find from (1.7),(1.8) and (3.2) that 
us ~ 3a xu + X33)*1 =", 


ul ~ S7(3x1 + ¥33)%] Se > (3.5) 





By symmetry U, Us = Uau, = 0. 

It is interesting to note that these fluctuation intensities increase as x, 
increases. The reason for this, of course, is that the increase in the scale 
of the primary causative turbulence more than compensates for the decrease 
in intensity. If we consider the variation on similarity surfaces x, < xj"*, 
instead of on planes x, = constant, we find that the fluctuation intensities 
are all proportional to x,’, as they should be. 

The equations for incompressible mean flow in the steady state are 





cu, 1oP ¢— e2U , 
U,— +-=— + — uu, —v-——— = 0. (3.4) 
Om pom o% ” * Cx), OX, 


* The actual surface should be parabolic, but this introduces considerable 
complication into the analysis and at large values of x, the difference should be slight. 
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If we have Reynolds number similarity (Townsend 1956) which we shall 
assume, the viscous term will be negligible. In the particular region of the 
flow in which we are interested we have seen that 





U, ~ Us; = u, uz = Ug; = 0, 


and we shall see later (3.9) that (¢/0x,)u,u. = 0; so it would appear, reason- 
able to write for the x,-component 























pOX, OX,’ 
or P+ pus = P,, (3.5) 
P.,. being the value of P at large values of x,, where u5 is very small. 
For the x,-component, the ony terms of any consequence will be 
ou _0U, 10P 2 ou? 
U,— +U,— +-—+—u4,4,+ = = 0. (3.6) 
~ OXg OX, pox, OX, ~ OX, 
Corrsin & Kistler (1954) have shown that, for irrotational flow 
Cu, u 1 cq? = 
il ma, (3.7) 
ON; 20x, 
._ 
whe re 7 Gh Uj Uj 
If mean quantities are independent of x5, (3.7) gives 
O1us Le... <> 
nea Be an ig — 2u5), 
OX LOX, 3 2Q 
srner ie (3.8) 
Cu,u, 1 (q?— 2u8) 
—— = ~— (q?—2u}). 
a , l 
COX ON 
In our case, where u5 = uj +uz = 4%, this yields 
CU, Uy 
—— = 0, (3.9) 
( vy 
OU, Uy d= <<; 2 
——— = (us — uz), (3.10) 
Oxy ax, 


— = —(- ui). (3.11) 





0. (3.12) 


We therefore see that the mean flow is unatfected by the irrotational 
fluctuations, not because the Reynolds stresses are zero, but because the 
gradient of the shear stress (¢/€x,)u,u, is exactly balanced by the gradient 


of the total normal stress (¢/0x,)(P/p+uj). It would appear that such a 
balance is the rule rather than the exception in irrotational fluctuating 
flows. 








a) oo i oe 8=6«6 
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The normal stress gradient tends to accelerate fluid in the direction of 
the mean flow. However, the shear stress transfers this acceleration to the 
retarded flow in the main body of the wake. ‘The whole process must be 
considered as a part of the mechanism by which this retarded fluid is 
accelerated towards U,, as the wake decays. 


4. MOTION BETWEEN BULGES OF TURBULENT FLUID 

There is a difference of opinion expressed in the literature as to the nature 
of the mean flow in that portion of the external, irrotational flow region 
which lies between out-thrusts or bulges of the turbulent fluid. Townsend 
(1949) states that “‘the non-turbulent fluid between successive jets is 
constrained by pressure gradients to move at the same mean velocity as the 
fluid in the adjacent jets’, and shows measurements (Townsend 1956, 
p. 163) to support this view. Corrsin & Kistler (1954), arguing that the 
mean flow in these regions is irrotational and therefore gradient free, say 
“‘the mean velocity everywhere in the potential field must be constant and 
equal to that at infinity” and publish measurements qualitatively sub- 
stantiating this statement. 

An examination of Phillips’s results may cast some light on this contro- 
versy. It is noteworthy that equations (1.1) exactly correspond to the 
classical theory of flow in gravity water waves (Lamb 1932, ch. 9) if we consider 
a situation where there is a continuous spectrum of wave frequencies and 
waves may be oriented in all possible directions. Now it is well known 
(Lamb 1932, p.418) that, although such waves are purely irrotational 
flows, if they have finite amplitude they are accompanied by a mean transport 
of fluid, and a resulting mean shear. ‘The reason for this apparent contra- 
diction is that they are bounded by a non-planar surface, and any plane 
passing between the level of maximum elevation and the level of maximum 
depressionincludes regions in which the equations of irrotational flow do 
not apply. ‘The simple case of a single wave train with a single frequency 
is sketched in figure 2 (a), which clearly illustrates the phenomenon. Every- 
where below the boundary the flow is irrotational. ‘The velocity potential is 

¢@ = —Bcoskx,e*”, (4.1) 
where B is a constant, and the velocity components are . 

u, = Bhoankx,e**, 

ue = Bkooskx, e—** 
For any x, > 0, 


“u,= >— | sin kx,e~*"* dx, = 0, (4.2) 
7 1) k 


but, if the boundary* is described by 
x, = —h(1+sinkx,), (4.3) 


* The fact that the actual boundary of a theoretical water wave is trochoidal does 


not concern us here. 
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for 0 > x, > —2h, we find that, 
i. Bke Kis -{sin-3(—1—2./h)Vk 


uy =- sin kx, dx 
: cos-'(—1—x,/h) J 24 —s 


~~ 2%, & 
~ cos“ —1— v,/h) JG 7 3). ae, 


This function is plotted in figure 2(b) for B= k=h=1., 








e ane -2x,-x3 


cos'(-x,-1) 
-2 -—— 











O =-B coslkx,e **2 


(a) (b) 
Figure 2. (a) Flow in a simple gravity wave advancing from left to right. The 
length of the arrows is proportional to the logarithm of the velocity at their 
centres ; (6) the mean horizontal velocity within the wave. 





Figure 3. Probable directions of potential flow near the boundary of 
turbulent fluid in a wake. 


The irrotational flow induced by free turbulence is another example of 
potential flow with an irregular boundary beyond which Laplace’s equation 
does not hold. It is quite probable that a similar phenomenon occurs, for 
the pressure forces are likely to make the mean direction of the potential 
flow between bulges the same as that of the bulges themselves. The 
example of gravity waves shows that such a motion is mot inconsistent with 
the equations of irrotational flow. If the flow is of this type, we can make 
some predictions concerning the direction of uv, in the non-turbulent region 
(figure 3). In a wake this component should consistently be negative 





he 
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just before the arrival of a bulge and positive just after its passage. Perhaps 
the similarity in appearance between the boundary of the turbulent fluid 
and the surface of a tumultuous sea is more than merely superficial. 

5. Jers 

In considering the two-dimensional wake we passed rather lightly over 
the fact that our similarity surfaces were not planes x, = constant but were 
parabolic. This was justified at sufficiently large values of x,, since these 
surfaces there approach closely to such planes. For jets, however, such 
an assumption becomes very questionable. It is found (Corrsin & Kistler 
1954) that in jets characteristic scales vary as x,, and characteristic velocities 
vary as x; | for circular jets and as x; '* for two-dimensional jets. In these 
flows it is obviously not possible to choose suitably a unique direction for 
xX) in cartesian coordinates in order to apply Phillips’s theory. ‘The appro- 
priate coordinate system will be spherical for the circular jet and cylindrical 
for the two-dimensional jet. 

For the circular jet we have boundary conditions defined on a cone, 
with vertex at a virtual origin close to the source of the jet and angle such 
that it just contains all the turbulent fluid (about 25° included angle 
according to experimental results). In spherical coordinates (r, «, 8) this 
cone is defined by x = a», and following Phillips’s general line of argument we 
define a potential ¢(r,,8) such that u, =(1/r)(¢¢/ex) is a statistically 
prescribed function of r and 8 on x= 4%. (Uy), =, Will be characterized 
by a typical scale which varies as r and a typical magnitude which varies as 
r1, Such a function can be accurately and conveniently described by 

(Ux)aeax = > | tk, A(k,,kg)r'*r* exp(ikgB) dk,, (5.1) 
where k, and kg are dimensionless wave-numbers, the summation is over 
kg which takes only integral values, and A is complex. he potential 
function everywhere outside the cone will then be 


wn 


@=)5 1 A(k,, kg)r'“rexp(ikgS)b(, k,, kg) dk, (5.2) 


where ¢ will be expressible in terms of Associated Legendre Functions of 
integral order k3, pure imaginary degree 7k,, and real argument cos«. 
The function % may be analytically different for the various combinations of 
positive and negative values of k, and kg and will depend to some extent 
on other boundary conditions such as those that are imposed, for example, 
if the jet issues from a hole in an infinite wall. 

Unfortunately there appears to be no published study of such functions, 
and it hardly seems worth conducting a detailed study for the present 
application. 

6. MEAN FLOW OUTSIDE A CIRCULAR JET 

A similar argument to that used to derive (2.1) shows that the mean 

flow component U, is given by 
UG. = Cr“, (6.1) 
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where C is a negative constant on the boundary of the turbulent region 
E = Ky 
A solution of Laplace’s equation for the mean velocity potential ®, 
with (6.1) as a boundary condition and symmetry in f, is 
® = A, log(rsinx) + A, log(tan 3), (6.2) 
A, cos% +A, = Csin %. 
Expression (6.2) leaves one arbitrary constant to be determined by other 
boundary conditions. 
If the jet issues from a hole in an infinite wall at « = }7, then 4, = 0 and 
® = Ctan x, log(r sin x), (6.3) 
which corresponds to a very simple uniform radial flow directed everywhere 
towards the axis of the jet, with speed inversely proportional to the distance 
from this axis. 





Figure 4+. Flow lines into a self-preserving circular jet issuing from a nozzle. 


If, as in Corrsin & Kistler’s (1954) experiments, the jet issues from a 
nozzle, the boundary condition may be approximated by 





om 
—- ==) 2f aor 
on 
and we have 
C sin x , 
@ = ——— log(rsinxtan }z). (6.4) 
cos %, + 1 ; 


In this case the flow lines are given by 
r(cos%+1) = constant (6.5) 
and are sketched in figure +. It will be noted that, even here, when close 
to the jet, the mean flow is directed essentially perpendicular to the axis 
of the jet. 
7. ‘TWO-DIMENSIONAL JET—-MEAN FLOW 
By a ‘two-dimensional jet’ we refer to the type of flow produced when 
fluid issues from a long uniform slot. ‘The mean flow is two-dimensional, 
but fluctuating motions occur in all three dimensions. In this type of 





on. 
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flow the similarity conditions require that the scales vary as r and the 
characteristic velocities as r4?. In cylindrical coordinates (r, «, 2) this 
requires that 
U, = Dr” (7.1) 
on the edge of the turbulent region, « = %. D is a negative constant if 
% IS positive. 
The velocity potential satisfying (7.1) and also the condition U, = 0 
at x = 47 (corresponding to a slot in an infinite wall) is 
® = Dr'{cos(}o9+42)}"! sin($a+ jz). (7.2) 
The corresponding flow lines, which are given by 


r}2cos(4}a+47) = constant, (7.3) 


are sketched in figure 5. 








Figure 5. Flow lines into a self-preserving two-dimensional 
jet issuing from an infinite wall. 


8. FLUCTUATING POTENTIAL FLOW OUTSIDE A TWO-DIMENSIONAL JET 

I have been unable to find, in three dimensions, a suitable function to 
describe the boundary condition for flow outside a two-dimensional jet. 
However, if we can imagine a flow in which the fluctuating motions are 
also two-dimensional the problem becomes rather simple, and the conclusions 
derived from a study of this hypothetical flow may, with some caution, be 
applied to the more complex real flows. fo 

A velocity potential in polar coordinates (r, «), for which u7 varies as 
r-! and typical scales vary as r, and which satisfies the boundary condition 
u, = 0, at «= 47 (corresponding to a slot in an infinite wall), can be 
constructed thus: 
d=r2| Ase-'*cos(B+klogr+}xa) +e" cos(B+k log r + 3a — }x)} dk, 

- 0 
(8.1) 


where 4 and B are real functions of R. 
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Now in the region of greatest interest, 7— x = 3, which is of the order 
of ten times as great as x. ‘Therefore, for k greater than about one, the 
second term in this integral is negligible compared with the first. For 
small k, however, the two terms approach equality in magnitude. 

For k < 1, then, the velocity components reduce to 

ad na 
= ae '2sin(ka+ 4m) | Acos(B+klogr+47) dk, 
2 /0 
ee +o 
= eee L'2 cos(4% + 477) , A cos(B+klogr +47) dk. | 


(8.2) 





u, and u, are therefore fully correlated, and the only motions at these very 
small wave-numbers are along streamlines given by 


r}2cos(}a+47) = constant. (8.3) 


It will be observed that (8.3) is identical with (7.3), the expression for the 
streamlines of the mean flow into a two-dimensional jet. We have therefore 
come to the interesting, if not particularly surprising, conclusion that the 
very small wave-number fluctuations follow the same streamlines as the mean 
flow, shown in figure 5. ‘This result is expected to be valid much more 
generally than just in this two-dimensional case, for temporal variations in 
the rate of flow from the jet will be described in just this way. It should 
be noted that although these fluctuations are characterized by very small 
wave-numbers, their frequency (in time) may be quite high. 


lor k > 1, the second term in (8.1) becomes negligible, and we have 
ups ro? | Ae!hcos(B+klogr+}a)— 
0 | 
J he ) 1 
—ksin(B+ logr+iz)! dk, | : 
) webct (8.4) 
uy —r? |) Ae} sin(B+klogr + 3a) 
0 


+kcos(B+klogr+4x)} dk. 
B(k) is a rapidly changing function of k satisfying the inequality 
—-7<B<z7, 


and if we average over different flow realizations we would expect that 





sin {B(k,) + B(ka)} = cos{B(k,) + B(ks)} = 0, 
unless kj = ky. B will take all values equally probably in different flow 


realizations, so 


w= =dr| A714R2%)c2dk | - 
J 0 (3 t (8.5) 


uu, = V. J 


Thus u? = uz, and wv, and u, are uncorrelated, and (8.5) is obviously the 
two-dimensional analogue in polar coordinates of Phillips’s results (1.1) 
and (1.2). 

The scale of a motion characterized by k = 1 will be such that logr 
changes by about $7, that is r changes by a factor of the order of 5. Thus 
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k ='1 is a very small wave-number indeed, and virtually all the energy of 
the velocity fluctuations may be in the range of wave-numbers for which 
(8.5) is applicable. 

Although a detailed analysis of the actual case of three-dimensional 
fluctuations in the two-dimensional jet has not been found possible, some of 
its probable characteristics can be shown. 

It can readily be demonstrated that for irrotational flow in which mean 
quantities are invariant in the z-direction the cylindrical coordinate equi- 
valents to equations (3.8) are 


Ci 2 
1 2 ’ x | 
= (29° — ur) oS ee ames] 
Cc i \ 
\ 
é a —— a——si| te 
1 72 rn . 
—(43q°—u-) = 2u,u,t+r—Uu,u,. | 
a, \29 x iad | or * & J 
Our similarity conditions permit us to write 
P= rtp, ue = r-ly?, (8.7 


where p and 7@, v, are independent of r. Equations (8.6) then become 
iV; ] 





(8.8) 


ase 4 


In view of (1.8) and (8.5), the most likely solution to these equations is 

v2 = 4p’, 0,0, = 0. (8.9) 
It would seem probable then that there is no Reynolds shear stress 

across surfaces x = constant in the two-dimensional jet. This is a result of 

the behaviour g? « r~! and would not be possible with any other dependence 

on r, unless vu, = 0. It is in contrast to the situation in the wake, where 


we saw that a shear stress does occur. 


9, CONCLUDING NOTE 

This study, as far as it has gone, has shown that the results of Phillips’s 
simple model of irrotational flow outside a turbulent region are not greatly 
affected by introducing simple inhomogeneities into the boundary conditions. 
It has been demonstrated that these irrotational fluctuations are capable in 
some circumstances of producing a Reynolds shear stress, but the basic 
conclusion that the velocity component directed away from the boundary of 
the turbulence contains about half the total energy in the fluctuations seems 
to be firmly established. Still unresolved then, is the contradiction pointed 
out by Townsend (1956, p. 191) between this kind of theory and the experi- 
mental results which indicate that in jets the downstream component is 
larger than this cross-stream component. Although it is not altogether 
satisfactory to rely on the extensibility of the results of the two-dimensional 
analysis of § 8 to real three-dimensional cases, flow similar to that described by 
(8.2) is intuitively likely and (8.5) is very similar to the three-dimensional 
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result (1.1). It seems reasonable to infer that in three-dimensional cases 
there will also be two basic modes to the irrotational fluctuations: very low 
wave-number motions along the mean flow lines, and higher wave-number 
motions similar in structure to those in gravity waves. Figures 4 and 5 and 
equations (1.1) and (8.9) show that both modes give greater contributions to 
the cross-stream than to the down-stream motions. 


This paper describes work done for the Defence Research Board of 
Canada. 
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The intensity of Aeolian tones 
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SUMMARY 

The generation of Aeolian tones is interpreted in terms of the 
theory of aerodynamic sound. ‘To do this, the intensity of the 
radiated sound is expressed in terms of the fluctuations in force 
upon a moving rigid body using the approximation that, at low 
Mach numbers, these forces can be calculated assuming incom- 
pressibility of the flow. The fluctuations in lift and drag upon a 
circular cylinder at Reynolds numbers between 40 and 160 are cal- 
culated by integrating the fluctuations in momentum in the eddying 
wake, using the experimental data of Kovasznay (1949). It is 
found that the fluctuating lift per unit length is approximately 

f, = 9-38pU7d cos 2rnt, 
where p is the density of the fluid and ¢ the time, and that the 
magnitude of the fluctuations in drag is about 10°, of this. 

The axial length scale of these force fluctuations was found 
by an auxiliary experiment to be very large for Reynolds numbers 
below 100. When 100 < R < 160, the length scale is approxi- 
mately 17d, the transition apparently occurring as a result of the 
instability to three-dimensional disturbances of the laminar 
eddying wake. Using this datum, the mean square acoustic 
pressure generated by the motion about the cylinder at these 
Reynolds numbers is found to be 

ar 97 cos2g Ph GUS® 

pr) ~ 0:27 cos '— 
where @ is the angle between the direction of observation and 
the incident stream, / the length of the cylinder and a the velocity 
of sound in the medium. Experiments undertaken to test this 
result directly, using sound intensity measurements from a whirling 
arm apparatus in an accoustically quiet room, gave very good 
agreement. 

At higher Reynolds numbers, when the cylinder wake is 
turbulent, the theory leads to a similar expression for p? but with 
a smaller numerical factor. Analysis of previous experimental 
data on this basis gives good agreement, with a value of about 
0-037 for the numerical constant. In view of this, it is likely 
that some of the earlier data which had appeared to give a Mach 
number dependence of M*, not M®, have been misinterpreted. 


607 





608 O. M. Phillips 


Finally, mention is made of the conditions when the frequency 
of vortex generation is the same as the natural frequency of the 
wire, and the greatly increased intensity of the sound is seen to 
be the result of the greatly increased axial length scale of the 
force fluctuations. 

The relative simplicity of this phenomenon makes possible 
perhaps the least ambiguous confirmation yet found of some of 
the essential ideas in the theory of aerodynamic sound. 


1. INTRODUCTION 

\eolian tones were known in antiquity. ‘The Aeolian harp is an old 
musical toy which, when placed in the wind, produces these sounds in 
strange successions and combinations. They are heard in the singing of 
the wind through telegraph wires and its sighing through the twigs and 
branches of a leafless tree. It was not until 1878 that the first scientific 
study of them was made by Strouhal, whose apparatus consisted of a 
cylinder mounted upon a whirling arm so that it moved through the air 
at constant speed. Already by 1896, Rayleigh had recognized that the 
production of Aeolian tones is connected with the instability of the vortex 
sheets in the cylinder wake, and most significantly, that it is not essential 
that the cylinder should partake of the vibration to generate these tones, 
although the intensity is greatly enhanced by resonance between the natural 
frequency of the wire and the frequency of the vortex instability. He 
further perceived that the frequency of the note (expressed non-dimensionally 
as the Strouhal number S) depends upon the Reynolds number of the flow 
past the cylinder, and until recently, much of the later work on this 
phenomenon was connected with the determination of this dependence, 
and with the relation between the frequency of the note and that of the 
vortex instability in the wake. Relf (1921) showed that the fundamental 
frequency of the note and the vortex frequency on each side of the cylinder 
wake are equal over a range of Reynolds numbers which has since been 
extended to values between 50 and 10! by the work of Kovasznay (1949), 
Roshko (1953) and Gerrard (1955). Measurements of the intensity of the 
sound under various conditions were made by Holle (1938) and Gerrard, 
who measured its directional distribution also. Comparison will later be 
made between the results of these intensity determinations and the 
theoretical predictions of this paper. In 1896, Rayleigh wrote concerning 
this problem that ‘‘a dynamical theory has yet to be given”’, and to some 
extent this statement remains true. ‘The purpose of this paper is to interpret 
the properties of these tones in terms of our empirical knowledge of the 
flow pattern; the prediction of the latter still lies beyond the attainments 


of present theory. 

Rayleigh’s important observation that Aeolian tones can be produced 
without vibration of the cylinder suggests that the phenomenon can be 
interpreted within the framework of the theory of aerodynamic sound. 
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However, to include the resonant case also in our considerations, the 
definition of this term can be extended slightly. It has previously been 
restricted to sound generated by the action of aerodynamic forces upon 
fixed surfaces, but since in many instances, these forces will themselves 
produce vibrations of the surfaces, the term can conveniently be extended 
to include those cases in which oscillations of the surface may occur as a 
result of the instability of the basic flow. A distinction is therefore drawn 
in what seems the natural place, between such situations and those in 
which the vibrations are caused by some external mechanism unrelated to 
the flow (a violin bow, perhaps). In the latter class, it is the vibrations 
themselves that set up the fluctuating force in the medium, whereas in 
aerodynamic sound the vibrations may occur as a result of the force 
fluctuations from the flow. 


‘THE SOUND FIELD 

The two well-known papers by Lighthill (1952, 1954) which laid the 
foundations of the subject of aerodynamic sound, were concerned with the 
sound field generated by the unsteady motion of an unbounded fluid. 
More recently, Curle (1955), using the same technique as Lighthill, was 
able to account for the effect of solid boundaries to the flow and showed 
that the sound can be considered to be produced by three separate 
mechanisms. ‘The first is due to the unsteadiness in the fluid itself, and 
as Lighthill showed, is equivalent acoustically to a distribution of quadrupoles 
throughout the fluid. Curle pointed out that, in general, the presence of 
boundaries introduces additional dipole contributions to the sound field 
from the fluctuating stresses acting upon the fluid at the boundary and 
also perhaps source contributions from dilatation of the boundaries. ‘lhe 
general expression for the density field p(x) involves terms which represent 
these three effects separately, and can be written as 





ieee | Piteede 
p(X) — py = 47a? ax, ax, | 7 ily) « (y) T 
1 ofl 1 (1 A(pz,) 
$+ —— — (pv, v; + 3) dS ( )- —- | - dS, Ay), (2.1) 
4ra2 Ox; ] r (Pests + Pig) @OAY 47a‘, | r at 


where py and a, are the density and velocity of sound in the undisturbed 
medium, 7; and p;; are the velocity and stress tensor at a point y in the 


fluid and r = |x—y|. The tensor 7;; is given by 
T 5 = pv; tj + Pi a, p 9535 ( 
and the integrands in (2.1) are taken at the retarded time t—r/dp. 
This gives the fluctuations in density in the medium in terms of the 
properties of the flow, which are assumed to be specified a priori. In 
general, of course, the density field itself must be included in the specification 
of the flow, and in practice (2.1) is meaningful only to the first order in the 
Mach number in which the direct dynamical interaction between the sound 
field and the velocity field can be ignored. ‘To obtain this first order term, 
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it is legitimate to substitute into the right-hand side of (2.1) the values 
that would obtain if the fluid were incompressible, and this, incom- 
pressibility approximation’ has invariably been made in the applications 
of the theory up to the present time. With its aid, retaining the leading 
term in each of the integrals, it is found that the radiated part of the sound 
field from the motion in the neighbourhood of a finite closed body is given 


by 
t <a 0 : l “fe 
) -p~pyp ~~ — —~ | — T., dV(y)+ — -s | —(por;7;+9,;) dSj(y 
p(X) PO Anat x3 or ij ¢ (y) 47a? x2 Ara ek Pi;) ¢ ly) 
1 Po ov ’ 2 
_—P0 | Ziggy), (2.3) 
4razx | ot’ Ay) (49) 


Curle (1955) simplifies this expression when the boundaries are perfectly 
rigid perfect acoustic reflectors—and at rest, so that U; 0 over Si(y), 
but it will be shown here that the same simplification is possible in rather 
less restrictive circumstances. A sufficient condition for the vanishing of 
the source contribution is simply that the body be incompressible, so that 
its dilatation is zero. “Vhe Reynolds stress term in the dipole contribution 
also vanishes if the surface moves uniformly as that of a rigid body. For 
then vy) = u;+<«,;,w;V, on Sy), where wu; is the velocity of, and w, the 
angular velocity about, the centre of mass, and on substitution into this 


term of the integral and application of Gauss’s theorem, it is found to 


vanish. ‘This second condition clearly embraces the first and (2.3) reduces 
to 
1 xx, Fer, 1 wdF 
ij re i i >) 
o(X)— py ~ —_~ | ——~ dV(y)-— —— —_—., (2.4) 
j i 47a x® | dt? : tra’ x? dt 
where F | p;; @S; is the total force with which the body acts upon 


the fluid. Finally, when the orders in the Mach number .V/ of the two 
terms on the right of this equation are examined, it is found that the dipole 
contribution is of order ./* whereas the quadrupole term is of order MW), 
that is, of the same order as the quantities already neglected in the in- 
compressible approximation to the dipole term. So the quadrupole term 
too can be neglected and, when is small, 

1 x,dF, 
~ Fait a” 

In the generation of Aeolian tones, the fluctuating force is distributed 
along the length of the cylinder, and at low Reynolds numbers, if the flow 
were accurately two-dimensional, the force fluctuations would be in phase 


(2. 


wi 





p(X) Pa ) 


at all points. However, there are many reasons for believing that in a 
practical situation with a long cylinder (say, a long, thin wire) this will not 
be so, but that the phase of the force fluctuations will vary slowly, probably 
randomly, along its length. One reason is that the inherent small variations 
in the cylinder diameter will result in small changes in the frequency of 
vortex generation at different points which, in a large number of cycles, 
produce finite variations in phase. The coupling between the motions in 
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planes separated by a large number of diameters cannot be presumed to 
be so tight that such phase differences will not occur. Again, even with 
a perfectly uniform cylinder, if there is a three-dimensional instability of 
any kind, phase differences assuredly will be found, the extreme case being 
given when the flow behind the cylinder is turbulent. This question will 
be considered later in some detail, but if, for whatever reason, a random 
phase variation is present, then the sound intensity radiated from a 
cvlinder of length / is approximately 


VX; ] i’ Of (%) Of (So T >) 








Be ep ee es Bald te Fk Sl 7" 
(x) 1672pa3x! ot ct d 
: Foe : 
= i a Ze 
1677pa*x4 ( ‘i ct ) sal iia 


where f,(z) is the fluctuation in force per unit length at the axial position 2, 
s is the dimensionless correlation distance along the axis (SO that the actual 
distance is sd) and the subscripts to p and a have been omitted. In (2.6) 
it is further assumed that / < sd, as seems likely in many practical cases, 
so that end effects can be neglected. ‘he two principal quantities to be 
determined on the right of this expression are first the mean square value 
ot cf,/ct and second the axial scale s. The first is to be found from the 
experimental results of Kovasznay (1949) for the eddying motion in the 
wake at a Reynolds number in the regime of two-dimensional flow. ‘This 
restricts our consideration initially to Reynolds numbers below those for 
which the motion is turbulent, although the theory is extended in a later 
section to conditions beyond this range. ‘The second quantity, the corre- 
lation distance s, is the subject of a subsidiary experiment described in §5. 


3. "THE FLUCTUATIONS IN LIFT AND DRAG UPON A CIRCULAR CYLINDER 

‘The fluctuations in lift and drag have been calculated by integration of 
the acceleration field of the eddying wake in the directions normal and 
parallel to the direction of the oncoming stream, using the very extensive 
experimental data given by Kovasznay (1949) for a Reynolds number of 50. 
He expressed the velocity fluctuations in the direction of the oncoming 
stream, in the two-dimensional motion, as 

u(x, vy, t) = $,(x, v)cos 27[¢,(x) + nt] + 4o(x, y)cos 4z[,(x) + nt], (3.1) 

where the amplitude functions ¢, and 4, vary slowly with x (the downstream 
coordinate), but rapidly with y, i.e. across the stream. ¢, is an odd function 
of y, and ¢, an even function. Here ¢,(x) and ¢,(x) are not quite linear 
functions of x, in order to account for the slightly non-uniform spacing of 
the vortices, and the frequency m is constant in time. Kovasznay measured 
$,(x, y) and ¢,(x, y) across the wake at various stations downstream as far 
as x/d = 40, and, by smoothing and interpolating between his points, 
these functions have been found for values of y within the wake and for 
values of x between 0 and 40d. Interpolated values of the functions ¢,(.x) 
and ¢,(x) have also been obtained from his measurements. 
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The fluctuations in drag per unit length can be found directly by 


integration of (3.1) and use of the equation 
d f/f 
j= pwr u dxdy, 
Ia p dt | | 7 


where the integration is over a plane normal to the cylinder axis. ‘lhe 
calculation presents no difficulties concerning convergence, since ¢,(x, y) 
is negligibly small for values of y beyond the centres of the vortices. ‘lhe 
function ¢,(x, y) is odd in y, so that the integral of the first term in (3.1) 
vanishes, and the fluctuations in drag have a frequency just twice the 
fundamental n. However, the calculation of the fluctuations in lift per 
unit length f, is not so straightforward, since the velocity fluctuations normal 
to the oncoming stream were not measured. ‘They could perhaps be 
calculated from the continuity condition 


A A 
( 


u ov 
-—~+ = =J0, 
Cx rom V4 
ry OU 
which gives v= | — iy. 
| Pe 


/ x 
‘The difficulty with this method is that measurements of u were not made 
beyond the vortical region of the wake, and to make a correct allowance 
for the virtual mass of the irrotationally moving fluid would indeed be 
complicated. An alternative method was used which required information 
about the vorticity field alone. 
It is shown in the Appendix that 


|| edxdy = — || ww drdy, (3.2) 


where the vorticity w, which vanishes outside the wake, is given by 





Cw Ou tu 
Te ~~) ~_ = 
ey ox® = dy? 
“y ofu c2u 
Hence w= | ( <r in =) dy, (3.3) 
oe” a] 


the integration being commenced at the edge of the wake. ‘The integral 
{{«dxdy is calculated from (3.2), (3.3) and (3.1), and some of the details 
of this, including the devices used to ensure more rapid convergence, are 
described in the Appendix. 

The result is found to be 


M = p|| v dxdy = 0-48pUd?cos 2ant, (3.4) 
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where U is the mean stream velocity, and thus the fluctuating lift per unit 
length is 


. dM ee 2 
f= — = 0:38pU"d cos 2znt, (3.5) 
dt 
(taking a new time origin), since 2 = 0-13U/d at R= 50. ‘The fluctuation 
in drag per unit length is found by direct integration, and its magnitude is 
very much smaller. In fact, 


f, = 0-04p U"d cos 4rnt. (3.6) 


Just as the fluctuations in drag are determined by the term in (3.1) which 
is an even function of y, so the fluctuating lift is determined by the odd part 
in y, as the result (3.5) implies. For, from (3.2), only the even part of » 
is relevant, or, from (3.3), only the odd part of uw. 

The mean drag per unit length at a Reynolds number of 50 is about 
0-75pL/*d, so that our calculation shows that the magnitude of the maximum 
lift is about 50°,, of the mean drag. 

The value 0:38 for the coefficient in (3.5) can be shown to be a not 
unreasonable figure by constructing a very crude upper limit to the 
amplitude of the fluctuating lift in the following way. Suppose that vortex 
sheets of strength U are generated at each side of the cylinder which, because 
of their instability, curl up into discrete vortices like a Karman street. 
Suppose further that all the vorticity is gathered into these vortices, none 
being lost by diffusion. ‘The separation of the vortices is approximately 
!U n, which, at a Reynolds number of 50 and Strouhal number of 0-13, 
is equal to nearly 4d. ‘The strength of the discrete vortices is therefore 4+Ud, 
and the amplitude of the fluctuating circulation about a fixed contour in the 
irrotational fluid is 2Ud, so that the amplitude of the fluctuating lift is 
pl .2Ud = 202d. ‘That this is so much greater than the actual value can 
presumably be ascribed to the fact that the process of viscous diffusion 
would lead to a reduction of the coefficient. 

‘he accuracy of the calculated value (3.5) is admittedly not very high. 
The numerical differentiations reduce the accuracy, only partial com- 
pensation being provided by the subsequent integrations, and inspection 
of the working suggests that an accuracy of 40°, is probably the best that 
can be claimed. However, this is ample for the present purpose in view 
of the uncertainties of the acoustical measurements themselves. 


4. "THE EFFECT OF REYNOLDS NUMBER 
‘he estimate (3.5) was obtained from consideration of the flow pattern 
at R = 50, and it would be very satisfying if similarly detailed measurements 
were available at other Reynolds numbers in the range of two-dimensional 
eddying motion. Unfortunately, they are not, and some general argument 
will have to be invoked concerning the effect of viscosity upon the flow 
pattern. Such an argument would not be required to give results that are 
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accurate in detail, but only to give a good indication of the variation of the 


amplitude of the lift fluctuations with Reynolds number. 

Let us consider the mechanism of the development of the vortices in 
the wake. ‘The initial instability introduces small fluctuations in the 
y-position of the vortex layer on each side of the wake, and these result in 
differential convection velocities of the vorticity and the eventual rolling 
up of the layer. The development of the vortices is opposed by the ditfusive 
action of viscosity, and the maximum concentration is reached when the 
two effects balance. ‘Thereafter the vortices decay by viscous ditfusion 
between those of opposite sign, but it is not difficult to show that, at R = 50, 
this process hardly commences within the first 50 diameters (this being the 
part of the wake upon which the lift and drag calculations of §3 depend). 
For, in time ¢, the vorticity diffuses over a distance of order (vt)!*, so that 
at a point k diameters downstream it has diffused a distance of order 
(vkd/ U)!? = (k/ R)' 7d. When R = 50, this is of order d (still much less 
than the vortex spacing) at 50 diameters from the cylinder; the details 
of the motion in the wake beyond such points having little influence upon 
the fluctuating forces on the cylinder. For the present purposes, then, a 
decrease in the viscosity of the fluid (or an increase in R) has two conse- 
quences ; it modifies the frequency of the initial instability and also postpones 
the attainment of the balance between viscous and convective vorticity 
transfer, so that the vortices become more concentrated. 

Regarding the second effect, (3.2) shows that, in discussing the 
fiuctuating lift, we are concerned not so much with the concentration of 
the fluctuations in vorticity as with their integrated magnitude or with the 
total strength of the vortices. If we assume as an approximation that the 
lateral spacing of the vortices is materially unaltered by a decrease in 
viscosity, the Reynolds numbe: effect must lie in the variation of the strength 
of the vortices. ‘This is indeed a consequence of the first effect, the modifi- 
cation of the frequency of the initial instability and so of the spacing of the 
vortices. ‘The strength of the vortices can be assumed to be proportional 
to their spacing, U’/n approximately, so that as the viscosity varies, the vortex 
strength varies as S~!, since S=ndUl. ‘Vherefore the momentum 
fluctuations normal to the direction of the oncoming stream are proportional 
to S“!, and on differentiation /, is independent of S and so of R; consequently 


ff, <= 0:38pU7d cos 2znt, (4.1) 
tor 40 R 160. It follows immediately that 


Cf, 


; - 2-4o0U3S cos 2znt. (4.2) 


The Strouhal number S is well defined empirically as a function of Reynolds 
number (Kovasznay 1949, Gerrard 1955) so that this equation gives 
approximately the dependence on Reynolds number of the term ¢f; ct 
in the expression (2.6). 
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5. "THE CORRELATION DISTANCE $ 

Very little published information is available concerning the scale of 
the variations in phase of the fluctuations along the axis of the cylinder. 
Both Roshko (1953) and Gerrard (1955), making observations at Reynolds 
numbers greater than about 150, found slow spanwise variation in phase 
along the cylinder, but insufficient data were available to estimate its length 
scale. On the other hand, Kovasznay (1949) found no measurable differences 
at R = 50 over distances of the order of 20 diameters, so that presumably 
the length scale is very much greater than this. An estimate of s is necessary 
for the application of (2.6), and so to obtain this, and to clarify the situation 
if possible, a supplementary experiment was made in a water tank. 

A hollow tube of outside diameter } in. was stopped up at one end 
and a row of small holes drilled down one side at | in. centres. ‘These holes 
were then filled with plasticine which could be removed easily from two 
of them to give open holes separated by any desired number of diameters 
from | to 30. ‘The tube was filled with dye and lowered vertically into a 
tank of still water. The tube was drawn through the water with a given 
velocity, and the vortex patterns were observed at two different axial 
positions by the dye traces at the two free holes. End ettects were avoided 
by making no observations within 10 diameters of the stopped end of the 
tube or within 4 diameters of the free surface. ‘The differences in phase 
were measured by looking vertically downwards at the dye patterns behind 
the cylinder, and, by repeating the experiment a large number of times 
for points at a given separation, the correlation coefficient of the fluctuations 
in lift was estimated. By repeating this procedure with holes at different 
separations, the length scale of the lift fluctuations was determined 
approximately. 

The results of these observations can be summarized as follows. For 
Reynolds numbers below about 80, very little phase difference was observed 
even for the largest available separation of 30 diameters, so that under these 
conditions the axial scale is certainly much greater than this. However, 
when the Reynolds number was between 100 and 160, the phase differences 
were clearly observed. At a separation of about 30 diameters, the phases 
appeared to be almost uncorrelated. ‘The scale s was estimated to be 
between 15 and 20 diameters, and did not appear to depend critically 
upon Reynolds number or upon the presence or absence of small disturbances 
in the water before the passage of the cylinder. At Reynolds numbers 
between these two ranges, i.e. between 80 and 100, it was noticed that 
there was little difference in phase if the tank was allowed to settle for a 
minute or two between successive trials, but that random variations were 
set up with about the same scale as at higher Reynolds numbers if two 
trials were carried out in rapid succession or if the water was disturbed 
slightly prior to an attempt. 

These observations suggest that, at Reynolds numbers between about 
40 and 80, the eddying motion is truly two-dimensional, and apparently 
stabie to small disturbances in the oncoming stream. As the Reynolds 
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number increases to a value between 80 and 100, however, the laminar 
motion becomes unstable to three-dimensional disturbances with a length 
scale along the axis of about 17d, and the flow becomes only approximately 
two-dimensional. ‘This new type of motion persists until about R = 160, 
when yet further modes are excited and the flow eventually becomes 
turbulent. ‘These observations are consistent with those already mentioned 
for Reynolds numbers of 50 and 150. 

Returning to our basic problem, it seems that a reasonable value for s 
in (2.6) is about 17 for Reynolds numbers between 100 and 160. Below 
this range, a much larger value can be expected, and its magnitude would 
be determined, not by the flow conditions, but by the small imperfections 
and irregularities in the diameter of the wire, as has been explained in § 2. 
‘These can hardly be considered to be under the control of the experimenter, 
and reasonable a priori estimates of the magnitude of the radiated sound 
are difhcult without detailed examination of each particular wire. One 
further point arises from these observations. ‘The fluctuating lift and drag 
were calculated on the assumption of a two-dimensional flow, and some 
error will be introduced into (4.1) by the lack of two-dimensionality at 
values of R between 100 and 160. However, since the scale of the variations 
along the axis of the cylinder is very much greater than that in the plane 
normal to the axis, such errors are expected to be quite small. 


6. COMPARISON WITH EXPERIMENT 
I'he radiated intensity can be found immediately from (2.6) and (4.1). 
Neglecting the relatively small fluctuations in drag, the radiated sound is 
generated by the fluctuating lift, and so consists of a dipole field with its 
axis normal to the direction of the oncoming stream. From (4.2), 


ee 
Z* 21 °S¢ 
(7) 4 ‘ 


so that the total radiated intensity at points distant r from a cylinder of 
length / is 
a : old U®S* 
I(r) ~ 0:27 cos?9-——. , (6.1) 
See 
a®) 
where 4 is the angle between the oncoming stream and the direction of the 
line joining the point of observation to the cylinder, and s has been taken 
as 17. ‘he mean square acoustic pressure is 
pr) = pal(r), 
pldU*s? 


~ ()-27 cos?4 serena a 
avr 





(6.2) 


Some experiments were undertaken to test this result, particularly the 
variation of p? as the sixth power of the Mach number. ‘The measurements 


of sound intensity, which were made at the Acoustics Laboratory of the 
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National Physical Laboratory, Teddington, by the kind permission of the 
Superintendent of the Physics Division and Mr N. Fleming, were executed 
by Mr W. C. Copeland. I am greatly indebted to them for their generous 
and able assistance. The apparatus consisted of two thin wires, 60 cm in 
length and 0-0123 cm in diameter, held vertically with just sufficient tension 
to prevent sagging, between two arms mounted on a spinning shaft of 
diameter 3 in.—the same type of arrangement as that used by Strouhal 
and Gerrard. The supporting arms were made of 3 in. dural sheet 
tapered and streamlined to reduce disturbances. ‘The shaft was belt-driven 
from an electric motor, and, to reduce the background noise level, the 
motor and all the ball-bearings were enclosed in felt-lined boxes. ‘The 
apparatus was mounted on an open steel frame and, when in position for 
the measurements, was placed in the acoustically quiet room at the N.P.L. 
with a microphone 115 cm from the axis of the shaft in the plane bisecting 
the wires. ‘The apparatus was tested by running it with the wires removed, 
and the noise level was found to be negligible in the frequency range at 
which the measurements were made. 

‘The microphone signal was amplified and passed through a variable 
frequency narrow-band analyser, and displayed on the screen of a cathode- 
ray oscillograph. ‘The modulation with time, corresponding in the whirling 
arm apparatus to the variation with 4, followed closely that indicated by 
the cos?é factor in (6.2). For the intensity measurements, the peak values 
of p>, corresponding to orientations 4 = 0), were observed at a number of 
frequencies within the narrow spectrum of the sound generated, from 
which the total radiation was calculated. ‘The speed of the wires was 
altered by variation of their radial position and of the speed of the motor. 
The results of these measurements are shown as the solid points in figure 1, 
together with the predicted relation (6.2), and the agreement is seen to be 
remarkably good. However, the constant in (6.2) is reliable only to within 
a factor of two, this range of uncertainty being approximately the same as 
the experimental uncertainty shown by the scatter of the points. For 
values of U(/dS?/r?)!® less than about 140, corresponding in these measure- 
ments to Reynolds numbers less than about 100, an intermittent signal of 
very much higher intensity was sometimes recorded, which was consistent 
with the establishment of the regime of truly two-dimensional flow over 
most of the length of the wires and the consequent enormous increase in 
the value of s. However, under the experimental conditions, it was not 
possible to produce a sufficiently stable signal to make systematic observa- 
tions, and indeed, even had this been possible, there remains the strong 
suspicion that the results would not be generally applicable, but would 
depend upon the imperfections of the particular wire used. 

These results offer good evidence for the validity of the expression (6.2) 
in the Reynolds number range from 100 to 160. ‘The variation of the 
sound intensity as the sixth power of the velocity seems to be confirmed 
well, and in the next section it will be shown that this remains true when 
the flow is turbulent at higher Reynolds numbers. 
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7. ‘TURBULENT FLOW 

Hitherto attention has been concentrated upon conditions at Reynolds 
numbers below that for the onset of turbulent flow. As R increases beyond 
this range, intermittent bursts of turbulence occur in the flow, as described 
by Roshko (1953), until at about R = 300 the wake behind the cylinder 
is fully turbulent. The expressions given for the fluctuating lift and axial 
length scale are no longer relevant but we can still predict their order of 
magnitude and dependence on Reynolds number and obtain, for these 
new conditions, an expression for the intensity of the radiated sound which 
is arbitrary only to the extent of an unknown constant factor. 

It is clear that the presence of a fluctuating lift depends upon the 
existence of large-scale axial eddies of alternating sign near the cylinder. 
Some measurements have recently been made by Mr I. C. 'T. Nisbett at 
the Cavendish Laboratory in the turbulent flow near the cylinder at 
Reynolds numbers of order 5 x 10°. He has established the presence of 
a surprisingly distinct large-scale structure, similar in many ways to that 


observed in the two-dimensional eddying flow at lower Reynolds numbers. 
These large eddies are approximately parallel to the cylinder axis, and 


contain an appreciable amount (about half) of the turbulent energy, the 
remainder being in the irregular smaller-scale fluctuations. ‘hey appear 
to be somewhat contorted in the axial direction, since simultaneous velocity 
correlations between similar points in planes normal to the axis are alter- 
nately positive and negative as the axial separation increases up to about 


?() ‘The integral length scale of these eddies along the axis is 


diameters. 
about three diameters. ‘These observations, in the light of our previous 
discussion and of such formal expressions as (3.2), suggest that the 
Huctuations in lift are still dominant even at these Reynolds numbers, and 
that their length scale is about 3d, although we can make no a priori statement 
ibout their magnitude. 

What of the variation of the magnitude of the fluctuating lift with 
Reynolds number? Here, again, the general arguments of §4 still apply 


vith equal force, and we have 


(f7)'? = constant pl*d. (7.1) 
Some further support to this relation is attorded by the principle of 
Reynolds number similarity (Townsend 1956). This principle expresses 


the fact that the large-scale components of the motion are not influenced 
by the fluid viscosity if it is sufficiently small, i.e. by the Reynolds number 
provided it is sufficiently large. We have seen that the fluctuating lift is 
a consequence of the large eddy motion in the wake, so that (7.1), which 
is correct dimensionally, is also independent of the (large) Reynolds number. 
Whether or not the principle of Reynolds number similarity holds in this 
situation tor values of R as low as 300 is perhaps open to question, but for 


‘higher values it can certainly be invoked to support (7.1). Using this 
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relation, then, we have from (2.6) that the radiated intensity from a cylinder 
of length / is 
oq plaUsS? = 
I(r) =. COS*U aa 9 (7.2) 
a®r* 
which differs from (6.1) only in the constant factor which is now represented 
by x. Again, as in (6.2), the mean square acoustic pressure is given by 
_ oq prldU®S? a 
pr) = «cos*¢ —_.—_, (7.3) 
a*r* 
where the Strouhal number S involves the frequency of ‘shedding’ of the 
large eddies, or the fundamental frequency of the sound generated. ‘lhe 
constant x can be expected to be considerably smaller than 0-27, because 
of the reduced length scale of the fluctuations in lift per unit length and 
also possibly a reduced intensity arising from the decreased coherence in 
the turbulent flow 

Detailed measurements of the sound intensity radiated at Reynolds 
numbers for which the flow is turbulent have been made by Holle (19358) 
and by Gerrard (1955). Now, Gerrard’s curves of intensity as a function 
of velocity at constant Reynolds number (obtained by keeping the product 
Ud tixed) appear to show that p? x U4, but this behaviour is inconsistent 
with the l° variation suggested by a rough application of the theory (as 
he himself points out) and his results have hitherto presented something 
of an anomaly. However, (7.3) shows that, with Ud fixed, we should 
expect that p? « U>, so that the anomaly is not as serious as might appear 
at first. A possible cause of the discrepancy lies in Gerrard’s having taken 
p? « Pr? when / > d in the analysis of his data, which implies that the 
Huctuations in lift are in phase along the entire length of the cylinder, over 
many hundreds of diameters. In turbulent flow, the case for assuming a 
random phase variation along the cylinder (leading to intensities propor- 
tional to the first power of the length as in (7.2)) is even stronger than at 
the lower Reynolds numbers discussed in the previous sections, yet 
independent confirmation that this is so even under these latter conditions 
has been given by the agreement between the absolute values of the 
measured and predicted sound intensities shown in the upper curve of 
figure 1*. Nisbett’s observations on the flow at high Reynolds numbers 
provide yet further strong evidence for random phase variations, so that 
it seems inevitable that, in fact, p? « /7? in the radiated field under the 
circumstances of Ge-rard’s measurements. 

In his figure +, some of Gerrard’s measured values of p? under a variety 
of conditions are plotted as a series of curves together with the relevant 
experimental parameters. Assuming that in fact P? x Il r*, these points 
have been replotted as a function of velocity and Reynolds number in 
the lower curve of figure 1 of this paper. ‘This diagram also shows a 


If, in (6.2), sd is replaced by the length /, the expression would give values 


tbout 250 times those measured in the experiments of § 5. 
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representative sample of Holle’s (1938) measurements”, and these agree well 
with Gerrard’s replotted points. ‘They all lie about a single straight line 
passing through the origin, as predicted by equation (7.3), and cover a 
fairly wide range of conditions. ‘The slope of the line is less than under 
conditions of non-turbulent flow and corresponds to a value of «x of 
approximately 0-037. The scatter of the points is not inconsiderable, 
but the linearity of the relation seems to be confirmed well, and provides 
good support for the result (7.3) which is based upon the ideas of the theory 
of aerodynamic sound. 
8. SOME FINAL REMARKS 

The satisfactory agreement between these various experimental results 
and the theory gives confidence that the latter is firmly founded. One 
factor ignored in our considerations, however, is the possibility of vibration 
of the wire. However it is not difficult to show that the amplitude of the 
vibrations is indeed small compared with the wire diameter provided the 
tension is such that the natural frequency of the wire is very far from the 
Vortex frequency, and the density of the medium is much less than that 
of the wire. Under these conditions, which have been satisfied in the 
experiments mentioned here, the vibrations of the wire can reasonably 
be neglected. On the other hand, if the frequency of vortex shedding 
coincides with the fundamental frequency of the wire, or with one of its 
low harmonics, then resonance will lead to vibrations of finite amplitude. 
Rayleigh’s observation that the acoustic intensity is greatly enhanced under 
these conditions was mentioned earlier in this paper, and there arises the 
question of accounting for this within the framework of the present theory. 

Now, the two basic variables which determine the intensity of the 
radiated sound are the magnitude of the fluctuations in lift, and the 
correlation distance s. The first quantity will clearly be modified by the 
presence of vibrations, both through the introduction into equations such 
as (3.2) of an additional term describing the motion of the boundary and 
through the modification of the flow pattern, but since the maximum 
velocity of the vibrating wire could not be expected to be much greater 
than the wind velocity, these could not reasonably be expected to contribute 
a factor of more than about 5 to the intensity of the sound radiation. A much 
more powerful contribution will be given by the increase in the value of s 
brought about by the equalization of the phase of the fluctuating lift along 
the wire. The vibrations, which may be in phase over the whole length 
and are at the same frequency as the vortex ‘shedding’, clearly provide a 
coupling between the motions in the fluid at different axial positions that is 
far more effective than any hydrodynamic coupling. ‘The phase of the 
large eddies loses its randomness along the cylinder, and is determined by 
the phase of the vibrations, so that for the fundamental the correlation 
distance sd now becomes the length /, which can represent an increase of 

* Incidentally, Holle found experimentally that p oc l/r? approximately, as in 
(7.3), and that the velocity dependence was as U’, nearly. But a reasonably fit is 
still obtained by plotting (p*)!® against U as in figure 1. 
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many hundreds of times. ‘lhe explanation of Rayleigh’s observation seems 
to lie in this effect rather than in any increase in the magnitude of /,. 
The value of this work may go beyond the interpretation and under- 
standing of the phenomenon of Aeolian tones. ‘The relative simplicity ot 
the fluid motion makes possible what is perhaps the least ambiguous 
confirmation yet found of some of the essential ideas in the theory ot 
aerodynamic sound, and gives even greater confidence in their application 


to more complex problems. 


APPENDIX. FLUCTUATIONS IN WAKE MOMENTUM 
From the vector identity 
x xw = V(x.v)—(v.V)x—(x.V)v, (A.1) 
it can readily be shown that 


v dx (x xw) dx+} |xx(dS xv), (A.2) 


vhere the integration is throughout the region bounded by the closed 


surface S. We apply this to the two-dimensional flow about a cylinde: 
centred on the s-axis, considering the volume bounded by the planes z = 0 
and : d, say. ‘Vhe cylinder is considered to be fixed, and the surtace of 
integration reduces to the two planes, so that (A.2) becomes 


- 4 
v dxd\ d {X <@) dxd\ ‘ i X (dS ; v) | 


| x x (dS » v) | . 


In the two-dimensional motion the fluctuations in velocity and vorticity 
are V = (u,z7,0) and w = (0,0,w), and on substituting these into the last 
expression, we obtain for the fluctuations in momentum in the y-direction 


normal to the oncoming stream, 
d | r dxd\ hd 1} Ww dxd\ t ld | ri dxdy, 
since dS = (0,0, —dxvdy) on s = d. Therefore 


v dxdy | | wx dxdy, (A.3) 


and similarly u dxdy | wy dxdy. (A.4) 


where the limits are understood to be taken in the sense 


lim( lim} | @ dyad). 
In the eddying wake behind a cylinder (of unit diameter) we shall 
neglect fluctuations in the vorticity of the boundary layer, since Kovasznay’s 


measurements indicate that the unsteady motion is small within the first 
few diameters of the centre of the cylinder. The fluctuations in vorticity 
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are found from equation (3.3), and on considering only that part of uw in (3.1) 
that is an odd function in y, 


24 _ (dd, 
=<>> 


oD) 


x? oxz 47°6,[C(x)]? > cos 2z[C,(x) + nt] — 


2») 


—< tn Ci (x) + 27d, C7 (x) > sin 27[C,(x) + nt], 


= A(x, v)cos 27[f,(«) + nt] — B(x, vy)sin 27[¢,(x) + nt], 
say. Therefore 


od , 
a= (: sc = )eos Zales) + nt] - Bsin 2z[C,(x) 7" nt}, 


where a(x, V) A(x, y) dy, Ax,y)= | B(x, y) dy, 

the upper limit being taken at y = 5 since it was found that the vorticity 
was zero along this line over the length of the wake to be considered. ‘The 
integration over v in (A.3) is most conveniently performed first; since ¢, is 


an odd function of vy, we have 


w dy = y,(x)cos 27[C,(x) + nt] — xo(x)sin 2z[C,(x) + nt], 
where Yi (4) | ady, Xo(*) B dy, 
and so | 
|| vdxdy =< | wxy,(x)cos27¢,(x)dx—| wxy,.(x)sin 27 C(x) dx >cos 2znt 
— 4) v ¥,(x)sin 27 C,(x) dx 4 | Vv xo(v)cos 27 C,(x) dx > sin 27nt. 


An integration by parts makes these integrals more rapidly convergent ; 


— || vdxdy = cos2znt | [6,(x)sin 27 C,(x) + 4,(x)cos 27 G(x)]-dx 4 


0 


+sin2rnt | [6,(x)cos 27 C(x) + 4,(x)sin 27 C,(x)] dx, (A.5) 





where 0,(x) = >> 


with a similar expression for 6,(x). 

The expression (A.5) is that from which the fluctuations in wake 
momentum were calculated. Kovasznay’s (1949) experimental results for 
¢,(x, y) were smoothed and interpolations were made so that a map of the 
function was obtained for the domain —5 < y <5, 0 <x < 40. The 
functions x(x, v) and §(x,v) were found numerically, and y,(x) and y,(x) 
calculated by integration. 6,(xv) and 6,(x) were thence obtained, and the 
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final result emerged from the integration of (A.5). It was that 
|| « dxdy = 0-48 Ud?cos 27nt, (A.6) 


giving the fluctuating lift as 

f, = 0-38pl 2d cos 2rnt, (A.7) 
since at this Reynolds number n = 0-13U/d. ‘These are the results quoted 
in (3.4) and (3.5). 
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SUMMARY 

This paper considers the flow due to a jet spreading out over 
a plane surface, either radially or in two dimensions. Solutions 
of the boundary layer equations are sought, according to which 
the form of the velocity distribution across the jet does not vary 
along its length. For laminar flow, such a similarity solution is 
obtained explicitly. For turbulent flow, an eddy viscosity is 
introduced, and it is eventually seen that complete similarity is 
not attainable, but that confident predictions can nevertheless be 
made about the nature of the velocity distribution and the rate of 
growth of the wall jet. 


1. INTRODUCTION 


When a jet of air strikes a surface at right angles and spreads out radially 
over it, it forms what will be termed a wall jet. Such a flow is produced by 
the downwards-directed jet from a vertical-take-off aircraft spreading out 
over the ground, or by a jet of water from a tap falling into a partly full sink 
and spreading out over the bottom. It should be noted that the flow due 
to a jet of water falling into an empty sink is in a quite different category, as 
then the condition at the free surface is that the pressure is constant; for 
a wall jet, as for a free jet, the corresponding condition is that the radial 
velocity component falls to zero at the outer edge of the jet. 

Plane wall jets may also arise. If the water levels are different in two 
sections of a canal, separated by a sluice, and if the sluice is slightly raised, 
the flow into the section with the lower water level will take the form of a 
wall jet. Alternatively, if a plane jet impinges on a fixed plate parallel to 
the direction of flow, a wall jet will be produced on each side of the plate. 

The theory of such wall jets, radial or plane, laminar or turbulent, forms 
the subject of this paper. It does not appear to have received attention 
previously. ‘The appropriate boundary layer equations are set up, and a 
search is made for a similarity solution, in which the form of the velocity 
distribution across the jet does not vary along its length. This is a familiar 
procedure in boundary layer theory. ‘Two similarity exponents, giving the 
variation with distance of the maximum velocity and the jet width, have 
to be determined, and one relation between them is obtained from the 
boundary layer equations themselves. In the corresponding problem for 
a free jet, a second relation is deduced from the constancy of momentum 
flux, and for a boundary layer the second relation is even simpler in form, 
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since the velocity must vary in the same manner as does the known external 
flow. For a wall jet, neither of the principles can be applied. It is shown 
in §2 that the solution depends on an eigenvalue problem, which for laminar 
flow is satisfactorily solved in §§2 and 3. The calculation of the jet profile 
is carried out in §4, and it is found possible to perform the integrations 
analytically. The same velocity distribution applies for both radial and 
plane wall jets, a fortunate circumstance which occurs in all the cases 
considered in this paper, and greatly shortens the work. 

In most practical examples, the wall jet will be turbulent. The attempts 
in §§5-9 at predicting its form are all based on the idea of replacing the 
molecular viscosity by an effective eddy viscosity, and requiring this to 
vary in a plausible manner. ‘The first attempt, in §5, makes use of the 
hypothesis for free turbulent flow due to Prandtl (1942), according to which 
the eddy viscosity is constant across the breadth of the jet. The velocity 
profile thus derived is shown to be identical with that in the laminar case, 
and further implications of this result are noted. Near the wall, this solution 
is not in accord with preliminary experiments (Bakke 1956), and in §6 
a second attempt is made, with an eddy viscosity as required to satisfy the 
law due to Blasius (1913) for flow ina pipe. The solution again is calculated 
without difficulty, and might be expected to be accurate in the region near 
the wall. A more realistic procedure is to postulate an eddy viscosity based 
on Blasius’s law near the wall and Prandtl’s hypothesis further out. This 
implies that complete similarity is no longer possible. Also, the eigenvalue 
problem to determine the similarity exponents, which in the earlier attempts 
was solved by a simple extension of the analysis for a laminar wall jet, has 
to be considered afresh. In spite of these difficulties, in §§7—9 an analysis 
is carried out which results in definite predictions as to the velocity dis- 
tribution and the values of the similarity exponents, as slowly varying 
functions of the Reynolds number. As in all problems involving Prandtl’s 
hypothesis, the value of a single constant has to be assigned from experimental 
observations. The final results are put forward with some confidence in the 
belief that they do correctly predict the main features of an actual turbulent 
wall jet. ‘ 

In what follows, the theory is developed for the case of a radial jet, as 
being of the greater practical interest. The modifications required to make 
the analysis applicable to a plane wall jet are only slight, and are noted at 
the end of each section. 


2. SIMILARITY SOLUTIONS 
We now consider the equations governing a radial laminar jet flowing 
over a plane wall. On the boundary layer approximation, the pressure is 
everywhere uniform, and the momentum equation is 


u— +7 = Vx) (2.1) 
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where x and y denote distances along and normal to the wall, x being 
measured from the jet axis, u and v the corresponding velocity components, 
and vy the kinematic viscosity. The equation of continuity is 
O(xu) (xv) 
eater ie (xv) ame 0, (2.2) 
Ox oy 


and hence there is a Stokes stream function i satisfying 


9 9 | 
ow ow 

xu= —, xwv=- = (2.3) 
Cy Ox 


The boundary conditions are 
“w=o=0 at y=0, u>Q asy> oo. (2.4) 
Consider the possibility that there shall be a similarity solution of these 
equations, with wu oc x“, and the jet thickness 6 oc x®, The two sides of (2.1) 
will vary with x in the same manner if 
a+2b=1. (2.5) 
in other boundary layer problems, a second relation between a and 5 is 
obtained from the general nature of the flow. ‘Thus for'a boundary layer 
in astream U = cx’, it is at once clear that a = m, and for a free jet, constancy 
of momentum flux gives the second relation. For a wall jet, there is no help 
from the boundary conditions, momentum is not conserved, and there is 
no obvious quantity to consider instead. Indeed % = 0 satisfies both 
equations and boundary conditions. ‘This suggests that only for particular 
values of a and 6 will non-trivial solutions arise. 
For convenience, we introduce non-dimensional variables by writing 
u=Uu, v=U0, x=vx/U, y=vy/U, p=ryp/U, (2.6) 


where U is a constant velocity. Using (2.5), we write 


where n = (2—b)yx-’. 
Then ui = (2—b)x*f'(n), 
and equation (2.1) becomes 


f" +if" +f? = 0, (2.7) 


where a = (2b—1)/(2—5). (2.8) 
The boundary conditions (2.4) require that 
f0)=F'(0)=0,  f'(~) =0. (2.9) 


For any «, there are solutions of (2.7) with f'(«)=0. We now 
investigate the possibility that, for such a solution, f and f’ shail both be 
zero at some finite y, so that (2.9) may be satisfied, with a suitably chosen 
zero of 7. Integrating (2.7) between the limits 7 and 0, we obtain 


f' + ff’ —(a—1)g = 0, (2.10) 


where a(n) = | f? dy. 
yy 
2T2 
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Note that g is certainly positive. Multiplying (2.10) by f’ and integrating 
again, we have 


Lf?—fe+(z—2) | f'g dy =0. (2.11) 


If, at the first zero of f’ to be encountered as 7 decreases, f is also zero, then 
since the integral in (2.11) is certainly positive, this equation requires that 
a= 2. Equations (2.5) and (2.8) now show that 

a= —3, b= 2. (2.12) 
Unique values have thus been found for the similarity exponents. The 
only assumption in the argument, that f’ is nowhere negative, 1.e. that 
there is no reversed flow, is clearly acceptable on physical grounds. 

For a plane wall jet, the equation of continuity cu/dx+dév/dy = 0 is 
satisfied by introducing a stream function y% such that u = cyp/dy, 
v= —dp/dx. On writing = %/v = x'-f(n), where 7» =(1—5)x-’, we 
obtain (2.7) precisely as before, where now « = (2b—1)/(1—6). The 
proof that « = 2 is unchanged, and hence, for a plane wall jet, 

a=-}, b = 3. (2.13) 
The equations being identical, the velocity distributions across the jet are 
the same in the radial and plane cases. 

Before proceeding further, we may note some important implications 
of (2.10). If f and f” are to vanish simultaneously, with f’ > 0, this 
equation shows that «= 1. ‘These are precisely the conditions which 
hold for a free jet (where du/ocy = 0 on the centre line y = 0), and (2.7) 
with « = 1 is indeed the equation for the velocity distribution in a plane: 
free jet, found by Schlichting (1933) and Bickley (1937). The above 
arguments show that this same velocity distribution applies also to the 
similar state of a radial free jet, as would occur when a radial wall jet had 
passed beyond the edge of a finite circular plate. For here, too, the equation 
must be of the same form as (2.7). The similarity exponents are found 
from (2.8) to be a= —1,6=1. These results for a radial free jet are in 
agreement with those of Squire (1955). A practical example of a radial 
free jet occurs, approximately, in the circular grilling burner of a gas cooker. 
The fact that, for a free jet, (2.10) leads to the same values of the similarity 
exponents as does a consideration of momentum flux, suggests that it may, 
after all, be possible to deduce the values for a wall jet by analogous 
considerations. ‘This is indeed the case, as will be shown in the next section. 


3. EXTERIOR MOMENTUM FLUX 
Our aim is to deduce an integral relation for a radial laminar wall jet, 
from the basic equations (2.1) to (2.4). Multiply (2. 1) by x and integrate 
with respect to y between the limits y and ©, using the condition that 
>Oasy— oo. By (2.2), 


(? ou ro O(xu) 





} xv>- dy = —xvut+ | u—— dy. 
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and so we obtain 

ee. ae 

= J, xu? dy — xvu + vx 5 0. (3.1) 

Multiplying by xu, and integrating with respect to y between the limits 0 

and o, we then have 
i [~ ‘ent [~ xu ay} dy — 4 sin % xu? ay} dy— 
GHiq” tte "i §6 fy Oe Use 

- | xvu? dy + : vata | =0. (3.2) 

0 


~ 0 


From the continuity equation, the second term of (3.2) is 


ro @ ro 
E | xu? dy | + | xvu2 dy. 
( 


) 40 


“9 

Now at y = 0, u =v = 0, and so (3.2) reduces to 

qr f pe ) 

—| xu: | xu dy> dy = 0, 

se sg ot te 
or ‘a Ew ) 

r= | xus | xu? dy ( dy = constant. (3.3) 
Yo lly 


Physically, this equation may perhaps be interpreted as saying that the 
flux of exterior momentum flux is constant, but this is hardly a familiar 
concept. The close similarity between the types of expression in equations 
(2.11) and (3.2) may be noted. In the case of a similarity solution, (3.3) 
shows at once that 3a+2b+2=0, (3.4) 


and, taken together with (2.5), this leads to the values of a and 6 already 
given in (2.12). If there is reversed flow, it is possible that F is zero, and 
in this case the deduction (3.4) cannot be made. So here, as in §2, there 
is the mathematical possibility of solutions with reversed flow, for some 
other values of the similarity exponents. 

The arguments leading to the constancy of F are not dependent on the 
velocity profiles being similar. The similarity solution requires an infinite 
velocity and zero jet width at x = 0, so in an actual wall jet there must 
be a region where the velocity distribution differs from its final similar 
state. Equation (3.3) gives valuable information about the state of affairs 
in this region. 

We may note that a rough estimate of the magnitude of F may be 
deduced from a knowledge of conditions in the impinging free jet, by 
writing a 

a!) rc ro 
F= | ‘id | xu? ay} dy = u* | ul | xU ay} dy 
| J ‘ ‘ 


~ 0 "7 


ll 
trol 
tJ 

7 
© es. 
ar 
| 
Ia 
cs 
™ 8 
& 
= 
& 
ue" 
a 
<= 


© 2 
= 1u{ xu ay} ‘ (3.5) 








630 M. B. Glauert 


where u* is a velocity typical of that in the jet. Thus 

F = }(typical velocity) x (volume flow per radian)’. 
Neither the volume flow nor the magnitude of the fluid velocity will alter 
greatly when the jet is deflected on striking the wall, so that the value of 
(3.5), estimated from conditions in the impinging jet, provides a measure 
of F in the wall jet. 

For plane flow, the analysis of this section is unchanged, except that 
the x’s are omitted and the final result, 

F= | us | wu? dy> dy = constant, (3.6) 
/0 L/y J 
leads to 
3a+2b = 0, (3.7) 
again agreeing with the previously obtained result. 

If the flow is turbulent, momentum flux arguments still hold for a free 
jet, but for a wall jet F will remain constant only under certain special 
assumptions about the nature of the turbulent stresses. ‘These matters 
will be discussed in later sections. 


4. VELOCITY DISTRIBUTION 
The expression for the stream function given in §2, with the value 
b = 2 found in (2.12), was = (v5x3/U)"4f(y), where » = 3(v3/Ux5)"4y, 


and hence u = }(v3/Ux*)!*f"(y). The equation for f was shown to be 
fae gee (4.1) 
with boundary conditions f(0) = f’(0) = 0, f’(«) = 0. 
Multiply the equation by f. ‘The first term may now be integrated by 
parts, and the remaining two form a perfect differential. Hence 


iia =S. (4.2) 
since f'(#) = 0. Multiply by f-**, and integrate again. Thus 
ff’ +3f?* = constant. (4.3) 


Now if f(y) is a solution of (4.1), so also is f,(y) = Afo(An), for any 
constant A, and f; satisfies the boundary conditions if f, does. But the 
effect of selecting f, instead of f, as the solution for substitution into the 
expressions given above for %, 7» and u is precisely the same as that of 
changing the value of the arbitrary constant velocity U to U/A*. Hence, 
without loss of generality, we may select the solution with f(.) = 1, and 
take the constant in (4.3) to have the value 3. Write f = g?; then f’ = 2gg’ 
and (4.3) becomes 


> 


Evie 3 
al ha 3 (4.4) 
which gives on integration 
(1l+¢+ g*) - \ 3g 
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The values of 7, f and f’ corresponding to given values of g may now 
be tabulated. The variations with » of f and f’ are illustrated in figure 1. 
The main properties of the velocity distribution are easily calculated 
directly. Thus, the shearing stress at the wall is related to f’(0) = , and 
the maximum value of f’ is 2~°% = 0-315, occurring at the point at which 
f =27-* = 0-397. 











Pe 4025 
/ tial 
: 4 1 1 I L 
o) 2 3 4 : 6 


Figure 1. Laminar wall jet. Variation of mass flux (f) and velocity (f’) with distance 
from the wall (7). 


The arbitrary velocity U may be eliminated by introducing the quantity 
F of §3. It is easily verified that 


pi el pee 2 6 
=), 04) Paypa= ge (4.6) 
Hence we may write 
0 Fyx3\14 . s 
b= (==) f(a), | 
a 

1I5F\12 ,, | - 
7" (8 f(t), } (4.7) 





The skin-friction 7, is given by 


Tm (cu _ (125F3\"4 
w= +(5) = (teat) (48) 
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For a plane wall jet, the expressions corresponding to (4.6), (4.7) and 
(4.8) are 


b = (40Fvx)!4f(n), 


u= (3) FO. (4.9) 


$Y 


T 1 /125F3\14 
p 7 3( =") : 


5. PRANDTL’S HYPOTHESIS 


Came 


In most practical situations, the wall jet will be turbulent, and we now 
investigate how the foregoing analysis must be modified to apply in this 
case. If we introduce the concept of an eddy viscosity e, the boundary 
layer equation (2.1) becomes 


ous ou é [ cu 
Cr Ue SS SER hs 
Cx oV cy oyv 


where u and vw now denote the components of the mean velocity. The 
continuity equation (2.2) and the boundary conditions (2.4) are unchanged. 
Some assumption must be made about the behaviour of e. The simplest 
and most convenient one, which has proved satisfactory in describing free 
turbulent boundary layer flows, is the hypothesis due to Prandtl (1942), 
that the eddy viscosity is constant across the layer, and is proportional to 
the product of the maximum mean velocity and the width of the layer. 
A wall jet would appear to be half-way between a free flow and a wall flow, 
so it is not unreasonable to adopt Prandtl’s hypothesis as a first attempt. 
If there is to be a similarity solution of (5.1), with u oc x7 and 6 « x”, 
we then require that € oc x“+°, and the two sides of (5.1) will vary with x 
in the same manner if 
b= 1, (5.2) 
In the arguments of §3 relating to exterior momentum flux, v has to be 
replaced by «, but since this is a function of x only, the analysis applies 
unchanged. Once again, F = constant, whether or not there is similarity. 
For a similarity solution, (3.4) holds as before, showing that 


i i 5 
a= “Tee (5.3) 
To study the velocity distribution, introduce the non-dimensional variables 
9) : : a a 
(2.6), and let e = Nx lsy, (5.4) 
where JA is a constant, to be determined later. Then if we write 
24 2y 








wa (5.5) 
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p= 2m), n= x 
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(5.1) shows that f satisfies f”+/f"+2f'? = 0, with boundary conditions 
(0) = f'(0) = 0, f’(«) = 0, precisely as for the laminar case. The solution 
of §4 is still applicable, and thus, on Prandtl’s hypothesis, the velocity 
profile is identical in the laminar and turbulent cases. This result follows 
at once from the integral arguments of § 2, as soon as it is established that the 
equation is of the same basic type in the two cases, and it applies equally 
well to a plane turbulent wall jet. 

For a plane free jet, it was shown by Gértler (1942) that, on the basis 
of Prandtl’s hypothesis, the velocity profile of a turbulent jet is the same 
as that of a laminar one. This again is an immediate deduction from the 
arguments of §2. Further, a radial free jet, as described in §2, would also 
have this same velocity profile in either laminar or turbulent flow. For all 
these jets, the use of Prandtl’s hypothesis is well justified. 

It would indeed be satisfactory if this simple approach sufficed to 
describe a turbulent wall jet. However, it is found by experiment (Bakke 
1956) that the velocity gradient near the wall is much greater than in the 
velocity distribution of figure 1. Consequently, we snall not proceed to 
set out further details of the solution based on Prandtl’s hypothesis. We 
merely note that, in this method, the only undetermined quantity is the 
constant A introduced in (5.4), and this would be chosen to give the best 
agreement with experiment, as in other problems where Prandtl’s hypothesis 
is used. 


6. EFFECT OF THE WALL 

Once it is conceded that there is large velocity gradient near the wall, 
the idea of a free turbulent flow is seen to be an inadequate one. It becomes 
clear that the friction at the wall is a decisive agency, and must be given a 
prominent place in any theory. We shall retain the concept of an eddy 
viscosity, but shall require it to vary in a suitable manner across the breadth 
of the wall jet. 

We may make use of the empirical formula due to Blasius (1913), based 
on a study of turbulent pipe flow, in the form 


- orf v \wA , 
tT, = 0:0225pL (Zz) : (6.1) 


where U is the maximum velocity and a the radius of the pipe. The formula 
is often written with the mean velocity and the pipe diameter in place 
of U and a; this merely involves a modification of the numerical coefficient. 
Since 7) is governed by conditions near the wall, it has been argued that 
U and a in (6.1) may be replaced by wu and y, so that the equation becomes 


y \u4 
— 0-0225pu2( ~) (6.2) 

\uy 
This may be expected to hold near the wall in any turbulent boundary layer 
flow, outside the viscous sublayer. A discussion of equation (6.2), first 
discovered by Prandtl, is given by Schlichting (1955, p. 404). He shows 





634 M. B. Glauert 


that the formula is adequate at Reynolds numbers up to 10°. This power 
law expression is easier to apply in our present theory than is the logarithmic 
relation which is theoretically more acceptable, and its accuracy should 
be quite adequate for our purposes. At higher Reynolds numbers a formula 
similar to (6.2), with a suitably modified exponent, might be used. 
Equation (6.2) implies that, near the wall, uo y'*. ‘The shearing 
stress t = € Ou/dy has a finite non-zero value at the wall and hence, for 
fixed x, 
ec yl? oc 8, (6.3) 


The assumption we shall now make is that e€ is proportional to wu®, in 
accordance with (6.3). At first we shall assume that this relation holds 
throughout the whole width of the jet, though later we shall apply it only 
in a limited region near the wall. Thus, we choose 
e = m(x)u®, (6.4). 

where m is as yet unspecified. ‘The main reason for choosing € to be a 
function of u rather than of y, is that in consequence y does not appear 
explicitly in the equations, but the behaviour is also satisfactory in that ¢€ 
will not continue to rise beyond the velocity maximum. An additional 
and rather amusing point is that the exterior momentum flux arguments 
of §3 still apply in this case. The last term of equation (3.1) becomes 
mxu® cu/dy, and hence the last term of (3.2) is replaced by }mx?u’, which 
vanishes at the wall and at the outer edge of the jet. Consequently 
equations (3.3) and (3.4) remain true. However, it is not really plausible 
to consider that this form for « is applicable in the outer parts of the jet, 
where the velocity is falling to zero again. 

Consider the conditions for a similarity solution of the boundary layer 
equation (5.1), in which u oc x“, 6 «< x’, Blasius’s formula (6.2) requires 
that 7, oc u748-!4, and hence 


€ oC 43/4934 oc xHa+oy4, (6.5) 
The two sides of (5.1) will depend on x in the same manner if 
a+5b = 4. (6.6) 


Introduce the non-dimensional variables (2.6) and, in accordance with the 
above results, let 





B= P44), n= ye, (6.7) 


and co Ay, (6.8) 
where A and A are constants. ‘The reasons for the introduction of the 


second constant will appear in the next section. Equation (5.1) then gives 


1 reg ” ro 
= (Af'*f") +f" + of? = 0, (6.9) 


where a = (5b- 4)/(5—40). (6.10) 
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The boundary conditions are f’ < 7!" as » >0, f'(y.) = 0, where 7. is 
the value of 7 at the outer edge of the jet. As a result of the peculiar choice 
of ¢€, 7,, is finite, as will be seen later. From (6.5) and (6.10), we have 


9a 445 

on gees a ee 6.11 

5+ 4a ++ 4a ( ) 

The integral arguments of §2 may be applied as before to show that 

“= 2, a result which also follows from the considerations of exterior 
momentum flux given above. With this value, (6.11) shows that 

a= —1-38, b = 1-08, (6.12) 


and equation (6.9) may be integrated directly. By a suitable choice of A, 
we may take A = 1, and the equation then becomes 


£ (ff) ff" +f? = 0. (6.13) 


dy 
Multiplying by f and integrating, we obtain 
PY -MO+PT = 0, (6.14) 
the constant being zero in view of the required behaviour of f near » = 0. 
Multiply (6.14) by 7f-1°’, and integrate again, 
f-*8f'7 + 2 f%8 = constant. (6.15) 


) 


Choose the constant so that f=1 when f’=0. Solutions for other 
values of the constant are then expressible as X7f(X~°n), for some constant 
X. We then have 


of? = fis — f2, (6.16, 


It is convenient to write g(£) = f78, where € = {(@)". Equation (6.15) 
then shows that 


g" = ]— 9? 
P “9g dg a 
and hence é= | Tae (6.17) 


This integral was computed numerically, and at each stage the corre- 
sponding values for f, f’ and were calculated. ‘The value of € at the edge 
of the jet, g = 1, is expressible as a beta function, and thus a useful check 
was made on the accuracy of the work. The variations with 7 of f and f’ 
are illustrated in figure 2. The solution is satisfactory up to the velocity 
maximum, but the rapid fall of velocity in the outer region, a consequence 
of the assumption € wu, is quite unacceptable. A more realistic assumption 
is that € is constant in the outer part of the jet, in accordance with Prandtl’s 
hypothesis. We now consider how this may be achieved. 
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Figure 2. Turbulent wall jet, eddy viscosity proportional to u®. Variation of mass 
flux (f) and velocity (f’) with distance from the wall (7). 


7. MATCHING OF SOLUTIONS 

The first thing to realize is that, if the inner part of the wall jet is to be 
governed by Blasius’s formula, and the outer part by Prandtl’s hypothesis, 
complete similarity is no longer possible. For, as functions of x, « oc R*4 
in the inner layer, as shown in (6.5), and « o R in the outer layer, where 
R = ud/v is the Reynolds number. With increasing x, u decreases rather 
faster than 6 increases, and hence R decreases slowly. Consequently the 
velocity profile must change gradually as the jet develops, the inner layer 
occupying a progressively larger part of the jet width. This is not dependent 
on the precise forms selected for the eddy viscosity in the parts of the jet 
near and far from the wall. Any method which relates the inner part to 
the presence of the wall, and the outer part to free turbulent flow, will give 
the same conclusion. 

However, the ratio of the eddy viscosities in the outer and inner parts 
of the layer varies only as the fourth root of the slowly changing Reynolds 
number, and we shall ignore this relative variation in seeking to derive a 
solution. A process of matching inner and outer solutions will then enable 
us to estimate the Reynolds number to which a particular matched profile 
is appropriate. It is hoped by this means to predict the main features of 
an actual turbulent wall jet. ‘The assumptions in regard to the eddy viscosity, 
on which the whole method depends, are themselves of such dubious validity 
that it is hardly justifiable to try to develop any more precise procedure at 
this stage. 

We now consider the equations for the inner and outer parts of the jet. 
For the inner layer, equations (6.5) to (6.11) continue to apply unchanged. 
In the outer layer we retain the same variables, but, in place of (6.8), we 


write € =Axe%y, (7.1) 
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which is proportional to R*4, instead of being proportional to R as indicated 
by Prandtl’s hypothesis. ‘This is the price that must be paid if we are 
still to have a solution in terms of the similarity variables. ‘The equations 


are aR . ; 7 
— (Afi Sf’) + Afi tof? = 0, (7.2) 


fo tho fo + af." = 0, (7.3) 


where the suffixes denote the inner and outer layers respectively. The 
values of a and 3, the similarity exponents, are given as functions of « by 
(6.11). As boundary conditions, we have fj7~-!?7— constant as 7-0, 
and f;(0o) = 0. Clearly certain further conditions must hold at the junction 
of the inner and outer layers. ‘The value « = 2 was determined in previous 
cases by exterior momentum flux arguments. These are no longer helpful, 
so « remains arbitrary. Indeed it is the variation of « along the jet length 
which represents the gradual change of jet profile, as different values of « 
will be appropriate at different Reynolds numbers. 

Details of the calculations involved in the integration of (7.2) and (7.3), 
for general values of «, are deferred to the next section. For the present, 
we merely assert that such solutions are available. In (7.2) it is convenient 


to write 
fa(n) = A’*A(m), (7.4) 


so that the equation becomes 


dy 


Se = ete 
a (forfe) thofs tafe? = 0. (7.5) 


Now if one solution of (7.5), with the correct behaviour near 7 = 0, has been 
found as f,(y), then, as in §6, the general solution can be written as 
X‘f,(X—n), for some constant X. The corresponding solution of (7.2) is 


ful) = AXP (Xn). (7.6) 

The conditions to be satisfied at the junction between the inner and 

outer layers will now be considered. Suppose that the join is to be at the 

velocity maximum, and here y% and u are to be continuous. ‘Thus at the 
velocity maximum 7 = 7, 


fo=fi=fe=9%  fo=fi=AMX fy fo =fi= AWX'F2. (7.7) 
Since the solutions f, and f, are known, the values of A and X can be found 
at once. ‘The eddy viscosity is discontinuous at 7 =7,,, which is an 
apparent drawback. However, the shearing stress is continuous, being 
zero on each side of the join. 

The eddy viscosity in the outer layer may now be determined as follows. 
Blasius’s formula (6.2) shows that, near the wall, 


T 0 1 se »\1/4 
7 sex = 0-0225u( =) ; (7.8) 
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With the expressions of $6, this takes the form 


AM Af (nfi)"* = 0-0225 (7.9) 
for small 7. If f{ 7+’ > D as 7 + 0, (7.9) gives 
AAD™4 = 0-1575, (7.10) 


and if f37 1‘ > E as 7 0, (7.6) shows that D = A-)°X®*E, and hence 
(7.10) becomes 


A= 0-1575A} 20 27 4-21 4 


5/4 f! 
= 0-1575E na( Ze) “fom (7.11) 
Som Jom 

from (7.6), where the suffix m denotes the value at the velocity maximum. 
All the terms of the right-hand side of (7.11) are supposed known, for any «, 
and thus A is determined. Now, according to Prandtl’s hypothesis, in the 
outer layer, 

€ = KUy, Oy, (7.12) 
where u,, is the maximum velocity, 6, a typical measure of the jet width, 
and « a universal constant; i.e. one which should not depend on z The 
jet Reynolds number may be defined, more precisely than previously, as 
R = u,,8,/v, and it now follows from (6.7) and (7.1) that 





d= Kn fon)®4R', (7.13) 
where 7, is the value of » corresponding to 6,. Hence, from (7.11), 
R14 — 575 F-21/4 fam . * fom ' \-3/4 9 
KR4 = 0-1575E a ii inf) (7.14) 
J Om J2m 


8. INTEGRATIONS AND RESULTS 
The integration of equations (7.3) and (7.5), for general values of x, 
will now be considered. As in §4, we may choose the solution of (7.3) for 
which f,( 00) = 1. A series expansion, valid for large values of », is readily 
obtained in the form 


fun) = 1—e" + A(L + ae ®”— (1 +9(5 + 4a)e-® + 


4+ (1 +2)(34453042122)e47-..., (8.1) 
where the coefficient of e~” has been arbitrarily taken as unity. This is 
equivalent to making a suitable choice of the zero of 7. From the series, 
fy and its derivatives were calculated at a suitable set of values of y, and 
the solution of (7.3) was then extended to smaller values of » by numerical 
integration, until the zero of fj was passed. The values of f,,,, fj, and 7, 
were then estimated, », being chosen as the interval between the point at 
which fj =f%,, and the point at which fj = 3fj,,. This procedure was 
carried out for « = 1-1, 1-2, 1-4, 1-6 and 1-8. The solution for « = 2 was 
found in $4, and, for « = 1, equation (7.3) is that corresponding to a plane 
free jet, with solution f, = tanh}. The values of fy,,, fo, and 7, thus 
determined are shown in figure 3, and lie on smooth curves, considered 
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as functions of « From this figure, the values corresponding to any 
required « in the range can be read off immediately, to sufficient accuracy. 

A point which emerges from the integrations, is the great insensitivity 
to the value of « of the velocity profile in the region beyond the velocity 
maximum. In figure 4, where fj/fj,, is plotted as a function of (y—7,,)/n, 
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Figure 3. Results for a turbulent wall jet. Values of mass flux (fj,) and velocity 
(fim) at the velocity maximum. Also thickness of the inner layer (7,,) and 
outer layer (7). 





Figure 4. Comparison of velocity profiles. ————— « = 1, ---~ a = 2. 


the curves for the extreme values « = 1 and « = 2 are seen to be so close to 
each other that no experiment could be expected to differentiate between 
them. Thus the outer part of the velocity profile of a laminar wall jet is 
virtually the same as that of a plane free jet. 
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For the inner part of the wall jet, we wish to determine the solution of 
(7.5) for which f3 4-!* > constant as 7 > 0. This implies that f, fo/fo? > 4 
as 7 >0. Also, at the velocity maximum, which is the farthest point to 
which the solution need be carried, fy = 0. Thus the second term of (7.5) 
is relatively unimportant compared with the third, throughout the range 
of integration. This suggests that a satisfactory approximation might be 
obtained by replacing the term f, ff in (7.5) by f;°, where @ is a small positive 
constant, not depending on x. The equation then becomes 


t Off) +(a+ 0)fo* = 0. (8.2) 


i; 
On writing f(y) = B f(y), BS = (2+ 6)/(~+9), (8.3) 
(8.2) becomes “4 (f35f3) +(2 + 6)f,? = 0, (8.4) 
an ~ 


in which « does not appear. The reason for the particular choice of the 
substitution (8.3) is that, by reversing the arguments given above, equation 
(8.4) may be replaced by (7.5) with « = 2. This is equation (6.13), of 
which the solution was found in §6. Denoting this known solution by 
f(y), we obtain the required solution of (7.5) from (8.3), once the value 
of @ has been chosen. 

The accuracy of this approach can be judged by a consideration of the 
extreme case x = 1. Equation (7.5) can in this case be integrated twice 


immediately, to give 
ft +48 = Cr, (8.5) 


where C is a constant. As previously, C may be chosen arbitrarily, and 
the solution of (8.5) for a general value of C is then expressible as Y‘f,( Y~°n) 
for some constant Y. With C = 4, equation (8.5) was integrated by a 
method of successive approximations. The following results were obtained. 
Ac seni cr i? " 
As . | 0, fe a hi { (8.6) 
When f3;=0, f,=0-188, f= 0-759, n = 0-269. } 
These exact values may be compared with those we should deduce 
on the basis of = approximation developed above, using the relation 
f3(n) = BY‘f.(Y~°n), and the following properties of f,(7) which were 
unk; in the integrations of § 6. 
As 7 > 0, fg? = 1-322. 
When fs=0, fs=0-479, f3=1:091, 1 =0-478. | 
With a suitable choice of the two constants B and Y, (8.7) should give 
values similar to those of (8.6) for the four quantities concerned. ‘The 


(8.7) 


value of » agrees if Y = 1-121, and the values of B then required to give 
agreement of f,, f; and lim(f3 7-?") are 1-140, 1-142 and 1-140 respectively. 
70 


This is a remarkable justification of the ideas we have used. Equation (8.3) 
now shows that @ = 0-07, a very acceptable value in view of the arguments 
leading to (8.2). 








bl till 
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Thus for the inner layer, as for the outer, the form of the velocity 
distribution is effectively independent of x However, the velocity 
distribution for the whole wall jet does vary with «, due to changes in 
the relative scaling factors in the two parts. The gradient near the wall 
becomes steeper as x decreases. A typical velocity profile, that corre- 
sponding to x = 1-4, is shown in figure 5. The width of the inner layer 





ee 


i 1 1 





Figure 5. Velocity profile for a turbulent wall jet, « 1-4. 
is proportional to »,, = 1-09f,,,/f),, and that of the outer layer to 7, and 
so the velocity profiles corresponding to other values of x may readily be 
deduced from the curves of figure 3, where both 7,, and 7, are given. 
The results for the inner layer, obtained above, may now be inserted 
in (7.14) to give 
0-0275 


«RM = ——— 


x. + (0-07 


> 


(Som) ; “(fom )! —? “’ (5.5) 





and on substituting the values for the outer layer, from figure 3, we obtain 
the results shown in table 1. 








KR : a ) 
kK 0-012) 
1-100 |) O45 1-4 - 10° 1630 | 1-005) | 
12 | O970 | 9 44 < 108 1096 | 1-010 | 
1-3 0-102 | 5200 1-139 1-015 
1-4 Q-O707 1200 1-178 1-019 
1-5 0-0530 380 1-214 1-023 
1-6 00-0417 150) 1-247 1-026 } 
1-8 0-0287 33 1-307 1-033 | 
20 0215 1-359 1-038 | 
Table 1 


The value of the constant « has to be determined empirically. 
Preliminary experiments (Bakke 1956), at a Reynolds number in the 


F.M. 2U 








642 M. B, Glauert 


neighbourhood of 5000, recorded a velocity profile in good agreement 
with that predicted for «= 1-3. This indicates that « = 0-012 is an 
appropriate choice, and the corresponding values of R are also shown in 
table 1. It may be noted that this value of « is approximately one-half 
of that found to be appropriate for a circular jet, and one-third of that for 
a plane free jet. More accurate experiments might indicate a modified 
value to be preferable. ‘The extreme insensitivity of « to the value of the 
Reynolds number powerfully supports the assumption of approximate 
similarity, on which our whole theory has been based. 

Summing up, table 1 represents our main conclusions in regard to the 
behaviour of a turbulent wall jet, and it may be hoped that these would 
not prove unduly sensitive to changes in our assumptions as to the nature 
of the turbulent shearing stress. ‘The predicted velocity distribution is 
given by figure 5, with suitable scaling of the inner and outer layers as 
discussed above. ‘The similarity exponents, specifying the ratio of change 
of the maximum velocity and the jet width, are given as functions of « by 
equation (6.11); a few comments on the accuracy of this equation are 
made in the next and final section. 


9, VALUES OF SIMILARITY EXPONENTS 

‘The one purely arbitrary feature of the matching procedure developed 
in $7, was the decision to link the coordinate system to the inner layer, 
and to ignore variations in the outer layer. An equally valid alternative 
would be to choose coordinates appropriate for the outer layer, and to 
assume the skin-friction to have the value required by Blasius’s formula 
at the particular Reynolds number. ‘The whole analysis would go through 
as before, and identical velocity profiles would be deduced. In place of 
(6.11) we should have, as in §5, 


a=-——, b=, (9.1) 


as is easily verified. ‘hus if x = 1-3, (6.11) gives 
a -1-15, b = 1-03, (9.2) 
while (9.1) gives 
a= —1:13, b= 1. (9.3) 
Although the difference between (9.2) and (9.3) is not great, it would be 
satisfying to see how it arises. 
For x= 1-3, equation (9.2) shows that Raw. If R changes 
from that corresponding to x = 1-3 at x = x4, to that corresponding to 
x = 1-4 at v=x,, then approximately, from table 1, 


7” 5 20(0\ 1/0712 
 .. (; ) ~ 200000. (9.4) 


ae 1200 


Incidentally, this figure strikingly illustrates the slowness of the change in 
velocity profile. Consider a wall jet which has developed from x, to xp. 
On the theory of §7, its profile at x, is still that appropriate to « = 1-3. 
However, in fact the profile should be that corresponding to « = 1-4, with 








> ~~ ~ 
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the scales of f and » suitably adjusted. If the maximum velocity and width 
of the outer layer are u,., and 5,., for x = 1-3, and u,., and 4,.4 for « = 1-4, 
then if the adjustment is made so that the mass flux, and also the velocity 
gradient in the inner layer (for which the solution is known to be accurate), 
are unchanged, a rough calculation shows that, approximately, 


y.¢/%y.¢ = 1:12 = (200000), 
§4.4/S..2 = 0-84 = (200000), 
Adding these contributions from the changing profile to the values of (9.2), 
we have as revised estimates of the similarity exponents 
a= —1-14, 6= 1-015. (9.6) 
‘The fact that these figures are precisely intermediate between (9.2) and (9.3) 
must be considered entirely fortuitous, in view of the roughness of the 
calculations, but the order of magnitude of the differences arising from the 
alternative choice of coordinates is satisfactorily explained. Physically, 
we may argue from (7.1), that, with « constant, the velocity gradient in the 
outer layer becomes progressively too small as x increases, and hence an 
adjustment qualitatively similar to (9.5) is required. For the alternative 
approach, leading to (9.1), the velocity gradient in the inner layer becomes 


(9.5) 


progressively too steep as x increases, and so necessitates an adjustment 
in the opposite direction. 

The discrepancies between (6.11) and (9.1) are probably no larger than 
other errors inherent in the whole method of analysis, but in practice the 
mean of the two might well be adopted as the best prediction of our theory. 
These mean values are also given in table 1. 

For a plane turbulent wall jet, the analysis of the last sections all applies 
unchanged, with appropriately modified values of the similarity exponents. 
Equation (6.11) is replaced by 








7 
54 4a’ 54a’ aa 
and (9.1) is replaced by 





b = | 2 9.5 
1 Ta ; ( §) 
The value of the constant « must be determined afresh, as there is no reason 
to suppose that it still has the same value as that for a radial wall jet. 
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SUMMARY 

The flow rate of liquid metals is commonly measured by 
electromagnetic flowmeters. In these the fluid moves through 
a region of transverse magnetic field, inducing a_ potential 
difference between two electrodes on the walls of the pipe. 
The ratio of signal to flow rate is dependent on the velocity 
profile, and this is affected by electromagnetic forces. 

In this paper the ultimate steady velocity profile and its 
associated pressure gradient and induced potential are calculated 
for the case of laminar flow in a circular pipe whose walls are 
conducting but without contact resistance. Laminar flow is 
encouraged by a transverse field. When the fluid conductivity 
and field strength are sufficiently high, boundary layers occur 
with a thickness inversely proportional to normal field intensity. 
The induced potential difference is then 0-926 of the value 
corresponding to the case of uniform velocity if the walls are 
non-conducting. 

‘The distance the fluid must travel after entering the transverse 
field before the steady state is reached is next estimated by a 
Rayleigh approximation. ‘The inlet velocity is taken to be uniform 
and ettects which occur at the edge of the field are neglected. 
‘The process falls into two stages, first a boundary-layer growth 
and then an adjustment of the velocity away from the walls, 
occupying a much greater length of pipe. ‘The entry length is 
shorter than it is in the case of flow in a rectangular pipe, but is 
still too long for appreciable distortion of the velocity profile to 
occur within practical flowmeters except at low flow rates. The 
pressure drop associated with the adjustment of the velocity 
profile is found to be independent of field strength, if this is 
high, and about one-eighth of the drop which occurs in the 
non-conducting case. 

Experiments are described in which steady-state pressure 
gradients and induced potential differences were measured in 
mercury flowing along Perspex pipes of 0-5 and 0-25 in. bore 
in transverse fields up to 14500 gauss. The results confirmed 


the steady-state theory within the limitations of experimental 
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accuracy and the assumption in the theory of high conductivity 
and an intense field. ‘The experiments also covered the entry 
region in many cases, and showed that the theoretical entry 
lengths were correct in order of magnitude but over-estimated. 
However, the exact entry condition was uncertain, and steady 
readings were difficult to obtain in the entry region. 


1. INTRODUCTION 

Ihe rate of flow of a liquid may be measured by passing it along a pipe 
under a transverse magnetic field. The potential difference induced between 
two electrodes situated at the ends of a diameter of the pipe perpendicular 
to the flow and to the field may then be used to indicate the flow rate. 

The ratio of the induced potential difference to the flow rate is affected 
by several factors, including the conductivity of the pipe walls and the form 
of the velocity profile. Ifthe pipe is circular, the induced potential difference 
for a given flow rate is insensitive to the form of the velocity profile, provided 
the velocity is dependent only on distance from the axis of the pipe (Kolin 
1945). ‘This is usually the case for laminar or turbulent flow in the absence 
of upstream disturbances. If the profile departs from radial symmetry, 
the ratio of potential difference to flow rate can vary very widely (Shercliff 
1954, 1955). 

If the liquid is a good conductor and the magnetic field is sufficiently 
strong, the velocity profile can be appreciably distorted by electromagnetic 
forces (Hartmann 1937, Shercliff 1953), and the induced potential seriously 
affected. ‘This can occur when electromagnetic flowmeters are used with 
mercury or liquid sodium. 

The full distortion of an entry profile by a transverse field takes place 
over an entry length, which may be large or small in comparison with the 
distance the fluid traverses between the onset of the transverse field and 
the electrodes. The transverse field cannot have an abrupt edge, but 
nevertheless may change rapidly from negligible proportions to a high 
uniform intensity. 

This paper examines theoretically the steady state of laminar flow in 
a circular pipe under a transverse field and, by means of a Rayleigh 
approximation, the entry length necessary to achieve it. Experiments are 
described which support the theoretical conclusions. 

There is considerable interest in laminar flow since it persists at high 
Reynolds numbers in the presence of a transverse field (Murgatroyd 1953, 
Lock 1955). The experiments described include cases of laminar and 
turbulent flow. 

In addition to the effect of a uniform transverse field on the velocity 
profile, there is another distorting effect due to secondary currents where 
the fluid enters the region of transverse field. This effect is ignored in the 


present analysis. 
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2. GOVERNING EQUATIONS 
As figure 1 shows, we select Cartesian axes (x, y, =) such that the 3-axis 
lies along the axis of the pipe, and the x-axis is parallel to the imposed 
uniform field H/, existing outside the pipe. 


y 





[ 

eg 

OT \ 
- Ho 
Figure 1. Axes and electrodes 

Various authors (e.g. Bullard 1955) have derived the basic equations 
for the motion of conducting fluids, pointing out that the displacement 
current and the convection of charges are negligible. In electromagnetic 

units Maxwell’s equations therefore become 
curlE = —ycH/ct 
and curlH = 4rj, 2) 


in which E, H and j are the electric tield, magnetic field and current density 


vectors respectively, and , is the permeability of the fluid. Ohm’s law 
requires that 
j = o(E + LV x H), (3) 
in which v is the fluid velocity and o is the fluid conductivity. ‘The 
dynamical equation is 
uj x H+7V*v— grad p = picv ct +(v. grad)v}, (4) 


in which p, p and 7 are the fluid pressure, density and viscosity respectively. 
We assume that p is constant, and , is unity in both the fluid and the pipe 
walls. Consequently 
divv=0, and divH = 0. (5) 
For the motion of fluid in a straight pipe under a uniform transverse 
field, where conditions may vary in time but not in the 3-direction (apart 
from the pressure gradient), the equations become linear in the unknown 
quantities. Differentiation of (4) with respect to z shows that grad(¢p/¢z) 
vanishes, and hence that cp/cz is a function of time only. It is compatible 
with the above equations, and with the boundary conditions, for v,,v 


T 


y 
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and H, to vanish and H, to have the constant value Hj. Equations (1) 
to (5) then reduce to two equations in the unknown components /Z. and v_: 

V?H,+47pu0H, ov,/0x = 40 cH, /ct (6) 
and nV2v, + (wH,/47) dH./ox = dp/dz t+ p ev,/et. (7) 
The Laplacian operators comprise two derivatives only. 

The estimation of the entry length for velocity-profile distortion to 
occur after entry into a transverse field is a non-linear problem, as it stands, 
There is, however, an obvious analogy with the above linear case, in which 
there is variation with time instead of with z. We should therefore expect 
results of the right order if we calculate settling times in the linear problem, 
and convert these to entry lengths by multiplying by vp, the mean velocity 
of flow in the pipe, which is necessarily constant. ‘The variables z and ¢ 
will be measured from the onset of the transverse field. 

The equations (6) and (7) can be expressed non-dimensionally by 
means of the substitutions 

U = U,/Vs h = H,/47v, (on), P = —(a?/nv,) cp/ez, 
A = 2/4, Y= vid, T = ynt/pa°*, 
M = pH, a(a'n)'*, and B = 4rpyo7n/p, 


where a is the internal radius of the pipe. We then have 
Ah+ Mev/0X = Bh (3) 
and Av + Moh/0X+ P= 2, (9) 


in which A = d2/0X2+02/AY?, and a dot denotes differentiation with 


ll 


respect to 7. 

The dimensionless fluid property 8 is of the order of 10-® for liquid 
metals. We shall therefore neglect the 8-term in (8). ‘This may be shown 
to be justified except at unusually large values of M, which is of the order 
of 10? in practical flowmeters. The role of / is now limited to that of a 
current stream-function, since 

jy, = —Up(on)'? dh/dx, and 7, = v(on)!20h/y. (10) 

In view of the interest in flowmeters, it is necessary to consider the 
electric potential V induced by the motion. This potential exists, because 
neglecting 4 implies that curl E is negligible. Then E = —gradV, and 
from (3) it follows that ; 

dV /oy = nH, v,—-},/¢. (11) 
It is convenient to express the flowmeter output Vyy, measured between 
the electrodes X and Y (figure 1), in terms of a sensitivity S, where 
Ss = V yy/2H, arp. (12) 
This takes the value unity for uniform or radially-symmetric velocity 
profiles in a non-conducting circular pipe (Kolin 1945). 
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Spatial boundary conditions 

One obvious condition is that v7, and v shall vanish at the wall. ‘The 
other condition concerns the magnetic field and potential. In the absence 
of any contact resistance between the fluid and the walls, the potential is 
everywhere continuous. Also H, and h must be continuous at the fluid-wall 
interface and zero at the exterior of the walls. ‘This secures continuity of 
current flow between the fluid and the walls, and precludes the flow of 
current outside the walls. Any voltage measuring device would take a 
negligibly small current. 

From (10) it follows that j,, the total wall current per unit length of pipe, 
is equal to (a)! *z, A, in which / is taken at the fluid boundary. If as and an 
indicate distances along the boundary and inward normal respectively, as 


shown in figure 2, and we assume that the wall has conductivity o, and 





thickness w (<a), then —jy/wog = (CV /cs)/a, V being in the wall or in 





the fluid at the wall. At the wall, v vanishes and —ograd V = curl H/4z, 
and hence cV /és (yn o)'*z,ch cn. Eliminating cV/cs and j,, we see 
that at the fluid boundary 
wo, oh 
r= = % (13) 
ac Cn 


the other spatial boundary condition, if contact resistance is negligible. We 
shall write c for wo,/ac, a dimensionless quantity. 

In all that follows, complete symmetry in the x and y axes will prevail, 
with the result that v is even in X and Y, and that A is odd in Y and even 
in Y’. 


3. "THE STEADY STATE SOLUTION 

In an earlier paper (Sherclitf 1953), (6) and (7) were solved for the case 
of steady flow in a rectangular channel with non-conducting walls. It 
was found that, when is large, a core of uniform velocity occurs together 
with boundary layers at the walls. We should therefore expect similar 
behaviour in a circular pipe with thin conducting walls, the case which is 
now analysed. 

We make the following assumptions. 


(i) M is large in comparison with unity. 
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(11) ‘The boundary layer thickness is much smaller than the pipe radius, 
so that all variables may be treated locally as functions only of the 
distance an from the wall.’ 

(iii) P is of the same order as Mz, as in the rectangular case, when 
M is large. P is here a constant, and z is of order unity. 

(iv) In the core, Az is negligible in (9), signifying that viscous forces 


are negligible. 


The boundary layer region 
The application of assumption (ii) to (8) and (9) for a point on a radius 
vector inclined at an angle @ to the x-axis (see figure 2) yields the equations 


h" —mzv' =v" —mh'+ P=0, (44) 


in which m = Mcos#, and dashes denote ditferentiations with respect to n. 
The boundary conditions are that 7 = 0 and h=ch’ when n= 0, and 
that v and A must become relatively steady outside the layer, approaching 
the local core values v7, and h,. The solution 1s 


v = v,{1—exp(—mn)}, (15) 


h = h(O)—v+ Pnjm. (16) 


The boundary-layer thickness is of the order of 1/m, which is proportional 
to the reciprocal of the normal field intensity. ‘The last term in (16) is 
negligible except where 6 approaches jz, by virtue of assumption (iii) 
and because ” is smaller than the boundary-layer thickness. ‘Thus the 
pressure forces in the boundary layers are negligible, except where 4 
approaches }z. Equation (16) now gives h = h(0) —v,{1—exp(—mn)}, and 
the condition (13) requires that A(0) cmz, v.cMcos@. It follows 
that the total currents towards .Y per unit length of pipe are (a7)! *z9(h, — Ay), 
which is equal to — (o7)! #7) v,, in the boundary layer, and — (oy)! *v,7v,.cM cos 6 
in the wall. ‘The current in the layer is proportional to the core velocity, 
because the wail shear stress, proportional to (core velocity)/(layer thickness) 
and hence to (core velocity) x (normal field intensity), is balanced by the 
magnetic force on the layer, proportional to (current in layer) x (normal 


field intensity). 


The core region 

From assumption (iv) we have ch,/°-X+P/M =0, and_ hence 
h.= — PX/M, since h, is odd in X. This indicates that the current 
density in the core is uniform and parallel to the y-axis, and equal to 
Pv,(on)'?/Ma. The electromagneti™force just balances the pressure 
gradient. Also, Ah, vanishes, and (8) shows that, in the core, v, = v,(Y), 
a function of Y only. ‘The current in the core must all return in the two 
boundary layers and walls, as shown in figure 3. Across any section AB, 
remote from X or Y, the total core current per unit length of pipe is 
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2X, Pu(on)!?)M, in which X, = (1— Y?)!?, if we neglect the boundary-laver 
thickness in comparison with the length 2X, of AB. Continuity of current 
flow, or of h, demands that PX,)M = v,(1+cMcos8@), and hence that 

Le x, 
Wil+cMX,)’ aia 
the required function of }. Equation (17) ceases to be valid as Y approaches 


}. 





The core velocity is not determined by local dynamical conditions, but 
indirectly by means of the current distribution. We may note in passing 
that when the walls are conducting, there is a downstream electromagnetic 
force on them with a corresponding net upstream force on the fluid. The 


pressure gradient has to counteract this as well as the viscous shear forces 
it the wall. When the walls do not conduct, there is no net electromagnetic 
force on the fluid. In neither case does the motion produce a force on the 


magnet supplying the transverse field, if edge etfects are excluded. 

It is now appropriate to examine some of the earlier assumptions. 
Equation (17) shows that Az, is of the order of P/V, and is negligible in (9), 
except when Y approaches +1. “Vhe boundary-layer thickness is small 
except where # approaches $7. Obviously, the approximate treatment of 
the core and boundary layers fails near X and Y, but we assume that the 
solution is approximately valid elsewhere when is large and the obscure 
regions are small. 

Equation (17) is compatible with assumption (iti) provided cM is of 
the order of unity or less, and c is of the order of 1/M or less. For sodium 
flowing in a stainless steel pipe with reasonably thin walls, this condition 
is satisfied. The quantity c is independent of the inside diameter of the 
pipe, and measures the relative conductivities of wall and boundary layer. 
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When c.W is larger, the pressure gradient becomes important in the boundary 
layers, and the wall and core currents increase until the normal currents 
in the boundary layers become comparable with the tangential currents. 
Taking / as a function of m alone implies neglecting the normal currents. 

Since the boundary layers are small, we may integrate (17) over the 
cross-section and relate P to the mean velocity vw». ‘This gives the relation 

at p eave »lon\l2 
- — a (1 +0-883cM —0-019(cM)?+...), (18) 
which is nearly linear in cM, 

The potential distribution may be found to the same accuracy from (11), 
in which nHyv.(j,/c) is of the order of WW. We therefore neglect j,/o and 
integrate along VY, assuming that no significant contribution to Vy, 
comes from the small obscure regions where v, and j, must be small. ‘The 
variation of S with cM is shown in figure +. When the walls are non- 
conducting, c vanishes and S = 377/32 = 0-926. Increasing the wall 





conductivity increases S, owing to the change in the velocity profile. Any 
decreasing tendency due to the short-circuit paths previded by the walls 
was neglected when j,/¢ was omitted from (11). 
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Figure +. The effect of wall conductivity on sensitivity. 


For low values of .17, (8) and (9) may easily be solved in polar coordinates 
to vield solutions in series of ascending powers of VM. In the steady state, 
the sensitivity is given by a series which begins with the terms 


S _ =a sake (19) 


” {ee Sais 


4. 'THE ENTRY PROBLEM 


The solution of (8) (less the 8-term) and (9) is now considered for the 
case of a circular pipe with non-conducting walls, where the initial velocity 
is uniform and equal to vp. It is not possible to impose an initial condition 
on h, since it is determined by (8). The spatial boundary conditions are 
that 7 and h/ shall vanish at the walls. 
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In order that the linear time-dependent case may be analogous to the 
original problem of entry into a region of transverse field, we must also 
impose the condition that v, and the integral of v over the pipe cross-section 
are constant. 

It has been shown (Shercliff 1956) that when is large and the fluid 
is flowing between parallel planes perpendicular to the field, a uniform 
entry-profile adjusts to its steady state in a dimensionless time T of the 
order of 1/M?. ‘The dimensional time is of the order of p/u?Hio, and the 
distance between the planes is irrelevant since the boundary-layer thicknesses 
are small in comparison. ‘The settling process in this case simply involves 
the growth of boundary layers, the core velocity being virtually unaffected. 

We should expect settling times for a circular pipe to be of the same order 
as for a square one, although somewhat less since the regions where the 
boundary layers are thick are smaller. For square pipes the dimensionless 
settling time is of the order of 1/.W (Sherclitf 1956). 

The change from a uniform initial velocity in a circular pipe to the 
final steady profile, as already discussed, would be expected to fall into 
two stages as follows. 

(a) During the first stage, lasting for a dimensionless time of the order 
of 1/M?*, thin boundary layers are created as in the parallel-plane case, the 
core velocity being barely affected. Initially, no currents flow until viscous 
forces alone have produced a very small boundary layer. Subsequently, 
electromagnetic forces restrict the development of the layer to a thickness 
proportional to the reciprocal of the normal field intensity. 

(6) In the second stage, occupying a much longer time of the order of 
1/M, the core profile is changed by electromagnetic forces to its final form, 
the boundary layer thickness being unchanged. 

The motion in the regions near X and Y where the simple boundary 
layer theory fails is again left obscure on the assumption that the results 
are not seriously attected. 

For both stages we shall make the following assumptions suggested by 
the steady state solution. 

(i) MW is large in comparison with unity. 

(ii) ‘he boundary layer thickness is always much smaller than a, so 
that all variables may be treated locally as functions only of 7 and n. 

(iii) A and P/M are, at most, of the order of unity throughout. P may 
be related to the shear stresses at the walls by applying Green’s theorem 
to the integral of (9) over the cross-section or the part in the first quadrant. 
This gives 


i" Ov Pr 


| dé = ; (20) 


Jo On 


(iv) In the core, Av, is negligible in (9), and Ah, is of the order of unity. 
‘The solution proves to be compatible with these assumptions except near X 
and Y. From (8) (less the §-term), it follows that év,/¢X is negligible and 


v, is a function of Y and T only. We may note that this is true initially 
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and finally. Equation (9) degenerates to 

M ch,/0X + P = a, (21) 
and we see that ch,/0X is also independent of X. Since h is odd in X, 
ch,|aX =h,|X = ho/Xo, the value hy referring to the edge of the thin 
boundary layer. Equation (21) becomes 

Mh,|X)+ P = v- 22) 
On the basis of assumption (iii), this equation implies that v, is of the order 
ot M at most. 


The first stage 

This is expected to last a time of the order of 1/M?, and so the change 
in v, will be of the order of 1/M and thus negligible. We shall therefore 
take v, as unity during this stage. 

Assumption (ii) applied to (8) and (9) for a point in the boundary layer 
on a radius vector inclined at an angle 6 to the x-axis yields the equations 

h" —mv' = v" —mh'-v+P=0. (23) 

Both h and wv vanish when n = 0, and, when a is large, must approach the 
local core values fy and 1 respectively. From (23), h’—mv = f(T), and 
hence v” — mv —wv = mf(T)— P (independent of n) = — m?, since when n is 
large, v = 1 and vw” and v vanish. In addition, v= 1 when 7 =0. The 
solution is 
exp(—m*T) 
2(nm?71)'2 ~ 


ore nf (mn—u\ — f — (mn+u/? 
| exp( i)| exp TT) exp( a) du, 


v = 1—exp(—mn)+ 


in which the integral is obviously less than unity. 

The form of the solution when 7 is a low multiple of 1/m? is given 
closely by the first two terms. The boundary layer then resembles the 
steady-state one in extent, but the velocity rises to v) across every section 
of the layer. The first stage can be considered complete in a time of the 
order of 1/M?, even though 1/m? is much larger in the obscure regions 
near X and Y. 

Since P, which is of the order of M, is negligible in comparison with m* 
except near XY and Y, then f(7) = —m, and so hf’ = m(v—1). This 
equation and the condition A = 0 at the wall determine the behaviour 
of h, hy and hence h, during the first stage. It can be shown that /, decreases 
steadily from zero, reaching —1 at the end of the first stage. 

Figure 5 shows the current distribution at the end of the first stage. 
‘The boundary-layer current per unit length of the pipe is (oy)! ?z, at all 
sections of the boundary layer. ‘There is singular behaviour near X and Y, 
where the currents issue into the core. 

During the first stage, 0v/én is of the order of m or M, except initially, 
when it is of the order of 7-1? and tends to infinity. Equation (20) shows 
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that P is, as assumed, of the order of M except when T is very small and 
M is an irrelevant variable. 








Figure 5. Current flow at the end of the first stage of the entrv process. 
ae 


The second stage 

We shall see that v7 in the boundary layers is now of the same order 
asv.. Hence vis, at most, of the order of .W, and, like P, is negligible in (23). 
The boundary-layer analysis then gives v = —h =7,{1—exp(—mn)}, in 
which @, varies with time, the extent of the boundary layers remaining 
constant. In particular hy -v, and cv/en = mv, at the wall. From (20) 
it follows that 


Cc 


“hn “1 
P= : | M cos6v,d@ = ae | vay. (24) 
“J0 = we 
Again P is seen to be of the order of 7. From (22) we mav now deduce 
the linear equation 
it Ne! ae 9) $ 5 
Avi x. | Rate ’ (2: } 
which governs the second stage of the motion. It is obvious that its duration 
will be of the order of 1/4. The initial condition is v, = 1, and no spatial 
boundary condition is necessary. 


v v, 4/3 


In calculating S we use (11) applied in the core, the term j,/o again 
being neglected. If no significant contribution to Vy; comes from the 
small obscure regions, it follows that 


S=]| %, ay - (26) 


and we observe that P and S are simply proportional during the second 
stage. Equation (25) may be rewritten 
a v + 


fe 4% Sigur 27 
g°*s.*3° vig 
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which has the integrating factor exp(/7).X,) and the solution 
r -MT 
v, =exp(— MT/X,)< 1+ : | S(ujexp(u/Xy) du > 


4 ;MT 
=exp(—MT/X,)+— | S(ujexp|(u-MT)/X,' du. 


Tso 


Application of the condition (26) yields the integral equation 
-M1 
S(MT) = H(MT)+ : | S(u)H(MT—u) du, (28) 


“eog 


gi 
where H(p) = | exp(—p/X,) dY. 


J ¢ 

A numerical solution of (28) has been calculated, taking steps of 0-1 in MT 
and wu. Initially S is unity, the velocity profile being almost wholly radially- 
symmetric (Kolin 1945); and finally the steady value 377/32 is reached. 
The continuous curve in figure 6 shows the variation of Sor 
— {7a/4nH, v(on)'*} 0p/ez as a function of MT, or of Mz/Ra if we 
return to the non-linear entry problem by means of the Rayleigh 
approximation, setting z = vyt and taking R to be the Reynolds number 
pt, a/n based on pipe radius. 
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Figure 6. The entry process. S or dimensionless pressure gradient as a function 
of MJT or Mz/Ra. The experimental points are from measurements of S at 
different values of M, R and z. 


When T is small, the exact variation of S is not physically significant, 
because it is chiefly affected by the behaviour of «, in the obscure regions 
near X and Y, and also because the first stage of the motion 1s still occurring. 
Moreover, the transverse field is not entered instantaneously in practice. 
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From figure 6 it will be observed that S is within 1°, of its final steady 
value when z = Ra/M. This makes an interesting comparison with the 
result for a rectangular pipe (Shercliff 1956) that departures from the 
ultimate steady velocity profile decrease by a factor 1/e in an entry length 
of the order of Ra/M. The departure of S from its ultimate value would 
be expected to decay similarly; and for comparison a chain-dotted 
exponential, decreasing by 1/e in a length Ra/M but having the terminal 
values appropriate to the circular case, has been added to figure 6. ‘This 
emphasizes the fact that when M is large, entry lengths in rectangular 
pipes are considerably longer than in circular pipes. 

‘The analysis has not been extended to cases where the entry velocity 
is not uniform. ‘There is no reason for expecting the entry length to be 
very different. ‘The analysis could easily be adapted for any initial variation 
of velocity in the y-direction. If the core velocity varies in the x-direction, 
such variations would probably be eliminated in the first stage by eddy 
currents in the xy-plane. ‘This is suggested by the fact that for flow 
between parallel planes (Shercliff 1956), the entry length is of the order 
of Ra/M? when the core velocity varies in the x-direction. 

Solutions in the form of series in ascending powers of M, valid when V/ 
is small, could be found, but the interest in such cases is slight. 


Pressure loss at entry 

In addition to the pressure loss due to the steady pressure gradient, 
there is a pressure loss 5p associated with the profile adjustment from a 
uniform inlet velocity. No allowance will be made for the further large 
pressure loss due to eddy currents at the edge of the transverse field, as 
discussed by Hartmann (1937). 

The drop 5p may be found by integrating the excess of cp/cz above 
its steady state value from zero to large values of 3. No significant 
contribution comes from the very brief first stage of the motion, despite 
the singularity of P when 7'= 0. We can therefore use the second-stage 


result 
P=4MS/z. (29) 
From (27) we obtain 
O(X 9 Ve) 4+ _ 
oer +O Te, 
MT) 7 
which, when integrated from Y = 0 to Y = 1, shows that 
n 1 
— | XY v.dY¥Y =0 
Tir 


because of (26). Hence 


X,v.aY = const. = if 


initially, which expresses the equation of continuity since the boundary 
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layers are very thin. Also from (27), we get 


0(X; v,) : os 
aa OT Xo UV, = Xo 8. 
oMT) 7 
Integration of this equation from Y = 0 to Y = 1 gives 
( l ° e 8S 7 
Xez dY=— — 


and hence 





xr? Z Y [* 3x” 
oe (s- 7) AMT), 


FR - 
128 37) 6 

since v, = 37X,/8 finally, and the required integral has been determined. 
It follows that 


2 = 0-082, 





dp/ipv2 = 


which contrasts with the larger value 0-667 calculated by the Rayleigh 
approximation (Shercliff 1956) for flow into a circular pipe when magnetic 
effects are negligible. When & is large we see that dp is independent of M. 
A similar result would be expected for flow into rectangular pipes, in 
contrast with the result (Shercliff 1956) that 5p/$p75 = 1/M for flow between 
parallel planes when is large. In this case there is no second stage of the 
motion, and Sp is correspondingly smaller. 


5. EXPERIMENTAL WORK 

Experiments have been performed with the main object of verifying 
the theoretical values of S and cp/cz in the steady flow of mercury in a 
non-conducting circular pipe under a transverse magnetic field at large 
values of M. The available values of 17 were limited by the saturation 
of the electromagnet and by the use of a tube small enough to produce 
measurable pressure gradients. Previous experiments (Hartmann and 
Lazarus 1937) reached a value of M of 18 only, and did not consider entry 
effects. 

The apparatus also yielded approximate results for the variation of S 
and cp/dz during the entry process, since only at the lower flow rates and 
greatest field intensities could steady conditions be reached within the 
extent of field available. At the lower flow rates and field intensities, the 
measured potential and pressure differences were liable to stray thermal 
and other effects, but nevertheless significant evidence in support of the 
theoretical work was obtained. 


Experimental detai. 

Figure 7 illustrates the mercury flow circuit schematically. ‘The test- 
section indicated was either of two alternative Perspex pipes, installed in 
the gap of an electromagnet. One of 0-5 in. internal diameter with seven 
pairs of flowmeter-type clectrodes at 1-5 in. pitch was used for measure- 
ments of S, and one of 0-25 in. diameter with four flush tappings at 3 in. 


F.M, 2X 
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pitch for measurements of dp/dz. ‘The electrodes were in fact small flush 
tappings also, this type having proved very successful. ‘Tappings of x in. 
diameter were the smallest that could be satisfactorily cleared of air. In each 
case the first tapping was 1-5 in. downstream of the edge of the magnet 
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Figure 7. The fluid circuit (not to scale). 


Steady flow occurred through the test-section under a gravity head of 
mercury maintained by an electromagnetic flat linear induction pump in 
the lower limb of the circuit. ‘The cooler removed heat supplied by the 
pump. ‘The mercury entered the test-sections through a streamlined 
contraction from a settling chamber with wire mesh grilles and preceded 
by a honeycomb to eliminate swirl. ‘The velocity would presumably be 
nearly uniform after the contraction. Horizontal piping was avoided to 
aid the removal of air. 

The flow rate was measured by an electromagnetic flowmeter in the 
upper limit of the circuit. Disturbance of this meter by the electromagnet 
was eliminated by the mild-steel screens shown and, when necessary, by 
taking the means of readings with the electromagnet energized alternatively 


in either direction. 

The flowmeter was calibrated before and after the experiments with 
the screens in position using the apparatus illustrated in figure 8. The 
mercury flowed through the meter and the control cock under a constant 
head of 5-5 ft. Runs began when plug G was opened with plug H shut. 
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The mercury rose in the weir tank until a clock was started when contact 
was made with the electrode Z. This was surrounded by a sheath into 
which the mercury entered through a constriction at the bottom. The 
plug H was then opened, admitting mercury to the weigh tank until H was 
shut shortly before the header tank emptied. The clock stopped when 
contact was made at Z for the second time, there being the same amount 
of mercury in the weir tank at both timing instants. Thus the timed interval 
corresponded to the amount admitted to the weigh tank. After weighing, 
the mercury was blown back to the header tank through another pipe. 
Mercury has the advantage that electrical level indicators are easily arranged. 
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Figure 8. The flowmeter calibration rig (not to scale). 


‘The apparatus worked well, even for runs as short as 18sec. The 
flowmeter was found to give a relation between flow rate and output voltage 
which was linear over a wide range and reproducible within +0-5°, of the 
mean in all positions relative to the upstream bend. ‘This reliable 
calibration is to be attributed chiefly to the narrow fluid passage of the 
flowmeter, 1-10 in. deep and 0-04 in. in the field direction. 

The electromagnet had pole faces 5in. by 12 in. and a 1-25 in. gap. 
The water-cooled, low-voltage winding enabled fields of 14800 oersted 

2X2 
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to be attained. There was a sensibly single-valued relation between field 
intensity at the pole-face centre and exciting current for both directions 
of magnetization. Accordingly, the exciting current, measured in a 
temperature-insensitive shunt, was used to indicate field intensity after an 
initial calibration with a Cambridge fluxmeter and search coil had been 
performed. ‘The experimental values of H, are considered accurate to 
0-5°%,, except perhaps at low field intensities. 

The field intensity fell off near the edges of the poles, the fall reaching 
2°,, at the end electrode positions when the field was at its greatest value. 
A simple proportional correction was applied to voltages measured at the 
end electrodes, but not to the pressure differences. The effect of secondary 
currents at the field edges is difficult to allow for, but is small. 

The test-section diameters were measured to 0-5°,, either by travelling 
microscope or weighing with a known length full of mercury. Mercury 
properties were taken from the Liquid Metals Handbook (Lyon 1952) and 
corrected for temperature changes. 

All signals from the various electrodes were measured to 1 or 2 microvolts 
with a potentiometer, the readings at the higher flow rates or lower fields 
being more uncertain owing to increased unsteadiness of the flow. Care 
was necessary to ensure that the leads to the test-section electrodes did not 
link any of the main field as otherwise its small variations would produce 
large spurious signals. ‘This is a point in favour of permanent magnets 
for flowmeters. Slight thermo-electric effects were observed when the 
magnet was running warm. All runs involving small signals were therefore 
made before the magnet had become warm. 

Pressure gradients were measured with the sensitive manometer 
illustrated in figure 9. This instrument could record pressure differences 
to 0-001 in. of methylated spirits, a better fluid than water which is liable 
to wet non-uniformly. Owing to the presence of a warm magnet, however, 
it was impossible to eliminate stray errors of the order of 0-002 in., despite 
thorough lagging. 

The four pipes from the flush tappings in the test-section led to the 
eight pinchcocks, which enabled the pressure difference between any two 
tappings to be measured, and also the ‘zero’ reading of the manometer 
to be established at any time. ‘The manometer itself consisted of an 
inverted U-tube containing air above spirits. ‘The spirits/mercury menisci 
were kept at fixed levels in reservoirs of large area by taking all readings 
with the air/spirits menisci brought to fixed positions in the two sloping 
parts of the U-tube. ‘This also eliminated surface tension errors. It was 
achieved by adjusting the volume of air above the spirits and by altering 
the height of the right-hand sloping tube. Changes in the volume of the 
flexible tube below this had negligible effect, and thermal and other effects 
were reduced by the near-symmetry of the manometer system. A micro- 


meter screw recorded the change in height of the sloping tube. The ‘zero’ 
reading was taken repeatedly as it was found to drift during long runs, 
The manometer took about 15 min to reach a steady state. 
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As the micrometer travel was limited, to measure the larger pressure 
differences the manometer was converted to a spirits/mercury one by 
allowing the spirits to rise, expelling the air, until the mercury menisci 
were in the sloping tubes. Readings taken by both methods were consistent. 
The smallest heads that were measured were about 0:05 in., read to the 
nearest 0-001 in. Care was taken to keep the mercury in the circuit at 
ambient temperature to minimize hydrostatic errors in the inclined test 
section. 
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Figure 9. The manometer system, showing the air/methylated spirits arrangement 
(not to scale). 


Observations were taken at four electromagnet field intensities and a 
range of flow rates, each 50°,, higher than the previous one. ‘The upper 
limit was set by the available head, and the lower by the smallness of the 
signals to be measured by the potentiometer in the presence of stray 
potentials. The lowest signals were measurable only to 2°), accuracy. 


Experimental results: pressure gradients in the 0-25 in. tube 


Figure 10 shows the readings reduced to dimensionless form, the three 
measured mean gradients between tappings being expressed as a multiple 
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of the theoretical gradient 37H, v)(on)!?/8a, valid for steady conditions 
when M is very large. Flow rate is represented by R, the Reynolds number 
based on radius. ‘There are three zones in which the constant-/M curves 
are horizontal, rising slowly and rising precipitately. 

Where the experimental points are all marked as crosses, the three 
pressure differences showed a scatter but no clear falling tendency along 
the pipe. ‘The central pressure difference tended to be slightly high. 
Horizontal lines fit the points reasonably well, and we deduce that this 
zone corresponds to steady flow that is laminar, because the gradient 1s 


independent of R and density. 
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Figure 10. Experimental values of pressure gradient plotted non-dimensionally 
against flow rate at four field intensities. 


At higher flow rates the pressure gradient was observed to fall along the 
pipe as steady conditions were approached, and here the three successive 
points are marked as arrow tail, cross (air,spirits) or circle (spirits/mercury) 
and arrow head. ‘The points suggest that the flow was approaching a 
steady state with the same value of (cp/0z)/(37uH  vo(on)!? 8a) as for the 
lower flow rates. Chain-dotted curves connect the arrow heads approxi- 
mately. When M was 29, steady conditions were barely obtained even 
at the lowest flow rate. ‘The pairs of arrows linked by brackets represent 
the fastest air/spirits readings and the slowest spirits mercury readings, 
consistent within experimental error. During the runs where a steady 
state was still being approached, the readings both of pressure gradient 
and voltage were particularly prone to wander, even though the total flow 


rate was steady. 
Finally the pressure gradients are seen to rise rapidly at a value of R 
that increases with WM. Most of the experimental points are off the figure. 
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With M = 29 or 52 the gradients increased along the tube. There seems 
to be no doubt that a transition to turbulence was being observed. The 
rig did not permit an adequate study of the turbulent regime, and this was 
10t attempted. It appears that provided R<250M approximately, 
laminar flow occurred under the conditions of the experiment, and that, 
in particular, there was laminar flow in all the runs with the 0-5 in. tube, 
described in the next section, except perhaps the fastest runs at the lowest 
field strength. 





' 
- 
= = 
BH ye (en) * \ 
(FoR Sovane) 


OO Box 4 (a en es 
-Bast \ 
Srp y(n)" | 


\ CIRCULAR (FROM HARTMANN) 
(FoR CiRcULAR) 


+ * SOUARE 
: } (THEORETICAL) 


ey i 
CIRCULSR 
( 4 Tus) | 
~ 
+ VALUES FROM FIGURE IO 
+ + } > 
+ 


2 20 “0 60 &© 100 M 120 











0 


Figure 11. Ultimate pressure gradient plotted non-dimensionally against M for 
circular and square tubes. 


The four steady values of (ep/dz)/(374H,yv(on)'?/8a) taken from 
figure 10 are plotted against M in figure 11. It appears likely that the 
curve does approach the value unity asymptotically as M_ increases. 
Experiments at higher values of MW would almost inevitably involve the 
use of liquid sodium. Figure 11 also includes an approximate chain-dotted 
curve taken from the experiments of Hartmann & Lazarus (1937) for values 
of M up to 18. This is clearly compatible with the new results. 


Experimental results: potential differences across the 0-5 in. tube 

Values of S calculated from readings at the seven electrode pairs are 
plotted in figures 12(a) and 12(b) for values of M of 239 and 58. Similar 
intermediate results were obtained with M equal to 156 and 102. Curves 
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have been added to illustrate the trend of the results as R varies. ‘The 
caption ‘entry’ refers to the leading edge of the pole faces. 

Attention is directed first to figure 12(a), since the readings are most 
reliable at the highest field strength. It is apparent that at the lower flow 
rates steady conditions were reached within the test-section, and that 
S then assumed the approximate value 0-92, confirming theoretical 
expectations reasonably well. At higher flow rates, decreasing values 
of S closer to unity occurred, as expected, because the velocity profile was 
not fully distorted from its presumed initial axial symmetry in the length 
available. 
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Figure 12. The fall in sensitivity along the 0-5 in. pipe. 
(a) M 239, (6b) M 58. 

The accuracy fell as the field was reduced owing to the smallness of the 
quantities measured. Nevertheless figure 12 (4) 1s consistent with theoretical 
expectations. Another difficulty was the unsteadiness of the flow in runs 
where the settling process was protracted. It was impossible to achieve 
precise reproducible trends in these cases. 
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The entry process 

‘The scatter of the results makes a close experimental check of the entry 
theory dificult. It is also complicated by the uncertainty of the entry 
conditions. As the streamlined contraction terminated about 1-25 in. before 
the edge of the pole-faces, there would presumably be some distortion of the 
velocity profile caused by the fringe flux before the main field was entered, 
and moreover the secondary currents at the field edges tend to distort the 
velocity profile. Nevertheless, in order to check the entry theory (itself 
approximate) we shall assume conventionally that the velocity was uniform 
at the section opposite the edge of the pole faces, and measure z from this 
section. 

Experimental points for values of M of 239 and 58 have therefore been 
plotted in figure 6, = being measured from the edge of the pole faces. 
Dotted mean curves have been added. ‘The theoretical entry length is 
seen to be correct in order of magnitude, if considerably overestimated. 
The experimental points also tend to confirm that the sensitivity S started 
to fall from unity while z was still negative. 

A closer investigation of the entry process would require an elaborate 
technique to secure known, standard entry conditions. ‘The dependence 
of S on the velocity profile provides a useful technique for studying such 
phenomena as the pipe-entry process, unless the velocity profile is radially 
symmetric. 

Equation (28) indicates that 0p/oz and S should be mutually proportional 
at high values of WM (except perhaps when the entry length is so long that 
the first stage with its higher value of cpj¢z becomes important). ‘There 
is rough confirmation of this result in that ¢p/¢z and S both fall by amounts 
as large as 8°%, during the settling process. Figure 6 shows that S is usually 
within 1°,, of the steady state value when Mz/Ra=1. A dotted curve 
for which R/M = 72, the value of z/a which corresponds to a point 
midway between the last two pressure tappings in the 0-25 in. tube, has 
therefore been added to figure 10. It is found to connect points where 
op/ez is roughly 1°, above the steady state value. ‘Thus both sets of 
readings give entry lengths of the same order of magnitude. 

As M increases, the steady state value of S approaches its asymptotic 
value more rapidly than (0p/éz)/(37pH, vp(on)'*/8a) tends to unity, simply 
because S (unlike ep/ez) only varies through a small range over the whole 
range of M. It was clearly impossible to find the variation of S with M 
from the present experiments. Hartmann & Lazarus did not measure 
potentials. 

For comparison, a theoretical curve (Shercliff 1953) showing the 
variation of the steady (cp/ez)/(uHyvo(on)'?/a) against M for a square 
pipe of side 2a, with the field parallel to a side, has been added to figure 11. 
‘This curve approaches the value unity less rapidly than the experimental 
curve for the circular pipe, presumably because the regions where 
the boundary layers are thick are smaller in circular than in square 


pipes. 
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6. PRACTICAL IMPLICATIONS 


It has been shown that distortion of the velocity profile can cause the 
sensitivity of a flowmeter of circular section to fall from the Kolin value 
by about 8°,,, provided the fluid traverses the magnetic field for a distance 
at least souk to Ra/M before encountering the electrodes, the theoretical 
distance being somewhat larger than the observed ones. It is important 
to consider whether the fall will occur commonly in practice. ‘Typical 
figures for a small flowmeter of 1 in. bore bearing liquid sodium at 300° C 
at 100 cm/sec in a field of 3000 oersted are R = 3-3 « 104, M = 500 and 
R M = 66. Since the tield would probably extend for 2 or 3 rather than 
33 diameters upstream of the electrodes, we see that only at much lower 
flow rates would appreciable distortion of the velocity profile and fall in 
sensitivity occur. Variation of sensitivity due to upstream disturbances 
or edge effects is likely to be much more serious (Shercliff 1954, 1955). 

Further experiments to test the variation of S with wall conductivity 
would be interesting, but also very difficult owing to the importance of 
contact resistance. 


The experiments described above were carried out with the help of 
facilities provided at Cambridge University by the Atomic Energy Research 
Establishment, Harwell. 
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REVIEWS 


Engineering Fluid Mechanics, by C. Jarcer. (‘lranslated from the 
German by P. O. Wolf.) London: Blackie, 1956. 529 pp. 60s. 


In English-speaking countries the term hydraulics has long been under 
a cloud. ‘The disparagement from which it still suffers arose out of the 
difficulties of the subject that until relatively recently permitted investi- 
gators to give only their observations of depth, pressure and discharge. 
No analysis of the details of the motion or effective comparison with previous 
papers was possible, and all too often the records of their painstaking work 
terminated with nothing more than ‘the author’s formula’, applicable 
only to the exact conditions of the experiment. It is presumably for this 
reason that a misleading title heads this revised and enlarged translation 
of Dr Jaeger’s Technische Hydraulik, of which a French version has already 
appeared under the title Hydraulique Technique. ‘The book is concerned 
wholly with the motion of water. After a preliminary account of the 
physical basis of hydraulics and the fundamental equations the author 
turns to his three principal topics, steady and unsteady flow in open 
channels, unsteady flow in pipes due to surge and water hammer, and 
seepage, which between them occupy nearly three quarters of the book. 

For several reasons the practical problem of steady flow in a uniform 
open channel of uniform slope is far more extensive than that of steady flow 
in a circular pipe, and our knowledge of it is relatively scanty. ‘The channel 
may be of rectangular, trapezoidal, circular or egg-shaped cross-section 
with almost unlimited variation in the proportions of the cross-section, 
and the free surface destroys any possible symmetry in the velocity distri- 
bution over the cross-section. ‘I'wo types of flow can occur according as 
the velocity in the main stream is greater or less than that of a long wave 
moving under gravity on the surface. ‘They may be given the descriptions 
‘superundal’ and ‘ subundal’, which seem clearer than the adjectives 
‘shooting’ and ‘streaming’ which are more commonly used. They 
correspond to Froude numbers greater and less than unity, and an analogy 
exists with the flow of a gas which likewise is of two types according as the 
Mach number is greater or less than unity. In practice, friction effects 
are of vital importance because they determine the discharge through 
a channel of given shape and slope, but for a large canal the slope is so 
minute that it cannot easily be measured in existing works or correctly 
constructed in new works. Further complications ensue if the slope or 
cross-section of the channel is not uniform or if blockages, such as weirs 
are inserted. Unsteady flow may be due to the passage upstream of 
‘jumps’ of various kinds (such as may be seen in nature as bores in 
esturaries) or to the passage downstream of disturbances caused by floods 
or by the removal of obstructions. 

Unsteady flow in pipe-lines is of two kinds, dependent on the rapidity 
with which the control valve is moved. If the valve is shut suddenly, 
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a loud noise (water hammer) is heard which is familiar to many house- 
holders, and a pressure wave passes upstream, the intensity of which in 
pounds per square inch is roughly 60 times the original water velocity 
measured in feet per second. Now the velocity of sound in water in an 
almost rigid pipe is about 4500 ft/sec., and the original water velocity 
cannot be more than a very small fraction of this. Hence the aero- 
dynamicist would describe the shock wave as weak; not so the hydraulic 
engineer as he surveys his shattered pipe and damaged power station. 
Sut normally the valve admitting water from the pipe-line to the turbine 
is moved relatively slowly by the governor, which can also open a bypass 
for a short time. ‘lhe mass of water in the pipe-line is, however, so great 
that, to smooth out the resulting disturbances without excessive waste 
through the bypass, a surge tank is usually connected to the pipe-line. 
‘This is an.open vertical pipe, much larger in cross-section than the pipe- 
line and placed as close to the power station as circumstances permit. 
Under steady conditions the water level in the surge tank stands lower 
than that in the reservoir by an amount equal to the friction loss of head 
in the intervening pipe-line. ‘The installation can be regarded as a gigantic 
U-tube; one limb (the reservoir) is virtually of infinite cross-section, and 
from the base of the other (the surge tank) water is drawn off tothe turbine. 
Any change in this demand leads to damped oscillations of level in the 
surge tank which are governed by non-linear differential equations. If the 
demand increases, the surge tank serves as a source close at hand, and 
initially the level in it falls; if the demand decreases, part of the flow from 
the reservoir passes into the surge tank, the level in which begins to rise. 

Dr Jaeger deals with these subjects in great detail, and almost every 
possible case and sub-case is fully examined. ‘The treatment is entirely 
analytical supplemented, when the working grows too complex, by 
graphical methods. ‘he standard of mathematics that is employed is 
not high, but the algebra is often long and tedious. ‘The author evidently 
believes that there is no satisfactory means of escape from these laborious 
calculations, for alternative methods are cursorily dismissed. ‘Thus in the 
description of surge tanks six lines and five references suffice for an account 
of the importance and use of models and electrical analogues, as “ there 
is not, however, sufficient space to include a description of these techniques 
here’. (Additional references on the first of these matters have strayed 
into the chapter on water hammer.) 

The book does not contain examples to be worked out, and it is clearly 
meant, not for the general student, but for the designer and the advanced 
specialist. ‘lhe designer will not, however, find here all that he wants. 
For many of the problems alternative methods of solution are provided, 
but he is given little help in choosing which is the best; the treatment is 
comprehensive rather than critical. Moreover, he is made aware that the 
theories are inevitably based on a variety of assumptions, yet they are left 


unsupported by experimental evidence. ‘Thus a doubt may arise in his 
mind whether the lengthy computations are worth while—a doubt that will 
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not be lessened by the author’s remark that “it is not within the province 
of this book to investigate the usefulness and practical results of the methods 
of analysis developed in” the chapters on surges and water hammey. 
It is to be hoped that Dr Jaeger will soon provide this much wanted 
assessment. He remarks with great truth that observations of surge and 
water hammer are, regrettably, too rarely published, but in his list of 
references there is sufficient information for at least some comparison to be 
made. ‘This shortage of experimental evidence, which is seriously delaying 
the advance and spread of knowledge, seems to arise from financial pressure. 
A costly pipe-line must be put into commercial service the moment it has 
passed the simple tests required by the contract, and even a single day 
devoted to experiments would involve a serious loss of revenue. No lead 
in this matter has been given in the highest quarters, for in the course of 
a discussion on a paper on surge tanks read in 1933 before the Institution 
of Civil Engineers, a distinguished speaker expressed his hope of publishing 
surge observations on the Lochaber installation. ‘This hope is long 
deferred, for these readings have not yet been communicated to the 
institution. 

The book ends with appendices on empirical coefficients, flow over a 
movable bed, and density currents. A very valuable feature of the entire 
work is the enormous number of references extending up to 1955, which are 
presented in the form of footnotes. ‘These are accompanied by brief 
biographical notices of the more eminent authors. 

In his two prefaces Dr Jaeger describes his work as a handbook as well 
as a summary of his lectures at Ziirich, and in a very large handbook of 
this sort it is inevitable that the author’s chief interests are shown up by 
unevenness in the standard of treatment of the various topics. ‘Thus, in 
the chapter on the basic equations when he is dealing with the loss of 
energy at a sudden expansion in a pipe, he is content to reproduce the usual 
dubious analysis that survives in text books seemingly because it leads to a 
tidy result. It depends on the assumption that the pressure on the 
annular surface at the enlargement is the same as that in the smaller pipe. 
This assumption, he writes, ‘‘is supported by experimental evidence’’, 
but the references provided are scarcely convincing. ‘laylor’s name 
might have been included in Chapter II] amongst those who have 
gone deeply into that intractable subject, the theory of turbulence; and 
much more experimental confirmation of boundary-layer theory exists 
than can be found in Escande’s 1938 book on dams, which is the only 
reference provided on this point. Little use is made of Froude number to 
simplify the analysis of flow in open channels, and Chapter IV on steady 
flow in channels ends in confusion as the calculation of the critical depth 
for various shapes of cross-section is apparently extended to include free 
and aerated jets. ‘lhe detailed treatment of critical flow theory, which 
includes an interesting history of the problem, would have been improved 
by an account of the simple ideas introduced long ago by Hugoniot. ‘The 
wretched weir formula involving a dimensional coefficient of discharge, 
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which engineers are accustomed to use, is stated without protest; and the 
opportunity is missed of explaining the advantage of comparing the actual 
discharge with that over the theoretical broad-crested weir having the 
same width and head. Lemoine’s sinusoidal theory of undular jumps is 
noticed, but not its demolition by Benjamin and Lighthill. ‘The important 
advance made by Southwell and Vaisey in their application of relaxation 
methods to free-surface problems is mentioned only in a footnote to a 
remark on a problem which they did not treat. 

‘There is much to be said for the widely held view that the days are past 
when an author could produce a lengthy and balanced book covering the 
latest advances in many fields. It can no longer be expected that one man, 
no matter how scholarly and industrious he may be, can have everything 
at his fingers’ ends at the same time. Of recent years the stream of original 
papers has greatly increased in width, depth and velocity, and the only 
way of control is by division of the flood. ‘wo methods have been tried 
with success. ‘The first is the large book covering a number of topics, 
each with a chapter to itself written by a different yet expert hand. 
Conspicuous examples are Engineering Hydraulics and Modern Develop- 
ments in Fluid Dynamics under the respective editorship of Hunter Rouse 
and Goldstein, both of which have served as powerful springboards for 
further advances. A drawback of this method is the long delay that a 
big book encounters at the printer’s hands. Moreover, the editor, if 
conscientious, must expect to be brought by his labours almost to the 
grave. ‘lhe second method is the monograph of limited scope compiled 
and seen through the press by one man. ‘The final stages are relatively 
brief, and the book need not be seriously out of date as soon as it is 
published. ‘There is nothing to prevent a gifted author from producing a 
succession of monographs. In both methods it is unnecessary for the 
size of the book and therefore its price to be inflated by the inclusion of 
preliminary matter than can readily be found elsewhere. 


A. M. BINNIE 


CORRIGENDUM 


Equation (41) of my paper, “The influence of radiative transfer on 
cellular convection ”’ (7. Fluid Mech. 1, 1956, 424), should read as 


eee, (41) 


R. M. Goopy 
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